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An experiment is described which determines the variation in the number of extensive 
cosmic-ray showers per unit time with altitude from sea level up to 4300 m. The number of 
showers which should have been observed with the apparatus used is calculated on the basis of 
the cascade theory in which various power law energy distributions for the primary cosmic rays 
are assumed. Good agreement is obtained between the observations and the calculation at 
high altitudes for suitable choice of parameters; at low altitudes the observed excess is what 
would be expected from a mesotron component in the extensive showers. It is found that the 
same choice of parameters which will give good agreement with the extensive shower data will 
also describe cosmic-ray observations at much lower energies. It is pointed out that this offers 
strong support for the hypothesis that there is but one type of primary cosmic-ray particle. The 
effects of choosing the proton as this primary particle are discussed. 


INTRODUCTION 


awe prodigious energy associated with ex- 
tensive cosmic-ray showers is one of their 
most thought-provoking aspects. Auger and his 
co-workers! have estimated this energy to be 
above 10" ev. The entire energy of such a shower 
has been attributed to a single primary particle 
incident at the top of the atmosphere. The 
cascade theory?~* accounts for the development 
from a high energy primary of the very large 
number of ionizing particles observed in these 
showers. On the basis of scattering in the atmos- 
phere, Euler and Wergeland® have derived ex- 
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pressions which determine the spatial distribution 
of the shower particles about the axis determined 
by the direction of the incident primary particle. 
A fairly complete quantitative description is thus 
given of the number and spatial distribution of 
the ionizing particles in the shower at different 
depths in the atmosphere. Obviously, if the 
cascade theory is correct, the primary particles 
with energies greater than 10 ev must acquire 
their energies from some physical process hitherto 
unknown in kind or at least in order of magnitude. 
It seems important, therefore, to determine how 
well the cascade theory does describe the behavior 
of these extensive showers, and to determine the 
energy distribution of the incident primary 
particles as quantitatively as possible. 

The simplest attack on the problem is a 
determination of the variation in the number of 
extensive cosmic-ray showers as a function of 
altitude. If the cascade theory is adequate, the 
energy requisite for a primary incident on the top 
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Fic. 1. Counter arrangement under the roof of the station wagon for the measurement of extensive 
cosmic-ray showers. 


of the atmosphere to produce an observable 
shower at any given depth in the atmosphere is a 
sensitive function of that depth. Consequently, 
the number of extensive showers per unit time at 
any depth is a measure of the number of primary 
particles incident on top of the atmosphere with 
energies above a certain “cut-off” value. This 
latter is determined by the apparatus and by its 
depth in the atmosphere. Thus the variation of 
the number of showers with altitude should give 
definite information not only concerning the 
validity of the cascade theory but also a fairly 
quantitative determination of the energy distri- 
bution of the incident primary particles as well. 

Since the number of showers observable with a 
given apparatus increases rapidly with decreasing 
depth in the atmosphere, it is important that the 
observations be extended to as high an altitude as 
possible. With this in mind, an experiment was so 
designed that observations could be made on the 
number of extensive showers over a range of 
altitudes from sea level at Chicago to 4300 meters 
at Mt. Evans, Colorado. 


EXPERIMENT 


To measure the number of these showers per 
unit time at various altitudes, a set of four argon 
and petroleum-ether filled Geiger-Miiller tubes 
were mounted permanently just beneath the light 
wood and fabric roof of a Ford station wagon. All 
readings were taken in the open to eliminate any 
possibility of scattering from material above the 
counter set. Each counter was placed in an 
individual sheet-metal shielding can together 
with the first amplifier tube of the coincidence 
circuit. A one-half inch layer of Celotex was 
placed around each counter tube unit to furnish 
thermal insulation. One counter tube was 
mounted at the extreme forward end of the car, 
one at the extreme rear and the other two in 
vertical II-fold (twofold) coincidence midway 


between the extreme counters. The distance 
between the extreme counters was 2.5 meters 
(Fig. 1). Each counter had a diameter of 4 cm 
and an active length of 48 cm, giving an active 
area of 196 cm?. 

The counters were connected in fourfold 
coincidence by a conventional coincidence circuit 
and the usual precautions were taken to maintain 
constant working potentials on the G-M tubes 
and constant plate and grid potentials on the 
amplifier itself. 

The counter tubes were mounted permanently 

in the station wagon to guarantee constant 
geometry throughout the experiment and thus 
facilitate alternation in the measurements at the 
various altitudes. The observational procedure 
was to take a run at one altitude for a period of 
from twelve to twenty-four hours; next, to take a 
run at a second station for a similar period; 
finally, to return to the first station for a final run 
to complete the set. This continual alternation of 
the observations at different altitudes with a 
fixed apparatus tends to eliminate errors due to 
slow variations in the cosmic-ray intensity or in 
the circuit itself. Actually, the latter effect was 
extremely small for the observations made in 
Chicago at the start of the expedition differed 
from those made after the return by less than the 
standard statistical deviation for the entire 
Chicago sct. 
- The average depth below the top of the 
atmosphere for each observation station was 
determined by taking barometric readings for 
each run and averaging over the total counting 
interval. 

The accidental counting rate was determined 
from the onefold counting rates for each G-M 
tube as measured with a scale 16 scaling circuit, 
from the twofold counting rate for the central 
pair and the twofold counting rate for the 
extreme G-M pair. The twofold rate for the 
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extreme pair was assumed to be entirely acci- 
dental, and this together with the onefold rate 
for each extreme counter gave the observed 
resolving time for the circuit. This procedure 
gives a maximum value and in all cases the value 
of the resolving time so measured was less than 
10 ' sec. This resolving time was then used to 
compute the accidental counting rate for the 
present counter arrangement, which was treated 
as an effective threefold coincidence set with the 
onefold counting rates for the extreme counters 
and the vertical twofold rate of the central pair of 
counters as the onefold rate for a single effective 
central counter. In all cases the true fourfold 
accidental rate was completely negligible. The 
accidental rates so determined were, of course, 
also maximum values so that the accidental rate 
at Mt. Evans was less than the 0.1 count per hr. 
and that at Chicago was less than the 0.01 count 
per hr. computed for each respectively. 

Table I gives the observed counting rates with 
their standard statistical deviations obtained by 
the above procedure. 


COMPARISON WITH THE CASCADE THEORY AND 
DETERMINATION OF THE 
ENERGY DISTRIBUTION 


In order to interpret the experimental curve in 
terms of the validity of the cascade theory and to 
derive from it an expression for the energy 
distribution of the incident primary particles, it 
is necessary to carry through a complete calcula- 
tion of the cascade process for all showers which 
can be detected by the counter system used. 
Since the validity of any conclusions based upon 
such a calculation depend upon the way in which 
the calculation is carried out, the method 
followed in the computation will be outlined in 
detail. 

The computation is separable into three major 
parts: (1) the cascade computation itself for the 
development of the showers arising from a single 
incident primary particle of energy E, (2) the 
introduction of an energy distribution which will 
give the variation of counting rate with depth for 
vertically incident primary particles, and (3) the 
zenith angle calculation which then furnishes a 
quantitative description of the observed decrease 
of counting rate with depth in the atmosphere for 
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the 


parameters in 
arbitrarily chosen energy distribution expression. 


a suitable choice of the 


1. Cascade computation 

(a) The first step in the cascade computation 
for a single incident ionizing particle is the 
solution of the cascade equations to determine 
the total number of ionizing particles, m.,~ at a 
given depth z in the atmosphere, for a primary 
particle with a given incident energy /. For this 
purpose the cascade equations as given by 
Serber® and Snyder! were used. These equations 
express the number of ionizing particles as a 
function of the depth s below the surface of the 
atmosphere, the logarithm ¢ of the ratio of the 
incident energy to the critical energy in air, and 
an arbitrary parameter y. Values of the parameter 
€ were used corresponding to energies of the 
incident primaries of 1X10" ev, 5X10" ev, 
1X10" ev, «++, 1X10" ev. The depth equation 
was then solved for each of these values of «€, to 
determine the values of the remaining parameter, 
y, corresponding to the depths 630 g/cm*, 700 
g 800 g/cm*, «++, and 1100 g/cm*. The 
equation for the number of ionizing particles, 
when solved for the above values of the parame- 
ters y and e¢, vields the desired total numbers of 
ionizing particles, m.~ at cach of the chosen 
depths due to incident primaries of each of the 
chosen energy values. 

(b) Determination of the spatial distribution of 
the n. ¢ particles about the axis of the shower in 
the plane of observation. The axis of the shower 
is the direction of incidence of the primary 
particle. For this calculation the scattering 
formulas given by Euler® were used. While they 
are strictly valid only in the region of maximum 
development of the shower, the approximation is 
sufficiently good for the depths and energies here 


TABLE I. Extensive shower counting rate as a function of 
altitude.’ 


Deptu BELow 


Top or At- COUNTING 
ALTITUDE MOSPHERE RATE 

PLACE METERS IN G/CM? COUNTS/HR. 
Mt. Evans 4320 628 24.1 +0.3 
Summit Lake 3900 658 20.5 +0.5 
Echo Lake 3100 709 14.0 +0.5 
Idaho Springs 2190 795 7.3 +0.8 
Denver 1610 854 5.0 +0.4 
1025 1.47+0.15 
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considered. The values of the “half-width” of the 
shower Euler’s expression for the radial 
distribution were determined for cach of the 
chosen depths and the expressions normalized to 
yield the proper values of m. ¢ when integrated 
over the entire plane of observation. These 
expressions then give the mean values of the 
surface densities of ionizing particles, p.,z,,, at 
any distance r from the axis of a shower having a 
total number of particles at the depth z of n.,¢. 
Of course the actual distribution is statistical 
about this mean value. 

(c) Probabilities of observing extensive showers. 
The probability 'P,, ¢,, that a single G-M counter 
tube of active area A will be struck by an 
ionizing particle when the counter tube is at any 
given distance 7 from the axis of the shower can 
be calculated readily from the mean _ particle 
density p.,z,,. In the present case, the counter 
arrangement was not symmetrical about a 
vertical axis through its center. Consequently, to 
determine the probability of registering a three- 
fold coincidence ''P.¢,, it was necessary to 
average over all orientations ® of the counter set 
in the particle field. To do this the distance from 
the center of the counter arrangement to the axis 
of the shower was kept fixed at some value of 7 
and the distances to cach counter then de- 
termined for orientations of 0°, 30°, 60° and 90° 
between the radius vector from the counter set to 
the shower axis and the axes of the G-M tubes. 
The threefold coincidence probability was then 
calculated as: 


for each orientation, and the mean of these values 
(determined by a Simpson's rule calculation of the 
area under the ¢, curve when necessary ) 
gave the mean value of |"! P, x, ,. This orientation 
average was carried out in all cases for values of 
rof 3,1, 2, 3,4 and 5 meters. For greater values 
of the radius the effect of the orientation was in 
general negligible and except in a few cases no 
appreciable error was introduced by computing 
The orientation calculation 
was made for larger values of r in those cases 
where any sensible error would be introduced by 
the omission. The values of ''P., were ex- 
tended in this fashion for values of r of 10, 15, 
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20, +++ out to 100 m in the case of the high 
energy primaries. 

This means, of course, that there is a definite 
probability that a shower whose axis falls a 
distance r from the center of the counter arrange- 
ment will register a threefold coincidence. There 
is thus an “area of incidence’ at the top of the 
atmosphere over which an incident particle of 
energy £ will be effective in producing coinci- 
dences, and this area increases with the energy of 
the incident primary. 

If, now, there are incident vertically on the top 
of the atmosphere Nx particles of energy E per 
cm? per sec., their contribution to the threefold 
counting rate of this counter set when incident 
at a distance between r and r+dr will be 


Ne ¢.-X2ardr 


and the total contribution to the threefold 
counting rate will be 


0 


Since the integrand is an empirical function, 
this integration was carried out by a Simpson's 
rule determination of the area under the curve of 
IP. »,,2ar as a function of r. These numerical 
integrations were extended out to such values of r 
that contributions for still greater values of the 
radius were negligible. For the low energies the 
calculation could be stopped at 5 m or 10 m, 
while for the high energy primaries it was neces- 
sary to carry the calculation out to 90 m. 

This process was repeated for each of the 
values of £ at each depth z. The families of curves 


of ™'R. » as a function of / for each z, and of 


MMR, » asa function of zs for each £, were plotted 
as a check for possible arithmetical errors. Also, 


from these families of curves, values of MIR! 
were interpolated for values of E intermediate 
between the values for which the computation 
had been carried out. The nature of the curves 
was such as to make this readily possible and 
sufficiently accurate. Thus the cascade computa- 


tion has yielded a set of values of IR’ » for each 
of the energies (1, 2, 3, 4, 5, 6.25, 7.50, 8.75) 
X10!? ev, (1, ---)10'* ev (1, ---)10" ev, 
(1, ---)K10"% ev, 110'® ev at each of the 
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chosen depths. The above values of the inter- 
mediate energies were chosen to facilitate a 
Simpson’s rule integration over the energy 
intervals from 1X ev—5X10" ev, 5X10" ev 
—1X10" ev, ---5X10" ev—1X10'* ev at the 
next stage of the computation. 


2. Vertical counting rate vs. depth curves 


For this it is necessary to make some assump- 
tion concerning the nature of the energy distri- 
bution of the incident primary particles. The 
usual form was adopted, that is 


dNi, No. ai dE, 


where dN¢,q represents the number of particles 
per sec. with energies between E and E+dE 
incident vertically on each cm?® of surface at 
z=0, the top of the atmosphere. No’ and @ are 
arbitrary parameters to be chosen to give the 
best possible quantitative description of the 
experimental observations. Various values have 
been used for a@ in the literature under various 
circumstances. In the present calculation, values 
of dE, E* were computed for all of the above 
energies, including the interpolated values, for 
each value of a. The values of a chosen were 2.5, 
2.6, -3.0. 

The total vertical extensive shower counting 
rate at each depth is now 


0 


Again the integrand is an empirical function and 
the integration was carried out by using a 
Simpson’s rule determination of the area under 


the curve of ™'R! eNoeE-* as a function of F. 
The calculation was extended to a value of E 
of 5X10'® ev which was sufficiently large to 
make the contributions from still higher energy 
groups negligible for the present counter arrange- 
ment. 


3. Zenith angle computation 


The calculation to this point gives a set of 
relationships between the vertical counting rate, 


IR’ and the depth z, for each of the chosen 
values of the exponent a. These relationships do 
not represent the variation of total counting rate 
with depth, however. The above relationships 


show that the variation of counting rate with 
depth is extremely rapid, so that a given change 
of zenith angle at small depths will produce a 
much smaller change of counting rate than the 
same change of zenith angle at a relatively great 
depth. Consequently, it is essential to make a 
computation of the zenith angle effect. This 
calculation is complicated for the present counter 
arrangement since the entire upper hemisphere of 
incidence is not open to threefold coincidences 
due to the fact that unless a single particle can 
produce a twofold coincidence in the central pair 
a true fourfold rather than a threefold coinci- 
dence is required in order to register a count. The 
geometry is simple, however, so that it is possible 
to determine the effective depths corresponding 
to any zenith angle across any chosen zone on the 
hemisphere of incidence, and the fraction of each 
such total zone which is open to threefold 
coincidences. For the present counter arrange- 
ment and depths, the contributions due to 
fourfold coincidences were negligible from the 
portion of the hemisphere of incidence closed to 
threefolds. 

At the depth corresponding to the summit of 
Mt. Evans, the hemisphere of incidence was 
broken into zones at zenith angles of 0°, 21°6’, 
30°, 36°58’, 42°57’, 48°19’ and 53°16’. The depth 
corresponding to this last zenith angle is 1050 
g ‘cm? and the contribution from greater depths 
is negligible. Each of the first two zones was 
completely open to the threefold arrangement. 
Sixty-seven percent of the third zone was open, 
as was 48 percent of the fourth, 38 percent of the 
fifth and 31 percent of the sixth. At the 1000 
g/cm? depth the contributions from zones beyond 
the second were negligible. 

The contribution to the total counting rate 
from each zone was now determined from the 
vertical counting rate curves and the effective 
values of z sec @ across the zone. Again this was 
done by a Simpson's rule computation of the 
mean value. This mean value was multiplied by 
the value of the solid angle for that part of the 
zone open to the threefold arrangement and this 
gave the contribution from that zone to the 
total counting rate. The sum of these contri- 
butions from the various zones then gave the 
total counting rate at that depth. This, of course, 
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Fic. 2. The variation in the number of extensive cosmic- 
ray showers per unit time with depth in the atmosphere. 
The curves give the calculated variation for various 
choices of the power law exponent. 


assumes that the primary particles are iso- 
tropically distributed in angle at the top of the 
atmosphere, which should be justified for the 
energies here involved (>10" ev). 

Six such sets of values of counting rate as a 
function of depth were determined, one for each 
value of @ used. Each contains a single unde- 
termined constant, No.2, which now represents 
the number of particles per cm? per sec. per unit 
solid angle incident at the top of the atmosphere. 
This constant was then adjusted so that the 
counting rate at Mt. Evans was the experi- 
mentally observed value of 24.1 counts per hr. 


These curves for "'R,.. as a function of z for the 
various values of the exponent are shown in 
Fig. 2. 

DiscussION OF RESULTS 


It will be seen from Fig. 2 that there is excellent 
agreement between the cascade computation and 
the experimental observations at the high alti- 
tudes, that is at Summit Lake and at Echo Lake, 
for values of the exponent lying between 2.6 and 
2.8. At the lower altitudes, particularly at Chi- 
cago, there is definite disagreement for all values 
of the exponent. As pointed out in a recent note,’ 
this discrepancy can be explained qualitatively 
on the basis of a mesotron component ac- 
companyving these extensive showers. With this 
in mind, the choice of the best value of the 
exponent must be made from the high altitude 
observations. Here the electron-photon com- 
ponent is most fully developed and consequently 
the relative contribution of the mesotron com- 
ponent to the total spatial distribution of ionizing 
particles in the shower is very small. The Echo 
Lake observation is at an altitude sufficiently 
high for the mesotron correction to be small and 
yet far enough below Mt. Evans for the differ- 
ences between the calculations for the different 
values of the exponent to become appreciable. 
Basing the choice, therefore, on this point, the 
best value of the exponent would appear to be 
2.70 or somewhat larger, although any value from 
2.6 to 2.8 is permissible. 

The selection of a particular value of the 
exponent based on the present data alone is not 
as well defined as could be desired. On the other 
hand, the value of the coefficient which must be 
associated with any given choice of the exponent 
is much more precisely determined. This arises 
from two facts; first, that the Mt. Evans obser- 
vation is the best determined of the experimental 
counting rates; and second, that the cascade 
computation itself is most completely valid at 
high altitudes where the showers are most fully 
developed. The relationship between the coeffi- 
cient and the exponent as determined by the 
present data and computations is 


logio( No, a) = 


X (a—2.5)(a—2.4) 2.55a53.0. 
7N. Hilberry, Phys. Rev. 59, 763 (1941). 
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If numerical agreement with the extensive shower 
observations at the top of Mt. Evans is to be 
obtained, this expression must be satisfied for 
any choice of the exponent. To obtain good 
agreement with the Echo Lake observations a 
is limited to the range between 2.6 and 2.8. 

Figure 3 shows in the full curve the graph of 
d( NE) dE over the energy range up to 5X10" ev 
as predicted by this energy distribution ex- 
pression with a equal to 2.75. It should be 
emphasized that this expression has been derived 
from observations involving primary particles 
chiefly in the energy range from 5X10" ev to 
5X10" ev. In this case, the energy distribution 
expression becomes: 


dN=5.4;X10PE 


The dotted curve in Fig. 3 is that published by 
Bowen, Millikan and Neher® and the blocks 
represent their experimental observations. The 
broken curve in Fig. 3, for energy values less than 
1.0X10'® ev, has been obtained in the same 
manner as that of Bowen, Millikan and Neher by 
adjusting it so that the areas under the curve are 
kept equal to those of the corresponding blocks. 
As will be seen, the present curve fits the experi- 
mental observations remarkably well. ob- 
taining it, it has been assumed that over the 
major part of the observed energy range the 
particles are isotropically distributed at the top of 
the atmosphere. 

Integration of the energy distribution ex- 
pression to determine the energy per cm* per sec. 
carried into the top of the atmosphere from all 
directions by particles in various energy bands 
gives a value of 9.7X10*° ev for those particles 
with energies above 1.7X10'® ev as compared 
with the Bowen, Millikan and Neher observation 
of 9.4 10° ev. In the band from 0.67 X 10!” ev to 
1.710! ev, the present expression gives 9.8 X 108 
ev as compared with the experimental observa- 
tion of 8.7X10° ev. Since the observed distri- 
bution has already dropped well below the power 
law in the lower part of this energy range, this 
agreement is also satisfactory. 

*The value of the coefficient, 5X10" given in the 
previous note’ was taken from a linear logarithmic graph 
of coefficient against exponent. In re-evaluating the coef- 


ficient-exponent relationship, the more precise formulation 


seemed justified. 
8’ Bowen, Millikan and Neher, Phys. Rev. 53, 855 (1938). 
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If 2.74 is used for a, the energy which should 
be carried in by particles with energies above 
1.7X10'° ev drops to 9.0108 ev but it then 
becomes very difficult to fit a smooth curve to the 
experimental observations below 1.010!" ev. 
On the other hand, the choice of 2.76 for a makes 
possible a smooth fit below 1X10" ev but the 
energy which particles above 1.7 X 10"'° ev should 
carry in then becomes 10.5 10° ev and this is 
definitely in poorer agreement with the observa- 
tions than the result given by the expression 
with @=2.75. 

It would seem, therefore, that a single power 
law expression is capable of describing the 
distribution of the number of primary cosmic-ray 
particles per unit energy range from energies of 
10'® ev down to energies of 10" ev. 

The question immediately arises as to why a 
deviation from this power law should set in at 
1X 10'* ev, since this is a very high energy value 
to be affected by the solar magnetic field.’ A 
Letter to the Editor by Carlson and Schein'® now 
in press suggests that a large amount of energy is 
carried away by neutrinos in the disintegration of 
the mesotrons. If some such hypothesis is correct, 
the present power law may be valid down to 
energies even lower than 1X10" ev. Thus if the 
total energy carried into the top of the atmos- 
phere at Omaha be computed from the present 


Fic. 3. The energy distribution curve in the region of 
10" ev. The full curve gives the present power law ex- 
pression. The dotted curve is the curve of Bowen, Millikan 
and Neher and the blocks represent their experimental 
observations. The broken curve is the present power law 
distribution adjusted below 10" ev to fit the experimental 
observations. 


* P. Epstein, Phys. Rev. 53, 862 (1938). 
!” |]. F. Carlson and M. Schein, Phys. Rev. 59, 840 (1941). 
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expression, the value is 3.5410° ev per cm* 
per sec. from all directions as compared with the 
experimental observation of Bowen, Millikan and 
Neher of 2.25 10° ev per cm? per sec. There is 
then a discrepancy of 37 percent between the 
energy carried into the atmosphere as calculated 
and that observed, as compared with the loss of 
approximately 40 percent given by Carlson and 
Schein for the neutrino loss in the same energy 
range. On this basis it would seem possible that 
the power law given here may hold down to an 
energy considerably below 10! ev. 

The validity of a single energy distribution 
expression over an energy range from 10!® ev to 
10'° ev or lower indicates strongly that there is 
but one predominant type of primary cosmic-ray 
particle. The recent experiments of Schein, Jesse 
and Wollan' indicate that these particles are 
protons. This suggestion has been made previ- 
ously by R. Gunn” and W. F. G. Swann. 
Others" !® have emphasized protons as responsi- 
ble for the penetrating component. It must now 
be determined whether or not primary protons 
obeying the energy distribution given here are 
capable of accounting for the soft component of 
cosmic rays, the observed mesotron energy 
distribution and the ‘extensive shower phe- 
nomena. The recent discussions by Swann" and 
by Carlson and Schein!® describe the way in 
which the soft component can arise from such 
incident protons, and the mean energies used in 
their calculations are in good agreement with 
those computed from the present expression. 

As pointed out in the previous note,’ this 
energy distribution is capable of accounting 
qualitatively, likewise, for the mesotron distri- 
bution observed by V. C. Wilson'® in his deep 
mine measurements. 

The choice of the proton as the primary cosmic- 
ray particle has also certain advantages’ from the 
point of view of the mesotron component in 
extensive showers. On the other hand, for the 


"Schein, Jesse and Wollan, Phys. Rev. 59, 615 (1941). 

1 R. Gunn, Terr. Mag. 38, 247 (1933). 

183 W. F. G. Swann, Rev. Mod. Phys. 11, 251-254 (1939); 
Phys. Rev. 58, 200 (1940). 

4 A. H. Compton and H. A. Bethe, Nature 134, 734 
(1934) 


(1939). 
16V.C. Wilson, Phys. Rev. 53, 337 (1938); Rev. Mod. 
Phys. 11, 230 (1939). 


6 T. H. Johnson, Rev. Mod. Phys. 10 (1938); 11, 208. 


analysis presented here to furnish a valid de- 
scription of the primary energy distribution, it is 
essential that practically the entire energy of a 
primary proton in the range above 10" ev be 
transferred to electrons and photons very near 
the top of the atmosphere. It is possible that the 
radiation cross section of the proton in this 
energy range may prove to be sufficiently large to 
account for this energy transfer directly. If, 
however, the hypothesis of the explosive pro- 
duction of mesotrons by protons!” proves to be 
justified, it would seem possible that it is this 
process which forms the first step in the genesis 
of an extensive shower. In the latter case, a large 
multiplication with a consequent large division of 
energy occurs in a single event practically at the 
top of the atmosphere. This would leave a group 
of extremely high energy mesotrons very near the 
top of the atmosphere as the second stage in the 
development of an extensive shower. 

The next phase in the growth of the shower 
would depend on the behavior of the mesotron at 
these extremely high energies. Little is known 
with certainty from the theoretical standpoint 
concerning the nature of the mesotron. Experi- 
mentally, however, some information can be 
deduced from the absorption curves of V. C. 
Wilson.'® If the decrease in intensity of the 
mesotrons with depth is due mainly to ionization 
loss as is almost certainly true, Wilson’s absorp- 
tion curve gives the integrated energy distri- 
bution for the very energetic mesotrons and, as 
pointed out before, this agrees well with the 
present primary proton distribution assuming a 
fairly constant production cross section at these 
high energies. There is a sudden change in this 
mesotron ‘energy distribution, however, at a 
depth of about 250 m water equivalent, corre- 
sponding to an energy in the neighborhood of 
10" ev. For greater depths, the change of 
intensity with depth is much more rapid. The 
actual exponents for Wilson’s absorption curve 
are —1.77 corresponding to energies less than 
10" ev and —2.52 for energies greater than this 
value. This marked break in the mesotron energy 
distribution indicates either that many fewer 
mesotrons are formed above this energy value or 
else that those formed in this range lose energy 
rapidly until they fall below this critical value. 
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There appears to be no reason why mesotrons 
should not be formed at high energies unless the 
radiation cross section of the proton in this 
energy range becomes large in comparison with 
the mesotron production cross section, which 
itself must be very great. In such a case the 
radiation of the mesotron would not be required 
as an intermediate process in the transfer of the 
energy from the proton to the shower particles. 
The proton itself would then be the shower- 
producing particle. 

If the high energy mesotrons are formed 
abundantly, radiation from those with energies 
over 10'' ev must be assumed. This is in accord 
with the conclusion arrived at by Christy and 
Kusaka" in their analysis of bursts. The value of 
the radiation cross section must become large 
above approximately 10" ev to account for the 
magnitude of the break in Wilson’s energy 
distribution. On this basis, the third step in the 
formation of an extensive shower would be the 
transfer of energy from the energetic mesotrons 
to photons by successive radiation until the 
mesotron energy has dropped below the 10" ev 
value. Since the minimum initial energy of the 
mesotrons in the present case would have been 
above 5X10" ev this means that all but one or 
two percent of the energy of the softest mesotrons 
entering into extensive shower production would 
go into photons and electrons. This would mean 
that even for the lowest energy protons which 
make any appreciable contribution to the 
counting rates observed in the present experi- 
ment, all but some ten percent of the energy of 
the initial proton eventually would go into 


R. F. Christy and S. Kusaka, Phys. Rev. 59, 414 
(1941). 


typically developing shower particles. For the 
majority of the protons, the energy carried off by 
the mesotrons would be very much smaller. 

It is probably true that a certain fraction of the 
photon energy is also converted back into slow 
mesotrons,'**" but this again should have but a 
relatively small effect upon the total energy 
appearing in shower particles at high altitudes, 
and consequently in the total number of ionizing 
particles reaching the plane of observation. 

It would thus seem quite possible for a primary 
proton to generate an extensive shower which 
would differ from an electron generated shower of 
the same initial energy at depths of sixteen or 
more cascade radiation units below the top of the 
atmosphere by less than the errors in the present 
observations and those inherent in cascade 
computations of the kind here carried out. 
Consequently, it would still appear justifiable to 
assume that the expression here given does repre- 
sent approximately the energy distribution of the 
primary protons in the energy range from 10'* ev 
to 10" ev. 

In conclusion, the author wishes to express his 
indebtedness to Professor A. H. Compton for his 
advice and encouragement throughout the pro- 
gress of this work, and to thank Professor Bruno 
Rossi for his advice on the arrangement of the 
experiment and for discussion of results during 
the expedition. To Professor Joyce Stearns of 
Denver University, whose unfailing thought- 
fulness and assistance play such a large part in 
the success of any Mt. Evans expedition, the 
author is deeply grateful. 
18M. Schein and V. C. Wilson, Rev. Mod. Phys. LI 
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The angular distribution of the alpha-particles ejected in the reaction, Li7+H'—2 He’, is dis- 
cussed. It is assumed that the Li?-nucleus is odd, that the incident proton is in a P state and 
that a broad and a sharp nuclear resonance level participate in the reaction. The angular 
momentum 0 and 2 h are assigned to the broad and to the sharp levels, respectively. The con- 
stants entering into the theory could be determined only by extending the measurements on 
angular dependence to proton energies higher than 400 kev. The results are in qualitative 
agreement with the dependence on energy of the reaction yield and of the angular distribution 


observed so far. 


INTRODUCTION 


ECENT experiments! have revealed a 

marked departure from spherical symmetry 
in the distribution of alpha-particles created by 
the transmutation 


Li?+H'—2Het. (1) 


In these experiments a thin film of lithium was 
bombarded by protons of energy ranging from 
100 to 400 electron kilovolts. The emitted alpha- 
particles have a total of about 17 Mev kinetic 
energy. The relative number of alpha-particles 
was determined at angles of eight different 
cosines and for eight different energies of the 
incident protons. 

In the center of mass coordinate system the 
relative number of alpha-particles per unit solid 
angle is represented very well by the factor 
1+A() cos*@; here @ is the angle between the 
direction of the emerging particles and the 
direction of the incident proton beam. The 
coefficient which gives the departure from 
spherical symmetry, A(/2), is positive and small 
at 100 kev, rises to a value of § at 300 kev and 
then rises rapidly to about } at 350 and 400 kev. 
A plot of A(/2) is reproduced in Fig. 1. 

It is greatly to be desired that the experiments 
be carried out for proton energies greater than 
400 kev. 


* National Research Fellow. . 
tV. J. Young, A. Ellett and G. J. Plain, Phys. Rev. 58, 


498 (1940). 


As was pointed out by those who did the 
experimental work the behavior of A(/2) suggests 
an approach toa resonance level of the compound, 
Be’, nucleus. In the present paper the expected 
theoretical behavior of the angular distribution 
of alpha-particles in the neighborhood of such a 
resonance is examined in detail. It is assumed 
that the alpha-particles which are ejected at 
energies outside the apparent resonance region 
are due to a second, but very much broader, 
resonant level of the compound nucleus. 

The energy-dependence of the yield curve* 
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Fic. 1. The variation with bombarding energy of A() 
in the factor, 1+A(£) cos@, which represents the experi- 
mental results on the angular distribution of alpha- 
particles. 


2L. H. Rumbaugh, R. B. Roberts and L. R. Hafstad, 
Phys. Rev. 54, 657 (1938). 
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Fic. 2. The ratio of experimentally found yield to the 
calculated penetration factor as a function of the bom- 
barding energy. Experimental results for energies above 
350 kev were taken from reference (2), those between 150 
kev and 350 kev from the paper of Herb, Parkinson and 
Kerst, Phys. Rev. 48, 118 (1935), and those below 150 
kev from the letter of Heydenburg, Zahn and King, Phys. 
Rev. 49, 100 (1936). These three sets of results were 
joined to form a smooth curve on an arbitrary scale. The 
penetration factor used here is F?/p in the notation of 
reference (3). 


measured for alpha-particles ejected at right 
angles to the incident protons corroborates the 
presence of two processes giving rise to the 
emission of the alphas. Figure 2 shows the ratio 
of the yield and the penetration factor for 
protons (in arbitrary units). The penetration 
factor was calculated from the tables of Yost, 
Wheeler and Breit* by assuming a P wave and an 
effective collision radius of 4107" cm (the last 
two assumptions are of minor importance in 
determining the shape of the curve). The figure 
shows the presence of a process that depends on 
the energy practically only through the pene- 
tration factor and a second process which has a 
resonance around 400 kev. As will appear later we 
shall associate the former process with a J=0 
compound state, the latter with a J=2 state. The 
parallelism of the yield curve with the energy 
dependence of the angular distribution is striking. 
In view of our conclusions to be given below it 
would be more satisfactory if the maximum in 
Fig. 2 occurred at higher energy than _ the 
maximum deviation from spherical symmetry. 
But too much reliance must not be placed upon 
conclusions based on the ratio in Fig. 2. The 
resonance in this figure extends over so wide an 
energy range that the slowly varying influence of 


3F.L. Yost, J. A. Wheeler and G. Breit, Terr. Mag. 443 
(December, 1935). 


the energy on the yield may cause a difference 
between the position of the apparent maximum 
and the position of the resonance. 

The angular distribution of ejected alpha- 
particles will thus depend upon the simultaneous 
influence of two nuclear levels. Under this con- 
dition the relative phase of the particles coming 
from the two resonance levels will be of im- 
portance. In most of the considerations involving 
dispersion theory, intensities alone were im- 
portant; but in at least one example, the scatter- 
ing of neutrons by He, phase relations had to be 
taken into account.‘ In the subsequent discussion 
the superposition of outgoing alpha-rays with 
different phases will lead in general to a sigmoid 
dependence of the quantity A() (from the 
factor 1+A(£) cos*@) on the energy of incident 
protons. This sigmoid dependence is not unlike 
the dispersion curve in optics. The shape of both 
curves is due to essentially similar resonance 
denominators in the dispersion formula and the 
similarity of the curves helps to emphasize that 
the dispersion theory of nuclear reactions and the 
dispersion theory in optics have more in common 
than a name. 


STATISTICS AND PARITY 


The disintegration (1) in the energy range 
involved permits of a fairly simple theoretical 
discussion. Alpha-particles obey Einstein-Bose 
statistics so that a pair of them can be produced 
only in states of even orbital angular momentum 
and even parity ; since the spin is zero the allowed 
states may be designated 'D,'G- ++. Further- 
more the ground state of Li? is well known to have 
a nuclear spin of 3, whereas the parity of this 
state has not been experimentally determined. 
Theoretical predictions both from the individual 
particle model and the alpha-particle model are 
that the ground state of the Li? nucleus is odd. In 
this case only P, F, --- states of the incident 
proton will be able to form the even compound 
nucleus Be*. Computing the intensities of P and 
F “regular”’ solutions to the Schroedinger equa- 
tion in a Coulomb field’ it is found that the P 
wave is over 1000 times stronger than the F wave 
for particles of 200 kev. Even up to 1 Mev the 
F wave is inappreciable so if the parity of the 
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Li? nucleus is odd only the P wave need be taken 
into account. The relatively small yield of reac- 
tion (1) also points to the process being due to an 
incident proton in a P wave rather than in an 
S wave.® 

If it is assumed that the reaction (1) is pro- 
duced by protons in a P wave the experimental 
angular distribution, 1+A(£) cos*@, can be 
readily understood. Since the bombarded nuclei 
are completely unpolarized and since the spins of 
the incident protons are also unpolarized the 
emitted beam cannot show more complicated 
transformation properties under rotation of 
coordinates than those of the orbital motion in 
the incident beam. Thus the emitted beam may 
contain a part which behaves as a vector and 
gives a cos*@ term but no higher powers of cosé@ if 
the incident wave function transforms like a 
vector, i.e., if it is a P wave. The Bose statistics 
of the alpha-particles excludes odd powers of cos@ 
so that in fact 1+A(£) cos*@ is the most general 
angular dependence to which incident P waves 
can give rise. 

Although there are many indications that Li? 
nuclei are odd the possibility that they are even 
is not to be overlooked. In this case the incident 
proton must have an even wave function and an 
even angular momentum. In the angular distri- 
bution of the alpha-particles the term A (£) cos* 
might be explained as due mainly to an inter- 
ference term produced by an incident S wave and 
a less effective incident D wave. The square of 
the amplitude of the outgoing alpha-particle 
wave function produced by the incoming D wave 
contains in general a term proportional to cos*@. 
But if the incoming D wave is relatively ineffective 
in producing the reaction, the cos‘@ term may 
have too small a coefficient to be observed experi- 
mentally. The cos*@ term depending linearly on 
the small amplitudes produced by the incident D 
wave would in fact be much more noticeable than 
the cos‘@ term depending quadratically on the 
same amplitude. On the basis of present experi- 
mental evidence it is not quite impossible that Li? 
nuclei are even and the 1+A(£) cos*@ angular 
dependence is produced by joint action of inci- 
dent S and D waves. But in principle it is possible 
to establish the odd parity of the Li’ nucleus by 
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eliminating experimentally in the reaction prod- 
ucts the cos'@ terms to a high accuracy thus 
demonstrating that the angular dependence is 
not produced by an incident D wave. 

It will be assumed in the following calculations 
that the Li’ nucleus is odd and it will be shown 
that the experimental data obtained so far may 
be understood with this assumption. The incident 
proton P wave has a total angular momentum of 
3h at most and this combined with the 3h of the 
nucleus can give rise to all integral angular 
momenta from 0 to 34 for the compound nucleus. 
Outgoing alphas are therefore either in 'S or 1D 
states. From the shape of the experimental curve 
for A(E) one would judge that a resonance in one 
or the other of these states is approached as the 
proton energy increases. Below the apparent 
resonance, eé.g., at 100 kev, the angular distri- 
bution approaches spherical symmetry. Barring 
very special types of interaction which would 
make the angular distribution due to a compound 
nucleus of J=2 spherically symmetric it appears 
that the alpha-particles of low bombarding 
energies are ejected from a compound nucleus of 
J=0. 


DISPERSION FORMULA 


In view of the probable applicability of the 
simple conditions discussed above we have calcu- 
lated the expected angular distribution of alpha- 
particles in reaction (1). For this purpose the 
dispersion formula of wave mechanics as pre- 
sented by Breit and Wigner® and Bethe and 
Placzek’ has been used. According to the dis- 
persion theory as applied to the calculation of the 
cross section o”9 of the nuclear transformation 

|? 


Ep—E,+317,| 


a? 


In Eq. (2) 7”, is the matrix element of the 
Hamiltonian which couples the initial state P 
with the compound (resonant) state 7; [7g is the 
same thing for the final state Q. Both P and Q 
may be degenerate, and there are many levels r. 
Ep is the energy of the initial state, Z, the energy 
of the intermediate state and y, is the level 


6G. Breit and E. Wigner, Phys. Rev. 49, 519 (1936). 
7H. A. Bethe and G. Placzek, Phys. Rev. 51, 450 (1937). 
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width.’ In calculating y, one must take into 
account all possible modes of decay of the state r. 
In the present problem the most probable 
method is just the emission of two alpha- 
particles and in this case 


Le! (3) 


where the summation extends over the degenerate 
states Q. 

We are interested in intermediate states, 7, of 
two different angular momenta, 0 and 2h, which 
have resonances in the region of the experimental 
results and which are eigenstates of the angular 
momentum. The initial state P is a Li’ nucleus 
with a total angular momentum 3A plus an 
incident proton with an orbital angular mo- 
mentum / and a spin momentum 3h. Thus the 
initial state does not appear as an eigenstate of 
the total angular momentum. By proper linear 
combinations of the initial states with various 
orientations of angular momentum (including the 
three possible orientations of the P wave) one can 
construct eigenstates of the total angular mo- 
mentum which we shall call X,. The subscript r 
indicates that X, has the same transformation 
properties with regard to rotation as the inter- 
mediate state r. It will be seen that two functions 
can be constructed with angular momentum 


J=2. We shall call X. that linear combination of: 


them to which the intermediate state r with J =2 
responds most strongly. The orthogonal linear 
combination is then not coupled at all with the 
compound nucleus. The construction of the state 
Xo (which is invariant under rotations) is unique. 

The matrix element 77’, in the numerator of 
Eq. (2) may then be written in the form 


H?,= ¢'(E)a’ (P, X,). (4) 


Here (P, X,) is the transformation coefficient 
from the state P to the state X,. The matrix 
element connecting X, and r has been written as 
a product of two factors: the energy independent 
complex number a’, and the energy dependent 
y’(E). We shall assume ¢’(E£) to be the same 
function for all intermediate states involved. 


§ As has been shown by G. Gamow, Atomic Nuclei and 
Nuclear Transformations (Cambridge, 1937), p. 96, the 
energy of the unstable state r may be considered to be the 
complex value E,—}iy,; the }iy, gives the exponential 
decrease of the amplitude in that state with time. 


This assumption is probably approximately valid 
over the relatively narrow energy region of the 
resonance. Actually the most important reason 
for the energy dependence of ¢’(E£) is the varia- 
tion of the penetration factor. This variation will 
not depend greatly on the nature of the inter- 
mediate state (J =0 or 2) over a narrow resonance 
region as long as the angular momentum of the 
initial wave (which is a P wave) is the same. 

The final state Q is represented by two alpha- 
particles the wave function of which can be re- 
garded as plane waves which describe the ejection 
of those alphas in a definite direction character- 
ized by the angle @. This final state will again 
not be an eigenstate of the angular momentum. 
We may introduce an eigenstate Z which can be 
obtained by superposition of the final states 
corresponding to various directions 6. The matrix 
element J/’g can then be written 


Q), (5) 


where the symbols have a similar significance as 
in Eq. (4). The transformation coefficient (Z, Q) 
is proportional to the surface harmonic (nor- 
malized Legendre polynomials) Y,(@). Here the 
index 7 refers to the intermediate state and is 
meant to specify both the angular momentum J 
and its component, m, around the direction of the 
incident protons. The surface harmonic Y,(@) can 
be assumed to be normalized since we still carry 
the numerical coefficient a”’,. The energy depend- 
ence of ¢”(E) is negligible in the region of 
resonance due to the high energy of the ejected 
alpha-particles. It follows from the same reason 
that according to Eq. (3) we may replace y, by an 
energy-independent quantity I,. The equation 
for angular distribution then becomes 
a(P, X,) Y,(8) |? 


1+A(£) cos*@~ > ‘ 
Ep—E,+ | 


Here a’.a”’, has been set equal to a, and the 
factors y’(£)¢’’() common to all terms on the 
right-hand side have been dropped. 


Spin-OrBIT COUPLING 


The incident proton is assumed to have orbital 
angular momentum L=1 for the reasons given in 
the introduction. There is certainly no com- 
ponent of this orbital angular momentum about 
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the direction of bombardment so the states XY, 
cannot be obtained as a direct superposition of 
the initial states. From the j= 3 of the Li? nucleus 
and the spin of the incident proton we may 
compose two states: a threefold degenerate state, 
P3, in which the proton spin is opposite to the Li 
spin and a fivefold degenerate state, Ps, in which 
the proton spin is parallel to the Li spin. The 
plane wave describing the proton appears as a 
factor in both P; and P;. It is evident that only 
P; can form a compound nucleus of J=0 and 
contribute to the resulting S wave. Both P; and 
P; are able to form compound states with J= 2. 
Thus two intermediate states with J=2 can be 
constructed. 

If there is no significant spin-orbit interaction, 
i.c., if the splitting caused is smaller than the 
level width it would seem most likely from 
experimental results on magnetic moment that 
the ground state of Li? nuclei can be described as 
*P°,,;. In this case the spin of the odd nuclear 
particle is parallel to the orbital angular mo- 
mentum and then P; alone can contribute to the 
'D wave of alpha-particles; for the spin in the P; 
is parallel to that in *P°®,, and in the absence of 
spin-orbit interaction the compound nucleus 
would remain in the triplet state. The two 
resonant energy levels with J =2 coincide when 
there is no spin-orbit interaction but only Ps; 
gives rise to alpha-particles. 

We shall define strong spin-orbit interaction to 
be strong enough to separate the two levels with 
J=2 by an amount large compared with the 
level width. Resonance with one of these levels 
will be considered as important in the cnergy 
range around 400 kev. The more general case in 
which the levels are separated only by an energy 
comparable to their widths will not be presented. 
For the problem at hand the neglected level may 
be at higher energy than the one revealed so far 
in the experiments or it may be at negative 
bombarding energy and thus unattainable in the 
reaction. Further consideration will be limited to 
resonance with two compound states, one with 
J=0 which will be designated 79 and one with 
J=2 which will be called re. 

As has been stated above, the resonant state r2 
responds most strongly to the linear combination 
X>». To any orthogonal linear combination 72 will 
not respond at all. In the next section we shall 
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construct two states with J=2 which we shall 
call and X,". The state is constructed 
to be orthogonal to P; and X2 orthogonal to P3. 
We may express X»2 as a linear combination of 
these two states 

sing, 
where é is a real number.” The state 
sinE— Xe cosé 


will not be coupled with 72. It may be noted that 
the case of no spin-orbit coupling is represented 
by £=0. In this case a nuclear state interacting 
with XY, is present in the resonance region but 
this state does not disintegrate into two alpha- 
particles. For strong spin-orbit coupling —=0 is 
also possible but in this case a nuclear state 
interacting with X,, though it can disintegrate 
into two alphas, falls into a region sufficiently far 
removed from the one under consideration. 

Since P; alone is coupled with the nuclear 
state with J=0 one may substitute (P3, Yo) for 
(P, Xo) in Eq. (6). Similarly sing 
is alone responsible for the excitation of the 
nuclear state with J=2. Substituting this ex- 
pression in those terms of Eq. (6) which contain 
the second surface harmonic we obtain 


1+A(/) cos*é 


2 |a2(P;, X¥2°") ¥2"(0)|? 
m=—2] Ep— | 


~sin*é 


+cosé i. (7) 


where the magnetic quantum number 7 is made 
explicit. Since waves of different m values and Ps; 
and P; have no phase relation, mixed products 
containing different m values and P; and Ps are 
omitted. Terms (P3, Ys) and (P;, X2) have 
been omitted because of the assumed orthogo- 
nality relations. 

® There is no phase relation between the initial states P; 
and P; so that interference terms between P; and P; do 
not occur. Thus only absolute values of the coefficients of 


X, and X,® matter and we can use the sine and cosine 
of a real number for these coefficients without loss of 


generality. 


| 
1 Xo) Yo(@) 
Bp— Eo t 
( 


ANGULAR DISTRIBUTION OF ALPHA-PARTICLES 15 


all TABLE |. Transformation coefficients. s axis so that the only component of b which 
N29) | (Ps, Xe ANGUL: PENDENCE 
Incident waves P; and P; are then described as 
of 2 0 —1\/6 | (15/8)le?'* sinte follows: 
1 3,/2 6 (15/2)te'* sin@ cosé 
0 6 0 (5/4)1(1 —3 cos?@) Lays; }bo, 
by 2 | (15/2)'e~'* cosé 1 
CALCULATION OF MATRIX ELEMENTS v3a_ 158; —a_4s_, ]bo, 
” Let “a” represent the normalized wave func- 14540, 
d tion of the nucleus. It is fourfold degenerate IWv3aysy+as_) Joo, 
having components @;, @—;, Similarly let 
bo, describe the normalized incident P 
it 1 Oy 1 [a_ysy+ays_; (P35) 
1. wave of the proton and s;, s_; the spin state of the \2 
is incident proton. The subscripts give the s com- ; 
ce ponent of the angular momentum of each state. 
e We take the direction of bombardment as the A_148—4Do. 
The functions VY», Ye and XY. can be constructed from the same waves, a, > and s: 
r 
—} —v3a_ +v2a_ bos; +a_yhys_; —ayb_ — V2aybos_y + V3ayb_ 18-4}, 
2\3 
], 
- 
| la \la — Ja yb (6) i], 
2(6)! 
[ (6) *a_ + 2a_,b_ vla_ bos 2a,yb_ ], 
[ ~Qva Pos; +a 
23 
43 
1 
2v2 
2\3 


One easily verifies that (P3, X2)=(Ps5, Xo) =(X2, Xe) =0. The transformation coefficients 
(P3, Xo) and (P5, X2) and the normalized angular dependence of the surface harmonics for the 
different values of magnetic quantum number are given in Table 1. (Ps, Xo) = v3. 
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Substituting the results of Table I into Eq. (7) and summing up the coefficients of ao”, aoa2 and a,” 
one obtains 


1| ay” 2(5 2)}anae(1—3 cos*@) cosé 
() 
3 (Ep— (Ep— }tlo)(Ep— 


(5 2)as*[(1+3 cos? @) cos?é+3 sin®é(1—cos*@) ] 
(Ep—Es+hil2)? 


As has been anticipated in the introduction only the zeroth and second power of cos@ appear in this 
formula. 


THEORETICAL ANGULAR DEPENDENCE 


It is assumed for the reasons given above that the sharp rise of A(/) with increasing energy is due 
to an approach to resonance with a state of J=2. The spherically symmetrical contribution below 
resonance then comes from a broad resonance to a state of J=0. In Eq. (8) we accordingly assume 
Py>>Ep— Ey and neglect the dependence of the J=0 resonance on energy. The unknown quantities 
a2, a, Py and IP; can be lumped into one complex number, Be’? which is conveniently defined as 


(9) 


8 and éare real. Since 31, is unknown it is expedient to choose it as the energy unit and instead of the 
energy difference Ep— we shall use 


Be’ = (10) 


(10) 


e=(Ep—E») 


With the notations defined in Eqs. (9) and (10) we get 


1+A(e) cos*@~ 1+ €?+28(1—3 cos*4)(€ cosd+sind) cosé 
+cos*&B?(1+3 (11) 


from which follows at once 


Ale =3 cos2~—28(€ cosé6+sind) cosé 
28(€ cosd-+sin6) sin*é) 


The denominator of the right-hand side of Eq. (12) is positive and vanishes only under the particular 
conditions: sinE=0, 8 sind cos§ = —1. This means that the emission in the direction of the beam does 
not vanish at any bombarding energy unless these conditions are fulfilled or unless A = —1.'° In 
general A(e) has one real root which occurs for the value €= 9, with 


(13) 


€) = —tané+38 secé cos2é secé. 


At the energy corresponding to ¢) the angular distribution is spherically symmetrical. In terms of €, 
A(e) becomes 


€ 


(12’) 


=68 cosé cost 
cosé cos§+ 26? 


” It is interesting to note that if A(e) falls below —1 in a certain energy range there will be no emission at those angles 
6 for which cos? @=1/|A|. Thus in this range there will be a cone of zero alpha-particle intensity the opening angle of 
which rises from zero to a certain maximum value and falls back to zero again as the bombarding energy is increased. 
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The plot of A(e) against ¢€ is clearly a sigmoid 
curve asymptotic to the eaxis at large absolute 
values of € and crossing the axis at ¢€). There is 
one maximum and one minimum to this curve 
and these lie at €» 


Em = €o+ +14+28?)!. (14) 


It is characteristic of this theory that the maxi- 
mum and the minimum occur at equal energy 
differences from ¢€ . Substituting Eq. (14) into 
Eq. (12’) we find the extreme values of A (e) 


A(€m) = —38 cosé cost/(€m+ 8 cosé cost). (15) 


It will be noted that in general the maximum 
values of A are not equal in absolute value. If, for 
instance, the energy corresponding to €) does not 
lie in the resonance region (|€)|>>1), one of the 
extreme values of A will differ from zero insig- 
nificantly. From Eqs. (15) and (10) 


3/A max 3/A min 
(Emax — Emin) /3T 28 cosé cosé, (16) 


where Eyax and in are the E values at which 
the maximum and minimum occur. Thus the 
measured values of A and E at the maximum and 
minimum determine }8IT, cosé cosé at once. The 
resonant energy 2 is then related to the energy 
at which A =0, E°, by 


+4128 cosé cost[2+3/Amax+3 Amin]. (17) 


Actually it is not at all certain that the shape 
of the A(E) curve has a sigmoid appearance. As 
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has been pointed out above in one-half of the 
sigmoid curve the value of A(Z) may be very 
small. In fact if the energy corresponding to €¢ 
should happen to lie at negative energies of the 
incident protons one-half of the sigmoid curve is 
missing. Thus in appearance, or even in fact, the 
dependence of A on the energy may reduce to a 
simple hump. The observed dependence of the 
yield of alpha-particles on energy makes it more 
likely that A as a function of E has a simple 
maximum rather than a sigmoid behavior. In 
fact if A(E) would become negative above 400 
kev and yet within the resonance region one 
would expect the maximum in Fig. 2 to lie at 
higher energies than 400 kev. This is so because a 
negative value of A(£) within the resonance 
region would enhance particularly strongly the 
number of alpha-particles ejected at right angles 
to the incident beam and the experiments on 
which Fig. 2 is based were made with just such a 
geometrical arrangement. 

Though in the present calculations explicit 
assumptions about the nuclear model have been 
avoided some simplifications had to be intro- 
duced. In particular we assumed that only two 
compound nuclei are involved in the process and 
we also neglected the energy dependence of the 
matrix elements involved within the resonance 
range. It seems however that it is not worth while 
to make more refined calculations at present. 
Extension of the measurements on angular 
dependence to higher energies will probably 
indicate to what extent the basis of the present 
calculations should be modified or extended. 
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It has been shown by Fermi that the cross section for the scattering of slow neutrons by 
protons is four times greater for protons strongly bound to an infinite mass than it is for free 
protons initially at rest. A simple generalization of this factor of four can be given for the 
scattering of neutrons by molecular gases if the neutron energy is great as compared to the 
energy differences of the rotational levels of the scattering molecule but small compared to the 
quanta of molecular vibration. In this case, the proton may be replaced by a freely moving 
hypothetical mass point whose mass is a tensor depending on the mass and structure of the 
molecule. On this basis the total cross section for the scattering of slow neutrons by Hz and by 
CH,, NH; and H,0 at low temperatures has been calculated as a function of the ratio of the . 
neutron energy tothe thermal energy of the scattering gas. The possibility of applying the theory 
to a determination of the neutron-proton scattering cross section is discussed. 


1. INTRODUCTION 


HE problem of the scattering of slow 

neutrons by bound protons was_ first 
treated by Fermi! and then in somewhat more 
detail by Arley. Fermi showed that for suffi- 
ciently slow neutrons the scattering cross section 
o may be calculated by substituting a relatively 
wide and shallow neutron-proton interaction 
potential for the actual narrow and deep po- 
tential. Then the Born approximation may be 
applied to this shallow potential. The width of 
this potential must be small compared to the 
wave-length of the neutron (relative to the 
proton). The effect of the substitute potential is 
essentially the same as that of a delta-function 
interaction between neutron and proton: since 
the neutron wave function varies but little over 
its width. 

If the proton is elastically bound to an infinite 
mass the scattering cross section depends on the 
ratio of vibrational quantum to neutron energy. 
For neutrons of energy less than a vibrational 
quantum, only elastic collisions are possible and 
it has been shown'? that the cross section is 
nearly independent of the vibrational proper 
functions. If such be the case, then according 
to perturbation theory the transition probability 
depends only on the density of energy states 
into which the neutron can be scattered. This 
density is proportional to p*dp dE, where p is 


; * Now at Purdue University, Lafayette, Indiana. 
1E. Fermi, Ricerca Scient. 7, No. 2, 13 (1936). 
2.N. Arley, Kgl. Dansk. Vid. Selsk. 16, 1 (1938). 


the momentum of the neutron in the center of 
gravity coordinate system and £ is the relative 
kinetic energy. In terms of p'dp dE 
where y is the reduced mass of the neutron plus 
proton-system (i.e., the system of the proton 
and all the particles to which it is bound). If 
the proton is initially at rest, is related to the 
initial kinetic energy, /o, of the neutron in the 
laboratory system by E=(u,m)Eo. Thus the 
transition probability, and therefore the scat- 
tering cross section, is proportional to u*. Since 
we are considering a proton-system of infinite 
mass, w is equal to the mass, m, of the neutron. 
The cross section for this case will be called o,. 

For the scattering by a free proton at rest the 
cross section is again proportional to y®, but in 
this case the reduced mass is p=m 2 and the 
cross section is (})o,. Thus we are led to the 
well-known factor of | for the ratio of the cross 
section for a free proton to the cross section for 
a proton rigidly bound to an infinite mass. 

In the case of the scattering of neutrons by 
hydrogen-containing molecules (in the gaseous 
state) the proton may be treated as free if the 
neutron energy is very great compared to both 
rotational and vibrational quanta (i.c., h times 
rotational and vibrational frequencies) of the 
molecule. When the neutron energy is much 
smaller than both of these energies, only the 
translational motion of the molecule as a whole 
is affected by the neutron impact so the cross 
section may be calculated as above with yu 
determined by 1/u4=1/m+1/M where M is the 
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total mass of the molecule. The cross section 
is then (\/,m+A/)*e, (as long as interference 
effects may be neglected). 

In the intermediate case for which the neutron 
energy is great compared toa rotational quantum 
but too small to excite molecular vibrations, the 
situation is somewhat more complicated. For, 
if it is assumed that the rotational-translational 
motion of the molecule can be treated classically, 
the effective mass of a proton in the molecule 
during the collision depends on the direction of 
the initial and final momentum of the neutron. 
For example, a neutron in a head-on collision 
with one of the protons of the H» molecule will 
feel an effective proton mass of 2m if it is moving 
parallel to the H-H axis and a mass m if its 
path is perpendicular to this axis. In general, 
there will be a linear relation between the change 
in velocity of the proton and the change in 
momentum of the neutron due to the collision. 
Therefore we may associate with the proton a 
hypothetical mass point which moves with the 
velocity of the proton but whose mass is a tensor 
which depends on the mass and structure of the 
molecule. 

In order to calculate the transition probability 
for the scattering of neutrons with energy that 
is great compared to the rotational quanta of 
the molecule, the rotational wave function in the 
final state can be replaced by a rotating wave 
packet which is a combination of neighboring 
rotational states since many rotational quanta 
will be excited by the collision. By making the 
packet sufficiently small the rotational wave 
function within it may be approximated by a 
plane wave. The molecular wave function is a 
product of the rotational and translational wave 
functions of the molecule multiplied by a 
vibrational wave function which in this case 
serves only to limit the position of the proton 
with respect to the other atoms in the molecule. 
Taking into account the fact that the transla- 
tional wave function of the molecule is a plane 
wave, we see that the relevant part of the 
molecular wave function after collision may be 
treated as a plane wave which describes the 
properties of the wave packet under translation 
and infinitesimal rotation. 

Insofar as the calculation of the matrix 
clement for the transition is concerned, the 
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initial rotational state may also be replaced by 
a plane wave since its only contribution to this 
matrix element is in the region of the wave 
packet. Therefore, just as in the case of free 
protons, the matrix element serves only to 
guarantee the conservation of momentum and 
angular momentum and the transition proba- 
bility depends only on the density of final energy 
states into which the neutron can be scattered. 
By the conservation laws, the density of the 
neutron states is equal to the density of the 
permissible final states for the molecule. How- 
ever this density will now depend on the reduced 
mass that is associated with the mass of. the 
neutron and the mass tensor of the hypothetical 
mass point described above. 

Since the quantum calculations are essentially 
concerned with small wave packets moving 
under the influence of constant or slowly varying 
potentials and since the only quickly varying 
interaction—that between the neutron and 
proton—leads to matrix elements that merely 
represent conservation laws, the entire calcula- 
tion can be carried out according to classical 
mechanics. The meaning of the mass tensor as 
the linear relation between velocity and mo- 
mentum vectors is the same in the classical and 
quantum treatment. Instead of the density of 
energy levels the volume of momentum space 
per unit energy interval will enter into’ the 
classical theory. 

The condition that the energy of the neutron 
be greater than the separation between rotational 
states is roughly A<7o, where A is the neutron 
wave-length and ro is a measure of the linear 
dimensions of the molecule. Therefore, under 
this condition, interference between the scatter- 
ing from two different protons in the molecule 
may be neglected. 


2. Tne Mass TENSOR 


Since the hypothetical mass point has only 
three degrees of freedom and the molecule 
(treated as rigid) has five or six, the motion of 
the molecule cannot be described completely by 
the motion of the mass point. 

We separate the infinitesimal displacements 
of the molecule into the set of two or three 
displacements for which the proton remains fixed 
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and the set for which the angular momentum 
about the proton is zero. In a collision with a 
thermal (S) neutron the angular momentum 
around the proton is unchanged. Therefore the 
change in motion suffered by the molecule 
involves only displacements of the latter type. 
The velocities associated with these displace- 
ments are uniquely determined by the velocity 
components of the proton so the behavior of 


these displacements can be equivalently de- . 


scribed by the motion of the mass point to be 
associated with the proton. 

It will be shown below that the two sets of 
displacements described above are orthogonal 
in the sense that the kinetic energy of the 
molecule can be written as the sum of two 
terms depending on the first and second types of 
displacement, respectively. 

The molecular motion is usually described by 
means of its angular velocity w about the center 
of gravity, and the velocity V. of the center of 
gravity. In order to obtain for the displacements 
per unit time the two types described above we 
write w and V, as the sum of two terms: 


w=w+w"’, V.=V'+V.", 


where the set (w’, V,’) satisfies 


[w’Xr]+V.’=0 (1) 


and the set (w’’, V.”’) satisfies 


Iw” — ]=0, (2) 


if r is the vector drawn from the center of 
gravity of the molecule to the proton, MW» is 
the mass of the molecule, and I is the moment 
of inertia tensor of the molecule. Equation (1) 
corresponds to the condition that the proton be 
at rest and (2) to the condition that the angular 
momentum about the proton be zero. Thus we 
have written the original displacement described 
by w and V, as the sum of two displacements 
(w’, V.’) and (w”, V.”’) which satisfy the condi- 
tions for the displacements of the first and 
second type. 
The kinetic energy of the molecule is 


= (w’-Iw’) + MoV.2+(w" Iw") + MoV!” 
] 
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and the condition for orthogonality, 
=0. (3) 


To verify Eq. (3) we make use of Eq. (2) and 
find that 


(w’-Iw”) = Mo(w’- 
= M.([w’ Xr]-V.”). 


Then by Eq. (1) the condition (3) is satisfied. 

We now wish to describe the part of the 
motion corresponding to (w’’, V.”’) by the motion 
of the hypothetical mass point. Assume that the 
momentum of the mass point is equal to the 
momentum of the molecule, 


P= M.V.”. (4) 


Conservation of momentum will then be satisfied 
for the molecule if it is satisfied for the mass 
point since the motion (w’, V,’) is unaffected by 
the collision; and, for the same reason, the 
angular momentum around the proton will be 
conserved. 

Since the mass point is defined as having a 
velocity V equal to the velocity of the proton, 
the inverse mass tensor is defined by V=M™'P. 
Now V=[w”’ Xr]+V-”, so Eqs. (2) and (4) give 


(5) 


which is the desired linear relation between V 
and P. 

It immediately follows from Eqs. (2), (4), and 
(5) that the kinetic ei. orgy: of the mass point is 
equal to that part of the molecular kinetic energy 
due to (w”, V.”’); that is, 


3(P-V)=3(w"- Iw”) +3. MoV”. 


Finally, the expression (5) written in terms of 
components is 


P;——P;-—-—P,, 

I, I; Mo I; 
where the subscripts 7, 7, & refer to the directions 
of the principal axes of inertia and 7, re, 73 are 
the coordinates of the proton under consideration 
with respect to the center of gravity. Therefore, 
the inverse mass tensor is 


rie 1 ; 
I, I; Mo I 


It is seen that the tensor is symmetric, a fact 
which can be proved directly from Eq. (5); 
thus it can be transformed to diagonal form. 

It will be convenient to introduce the dimen- 
sionless tensor n= mM. In the following para- 
graphs we shall calculate the characteristic 
values 11, M2, n3 of n for a few cases. 


H; 


For the hydrogen molecule, the principal axes 
of inertia and mass coincide. As stated in 
Section 1, the mass perpendicular to the H-H 
axis is m and that parallel to the axis is 2m. 
Thus for both protons. 


CH; 


In methane the moment of inertia ellipsoid is 
a sphere so the axes of inertia may be chosen in 
such a way that one of them passes through the 
proton under consideration. These axes are then 
the principal axes of mass for this proton. The 
mass along the axis that passes through the 
proton and the center of gravity of the molecule 
is the total mass of the molecule; thus the 
corresponding component of n is m;= 7g. The 
other two components of n are n2=n3= 4%. 


NH; 


The principal axes of mass do not coincide 
with those of the moment of inertia for ammonia, 
but one of them is the line joining the given 
proton with the center of gravity; another is 
perpendicular to this and in the plane determined 
by the proton, the nitrogen atom, and the center 
of gravity; and the third is perpendicular to 
these two. The mass for the first of these direc- 
tions is again the total mass of the molecule so 
ny =1!,7(=0.0589). In sufficient approximation 
one can set the HNH angle in ammonia equal to 
the HCH angle in methane (109° 28’). Under 
this assumption, the other two characteristic 
values of n are given by n:=0.604 and n3;= 0.439. 


H:0 


One principal axis of the mass tensor for water 
is, as usual, the line joining the proton with the 
center of gravity. The other two are perpen- 
dicular to this (and to each other) and one of 
them lies in the plane of the molecule. If the 
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HOH angle is taken to be 106°, the results for 
the components of n are m,;=!;s(=0.0555), 
n2=0.535, n3=1. 


3. SCATTERING Cross SECTION 


The collision of a slow neutron with a proton 
in the molecule gives rise to changes in w and 
V. which are represented by vectors of the type 
(w’’, V.""). The density of final neutron states, 
and therefore the transition probability, is then 
independent of motions of the molecule of the 
type (w’, V.’) and depends only on the motion 
of the mass point. 

If the mass M of the hypothetical mass point 
has the characteristic values 1/;, Mes, M3, the 
position of the center of gravity of the neutron- 
mass point system is defined by (m+ 
=mx;+M,t;, t=1, 2, 3. Here, x; are the co- 
ordinates of the neutron and &; those of the 
mass point. The velocity, momentum and kinetic 
energy of the two particles in the center of 
gravity system may now be calculated in the 
usual way. The total kinetic energy Eo separates 
into the sum of the kinetic energy of the center 
of gravity E, and the relative kinetic energy E. 

The density of final states for the neutron is 
proportional to dr/dE, where + is the volume in 
momentum space contained within a surface of 
constant relative energy (£). That is, 


r= ferupaps 


E =const. 
where (p1, P2, Ps) is the final momentum of the 
neutron in the center of gravity coordinate 
system. E is the sum of the kinetic energy of 
the neutron and the kinetic energy of the mass 
point (both taken in the center of gravity 
system). The former energy is (p-p), 2m and 
the latter is (3)(p-M~'p) since p=(pi, Ps) is 
also the negative momentum of the mass point 
in the C.G. system. Therefore the surface of 
constant energy is given by 
2mE = (p-p)+(p-np) = (p- =const., 

where p-'=1-+n. The tensor wu can be considered 
as the reduced mass tensor measured in units 


of the mass of the neutron. When yp is in diagonal 
form, the energy surface is the ellipsoid 


2mE= pr 
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TABLE I. The limtts o fT of Eq. (7) are given by 
En>Erot, En>(ErotkT)', and En<Eyin where E,, ts the 
neutron and T is the of the gas. 


MOLECULE IIe CHa NHs 
(ev) 0.015 0.001 0.002 0.006 
Ey ip (ev 0.5 0. 1 5 0.1 0.2 


and the volume, 7, is 
so that 


dr dk= 


The transition probability is therefore pro- 
portional to (uie2u3)'/?. Now, by conservation 
of energy, E=Ey—E, where Ey and E are the 
total energy and the energy of the center of 
gravity before collision. If P, and P are the 
momenta of the neutron and mass point in the 
laboratory system before collision, 


2mky=P,2+(P-nP) 


and 


2mE=(P-nuP), 


where P=P+P, is the total momentum. Taking 
the difference between these equations we obtain 


=(w(P,,—nP) 


Thus in terms of the components of momentum 
along the principal axes of mass of the molecule 
the transition probability is proportional to 
1 

The scattering cross section is given by the 
transition probability divided by the initial 
velocity of the neutron in the laboratory system. 
Consequently 


I 


. 


— 

where ¢ is the constant of proportionality. When 

the proton is rigidly bound to an infinite mass, 

Ny =Ng=n3=0 and pi =p2=pn3=1, so the cross 
section is o,=¢. 

If the vector with components 1;P; has the 
magnitude p and the direction cosines B;, then 
the energy of the mass point is given by 2m, 
=>) The cross section ¢ may be 
expressed in terms of p and £;, and in terms of 
the magnitude P,, and the direction cosines a; 


of the initial momentum of the neutron by the 
equation 


AND E. 


TELLER 
ELLE 


TABLE II. Scattering cross sections per proton in the given 
molecules as a function of neutron energy E,, and the temper- 
ature T of the molecular gas. The last column gives the region 
over which the formulae are valid. In addition E, must 
satisfy the conditions given in Table I. The total cross section 
per molecule is obtained by multiplying & by the number of 
protons in the molecule and adding G;, (Eq. (14)). 


MoLecuLe CROss SECTION PER PROTON Limits oF VALIDILY 


O=kT/E, 
O=kT/E, <0.1 
O=kT/E,, <0.05 


CH, G/on=2. 7 +0. 42kT/E,, 
NHs G/on=2.18 +0.39kT/E,, 
H.O 


This equation is then the generalization of the 
proportionality of the cross section to the square 
of the reduced mass. For use in specific problems 
such as the calculation of the mean scattering 
cross section or the mean free path, expressions 
obtained from (7) have to be averaged over the 
angles a; and 8;, over the velocity distribution 
of the molecules, and often over the neutron 
velocities. In order to carry out these calcula- 
tions, it is necessary to know the experimental 
conditions. We shall first calculate the mean 
cross section for the scattering of monochromatic 
neutrons. This quantity is connected with the 
diminution in intensity of a neutron beam in 
passing through a sufficiently small amount of 
gas. The effective cross section will also be 
calculated for a neutron current distribution of 
the form exp(—P,2 2mkT,,)dP,. Although 
this distribution is difficult to realize experi- 
mentally, we will make use of it since it leads 
to relatively simple results. 

It has been pointed out in the introduction 
that the formula (7) can be expected to hold 
only if the neutron energy is large compared to 
the separation between rotational states and 
small compared to any vibrational quantum. 
In Table I are given the values of Eyot=h? Tin 
and Eyjn=hymin for those molecules which will 
be considered below. 7,,;, is the smallest moment 
of inertia of the molecule and »,,jn the smallest 
vibrational frequency. The condition that the 
neutron energy be large compared the 
separation between rotational states is then: 
En>(ErotkT)', where T is the temper- 
ature of the molecular gas under consideration. 
In order that £,, be small compared to a vibra- 
tional quantum we have < Lyin. 


) 
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4. MEAN Cross SECTION FOR MONOCHROMATIC NEUTRONS 


The direction cosines a;, 8; are given with respect to axes fixed in the molecule. Therefore, in 
order to obtain the mean cross section, the expression (7) must be averaged over all possible orienta- 
tions a;, B; of the molecule. For the purpose of carrying out this average, let 


) 
Ay 


p 2 
a;-—n} 
(« B Y 


where the average is taken over a@;, 8;. Then ¢ is given by 
If the quantity Y=S0; 4,N; ¥. yi, is less than one, the square root may be expanded in powers 
of VN and then averaged. Thus 


and 


The conditions under which this expansion is valid will be discussed below. The ranges of values of 
kT FE, for which it holds are given in the last column of Table If. 
Now making use of the averages 


1 1 1 
2n+1 (2n-+1)(2n+3) 105 
we find 
Li —p)p*/P,,’, =9, 
(Si 4451011 +a2p?/P,*}, 
where 
3D >; 
(9) 


In obtaining these expressions, use has been made of the relation 1 y,=",;+1. With a=1 4-1, 
the average of the cross section over a;, 8; becomes 
| 
— +b2(1 — p?/ P..*)(p?/ 


On = Ox 1— 


(10) 


If the scattering molecules are at a temperature* 7, their energy distribution is proportional to 
exp(—Ey kT) times the element of volume in momentum space. The kinetic energy Ey of the 


8 It is assumed throughout the paper that the temperature 7 is not high enough to excite molecular vibrations, 
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molecule is the sum of £, (kinetic energy of the mass point) and a kinetic energy that is independent 
of the motion of the mass point. The element of volume in the molecular momentum space is the 
product of an element of volume in the momentum space of the mass point and a second volume 
element that is independent of the momentum of the mass point. Therefore the number of molecules 
for which (2mE,)' lies between the values p and p+dp is proportional to exp(—p? 2mkT)p*dp. The 
expression (10) must now be averaged over this distribution. If we let n=p P,, this average is 


(11) 


where 


0 


1 2 
-—f (1+an?)-! exp[ Yan? f expl —(E,,/kT)n? ]n*dn, 
0 


1 Zz 
(1+an*)-*? —(E, /kT) ] 


X (Din? f exp[ —(E,./RkT)? ]n*dn. 
0 


These integrals may be put into a convenient form by writing 


43) 


|] (14+2")!x? exp[ — dx? Jdx 
and 
Vivo 
where 


v=E,,/akT. 


Then the terms 7», 7), 72 become 
To= To, 
T2= 22/a — bel 23/a? — 24/a*) 


(12) 


Unfortunately, the integrals Jo, Z,, could not be carried out explicitly so that it has been necessary 
to evaluate them numerically. A rough evaluation of J) has been carried out as follows: In the 
range 0=x=0.5, the expression (1+)! may be expanded in power series and the resulting expression 
integrated. Similarly, from x=2 to «, (1+1/x*)! may be expanded and the direct integration 
carried out. Then in the range from 0.5 to 2 the integration has been done numerically. 

The integrals J;, have been treated similarly, except that in this case the ranges chosen were 
0=x=0.1 and 10=x= ~~, because of the slower convergence of (1+2*)~!. Only upper limits have 
been obtained for J:, in order to check the convergence of the expression for the scattering cross 
section. 

The values of n; may be taken from Section 2 and the corresponding values of y;, wu, a, @1, de, di, be 
may be calculated from them by Eqs. (9). The results for the 7; are then given by Eq. (12) and 
é may finally be obtained from Eq. (11). 
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The final result for the scattering of slow neutrons with energy E, by one of the protons in Hy 
is given in Fig. 1 for 0.1=£,/k7=10. The ordinate in this figure is ¢ divided by the cross section, 
on, of a free proton at rest (oy =je,). The abscissa is log, /kT where T is the temperature of the 
molecules. In computing these curves 7, and all higher terms have been neglected, which is a rather 
good approximation for Hy. The broken part of the curve was interpolated with the help of rough 
integrations of J») and Ip,. 

The validity of this curve depends upon the convergence of the expansion (8). In the case of He, 
it turns out that for all values of E,,E, (E, is the energy of the mass point), the quantity X in 
Eq. (8) is less than one in absolute value so the expansion can be expected to converge in this case. 
However, for other molecules such as CH,y, NHs3, and H.O, X is less than one only for a restricted 
range of values of E,/E,. In these regions, the quantities ¢/oy are, to a very good approximation, 
linear functions of kT E,. Formulae for the scattering cross sections per proton in the molecules 
CH,4, NH3 and H,O and their limits of validity are given in Table II. It is to be noted that the 
conditions given in the last column of Table II indicate that these results may only be applied to 
gases at low temperatures if the limits on £, given in Table I are taken into account. 

The scattering cross section per proton for H, in the range O=k7) E,,=0.1 is also an approxi- 
mately linear function of k7/E,. It may be expressed by the equation 


/on=1.22+0.86kT, E,. (13) 


The total cross section per molecule is given by ¢ multiplied by the number of protons in the 
molecule plus the cross section for the heavy nucleus in the molecule. This latter cross section 
may again be computed by associating a mass tensor with the motion of the heavy nucleus. The 
characteristic values of the mass tensor will lie between the mass of the molecule and the mass of 
the heavy nucleus. In CH,, the symmetry of the molecule causes the mass tensor to be spherically 
symmetrical and to be equal to the molecular mass. For HxO and NH; the mass tensors are still 
nearly spherical and they may be sufficiently approximated by using a spherical mass tensor the 
characteristic values of which are chosen equal to the average of the characteristic values of the 
correct mass tensor. This gives an effective mass of 16 m for the nucleus C™ in CH,, an effective 
mass 16.6 m for the nucleus N' in NH3, and an effective mass of 17.06 m for the nucleus O" in 
H,O. If we designate the actual mass of the nucleus by M and its effective mass by M,, then inte- 
gration over the motion of the mass point representing the heavy nucleus yields‘ 


1 
M.+m/ \M 


and f exp[ —x* Jdx. 
0 


where 


ox is the scattering cross section for the heavy nucleus when it is initially at rest. 
The total cross sections (¢7) per molecule are given in Table III for very low temperatures, 


that is, for k7/E,, ~0. 


5. EFFECTIVE Cross SECTION FOR NEUTRONS WITH A MAXWELL-LIKE DISTRIBUTION 


We now assume that the number of neutrons with momentum between P, and P,+dP,, passing 
through unit surface per unit time is proportional to P,? exp(—P,?/2mkT,,)dP,,, where the effective 
temperature 7, is not in general equal to the temperature T of the scattering gas. Then the expression 
(10) for os, must be averaged over this distribution as well as the velocity distribution of the mass 
point. When the expansion (8) is valid (cf. Table II) the average cross section may again be written 


* Compare J. Schwinger, Phys. Rev. 58, 1004 (1940). 
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in the form 


with 


AND €. 


TELLER 


| f (1+ap?/P,")) ( + ) / 


1 


ete., where J), Js, J; are integrals of the form 


p 


f f ( ) (12’) 
0 0 2mk 


16 1 
0 2mkT,, 


These integrals can be evaluated explicitly and the result’ for \>a? is 


a 
m(A—a*) 
JitJs= — 15)  arecot(A/a? —1)3 
1 
+ + 8A" + — 
ad? 
a 
Jo= |  arccot(A/a® —1)!——(134 2a”) }. 
m(A—a*)* d 


For \<a? all of the terms (A—a?) ! arecot(A/a?—1)~) are to be replaced by 


(a?— arccoth(1 —A/a*) 


6. APPLICATION TO EXPERIMENT 


In applying the foregoing calculations to 
experiments a number of limitations, most of 
which have been previously mentioned, must be 
borne in mind. First, the entire theory is limited 
to neutron energies satisfying the conditions 
E,>Ert, En>(ErtkT)', and E,<Eyi, where 
Eo, and Eyi, are given in Table I. In addition, 
the particular approximate methods applied in 
Sections 4 and 5 are limited in every case except 
that of He by the conditions given in the last 
column of Table II. In the particular case of 
H. the cross sections have been roughly inter- 
polated in the part of Fig. 1 which is indicated 
by a broken curve. 


5 Cf. N. Arley, reference 2 (note 6), for a discussion of 
one of the integrals of this type. 


The greatest difficulty arises from the fact 
that the neutron current distribution is not 
known. It has been assumed in Section 5 that 
it is of the form P,”? exp(—F, kT,)dP,. This is 
the distribution that one might expect if the 
neutrons leave the free surface of a slowing down 
material in which the time required for the 
neutrons to reach the Maxwell distribution is 
small compared to the lifetime of the neutrons 
in the material. On the other hand, a distribution 
of the form P,’ exp(—E,/kT,) would be ex- 
pected if the neutrons behave like a gas at 
thermal equilibrium escaping through an opening 
of linear dimensions small compared to the mean 
free path. It is probable that the usual experi- 
mental situation corresponds to a condition that 
is between the two. However only the first case 
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Fic. 1. Cross section for the scattering of neutrons with 
energy £, by one of the protons in Hy» gas at temperature 
1. ow is the cross section for a free proton initially at rest. 


The broken part of the curve was interpolated with the 
help of rough calculations. In the range O=k7/E,=0.1 
the cross section is given by Eq. (13). 


was treated here because the second apparently 
involves a double numerical integration. 

Preliminary experiments on the scattering of 
300°KK neutrons by Hz and CH, at room temper- 
ature have been carried out by Carroll and 
Dunning.® According to Table I], we see that 
the present theory cannot be applied to the 
experiments on CH, since k7~1. In the 
case of He, the neutron energy, 3k7,, is approxi- 
mately equal to the limiting energy /yo. given 
in Table I, so we can only expect a rough 
agreement. 

This experimental value of the cross section 
per proton for scattering by He is given in the 
first column of Table IV. In the second column 
the theoretical value of @ oy is given on the 
assumption that the neutrons have a Maxwell- 
like distribution (Section 5). The third column 
contains the value of oy obtained from the first 
two columns. In the fourth column is the 
experimental value of oy, found by Cohen, 
Goldsmith and Schwinger? and by Hanstein,* 


®H. Carroll and J. R. Dunning, Phys. Rev. 54, 541 
938) 


7V. W. Cohen, H. H. Goldsmith and J. Schwinger, 
Phys. Rev. 55, 106 (1939). 
P 044) B. Hanstein, Phys. Rev. 57, 1045 (1940); 59, 489 
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TABLE III. o7 is the total cross section per molecule for 
the scattering of monochromatic neutrons if the thermal 
energy of the molecules is much smaller than the energy of 
the neutrons (RT /E,, ~0). on, oc, on, and ag are the scattering 
cross sections for the free nuclei H, C, N and O when they 
are tnittally at rest. 


MOLECULE oT MOLECULE oT 


2.440n NH; 6.550 +1.020y 
CH, 10.867 +1.040¢ H,0 


Taste IV. The first column contains the experimental 
values of & obtained by Carroll and Dunning,® the second 
gives the theoretical value of €/0n assuming a neutron current 
distribution of the form P,? exp(—P,2/2mkT,)dP,,. on has 
been calculated from these two. on ts the value given by 
Cohen, Goldsmith and Schwinger’ and by Hanstein,® and 
on ts the value given by Schwinger.** 


THEO- 
EXPERIMENTAL RETICAL EXPERIMENTAL 
o/oHn on on” 


32X10 cm? 1,69 18.9X10-% 21%K10° 16.6 


and the last column gives the value obtained by 
Schwinger’ from the experiments of Alvarez 
and Pitzer® with very slow (20°KK) neutrons. 

If it had been assumed that the neutron current 
distribution was of the form 
rather than P,Zexp(/, k7)dP,, then greater 
weight would have been given to higher values 
of £, in calculating the theoretical cross section. 
It can be seen from Fig. 1 that this would tend 
to decrease the theoretical value of @ oy and 
thus increase the value of oy given in Table IV. 

It would be particularly interesting to measure 
the cross section for the scattering by He of 
monochromatic neutrons with an energy some- 
what greater than £,.: (Table 1). If this measure- 
ment were carried out at different temperatures 
of the He gas it would be a means of checking 
the theory and should yield a fairly accurate 
value of oy. 

The authors take pleasure in thanking Pro- 
fessor Placzek for valuable discussions. 


9L. W. Alvarez and K. S. Pitzer, Phys. Rev. 58, 1003 
(1940). 
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Corrections to the nuclear radius and particle density arising from the neutron excess N—Z, 
from the electrostatic repulsion of the protons and from the surface tension forces, are calculated 
on the basis of the statistical model. The density throughout the interior of the heaviest nuclei 
is found to be about one-third less than in medium light nuclei (A ~50). The trend of ro=R/A! 
from the lightest to the heaviest elements is discussed and compared with experiment. 


INTRODUCTION 


HE law that the nuclear volume is propor- 
tional to the total number of particles in 
the nucleus is an approximate theoretical conse- 
quence of the statistical treatment of nuclei with 
exchange-saturation type forces acting between 
the particles. Let R denote the radius of the 
sphere into which the particles are packed with 
constant density; then R=7,A! where A is the 
mass number and fo, representing the inter- 
particle distance, depends only on the ratio N/Z 
of neutrons to protons. The monotonic increase 
of this ratio along the periodic table introduces 
a correction term to the density which, although 
constant within a given nucleus, will decrease 
slowly with increasing mass number. The kinetic 
energy part of this correction is due essentially 
to the exclusion principle and the potential 
energy part to the symmetry and saturation 
properties of the Hamiltonian, and both parts 
raise the energy and diminish the density as N /Z 
increases. The correction to the density depends 
on the nuclear compressibility which can be 
explicitly evaluated from the statistical model. 
Through the introduction of surface effects due 
partly to the desaturation of the short-range 
forces near the boundary of the nucleus and 
partly to the diffuseness of the boundary associ- 
ated with the uncertainty principle, the density 
throughout the nucleus will be altered. Omitting 
from consideration the lighter nuclei (A50), 
we can assume the same thickness of surface 
layer for all nuclei.! The density falls off con- 
tinuously from the constant value in the interior 


* On leave of absence from Purdue University. 

14~50 marks the minimum of the packing fractions 
curve, where surface effects have become small compared 
to volume effects. 
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over an effective distance approximately equal 
to the range a of nuclear forces. The ‘‘geometrical 
radius” of the nucleus, including the surface 
region, will differ from the ‘‘constant density 
radius” by approximately a,2~X~10-" cm. 
The surface tension forces uniformly contract the 
compressible nuclear matter and the resulting 
density in the interior is greater and the radius 
R smaller than their Majorana values. The 
correction to the density decreases with in- 
creasing mass number since the surface pressure 
is inversely proportional to the radius. Conse- 
quently, the density in the interior diminishes as 
A increases. For light nuclei (A+50) the surface 
region can no longer be properly differentiated 
from the interior. If the density of particles be 
represented by a Gaussian distribution about the 
center of the nucleus, the nuclear radius can be 
defined in terms of the breadth of this distribu- 
tion and the A} law will still hold approximately 
true. In the transition to light nuclei, the in- 
creased density in the interior due to surface 
tension is compensated by the reduced density in 
the surface layer. Accordingly, one can even 
expect that the average density is less and the 
average ro (to be denoted by re) greater as one 
goes to smaller mass numbers. 

A third correction to the A! law arises from 
the mutual electrostatic repulsion of the protons. 
This effect, which has already been considered 
by several authors,*-> obviously leads to an 
increase of ro with increasing Z. An additional 
effect of the Coulomb forces is to distort the 


2S. Fliigge, Zeits. f. Physik 96, 459 (1935). 

3K. Nakabayasi, Zeits. f. Physik 97, 211 (1935). 

*H. A. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 
82 (1936). 

5 E. Feenberg, Phys. Rev. 59, 149 (1941). 
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constant density distribution : the proton density 
tends to increase toward the boundary of the 
nucleus and the correlated neutron density varies 
in a similar manner. This effect has been dis- 
cussed by Feenberg® who has shown that al- 
though the departure from uniform density is 
appreciable, the associated energy correction is 
small and the correction to the radius negligible. 


Save for light nuclei, all three corrections tend — 


to increase 79 with increasing A. The calculations 
to be described give for the constant density 
radii: 


(1) 


where the corrections in order represent the 
effects of NZ, of surface tension and of electro- 
static repulsion. The factor in brackets increases 
by 9 percent from Fe®® to Pb**?, The Gaussian 
radii for light nuclei are given by: 


(2) 


The factor in brackets has a 7 percent increase 
from Sc" to N“; when the effects of neutron 
excess and Coulomb repulsion are included, this 
is reduced to 5 percent. Values of r*) and r*_ may 
be determined empirically. 


CALCULATIONS 


Bethe and Bacher‘ have outlined the calcu- 
lationsof the statistical model, assuming Gaussian 
functions for the intranuclear forces. Their nota- 
tion will be used as far as possible. The total 
energy E of the nucleus in its ground state is 
expressed as a function of the adjustable param- 
eter x which is proportional to the cube root of 
the density : 


E(x) 
(3) 


The first term gives the main binding energy 
when N=Z, the second term the correction when 
N#Z, the third term the surface energy, and the 
fourth term the classical Coulomb energy of the 
protons.? The A! law for the nuclear radius is 
obtained when only the first term is taken; 


°E. Feenberg, Phys. Rev. 59, 593 (1941). 
’ The exchange part of the Coulomb energy has been 
omitted here, since its effect on the results is negligible. 
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when this is minimized with respect to x, we 
obtain the standard density xo*. The inclusion of 
the three additional terms shifts the minimum 
of the energy from x» to x’. Let 


(4) 


where the c’s are all small. Since a’(x’») 
= (x’y—xy)a’’ (x9), we obtain: 


(xo) (x0) 
a’’ (x9) a" (xy) 
(5) 
5’ (x0) 
tog 
a’’(x») 


The quantity x?a’’(xo) is inversely proportional 
to the compressibility of nuclear matter; in fact, 
we have for the compressibility 


k= (x0). (6) 


From the statistical model with Gaussian poten- 
tials of interaction radius a, the following ex- 
pressions are obtained: 


a(x) = (3/10) Tox? —9A*x-*[.2 — 3x? 
+ (x*—2) exp(—x*) +2!x*b(x) J, 
B(x) = (1/6) Tox? + (A*/3x) 
[7x2 -9+(2x°+9) exp(—x*) ], 
x(x) =(9/8)(3/ 4) "a Tox (7) 
X [xt —4x2+6— (2x°+6) exp(—x*) ], 
5(x) =(2 5)(3,)e?/a)x, 


where 
x To=h?/Ma?; A*=A/8r' 


and A is the depth of the Gaussian well potential 
—A exp[—(r/a)*] for the ground state of the 
deuteron. (x) represents the error integral. The 
symmetrical Hamiltonian with the Volz* condi- 
tions has been used in the expressions for a(x) 
and B(x). The formula for y(x) has been given 
for the case that the auxiliary potential falls off 
linearly in the surface region over a distance ¢: 
the last two terms represent the corrections to 
the kinetic and potential energies due to the 


*H. Volz, Zeits. f. Physik 105, 537 (1937). The result 
for cg depends only slightly on the special choice of con- 
stants in the Hamiltonian. 
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redistribution of the particles and the first term® 
corrects for the fact that. the statistical model 
does not give an infinite kinetic energy for a 
sharp nuclear boundary. y(x) has to be mini- 
mized with respect to ¢. 

The ad hoc constants of the statistical model 
(reference 4, Eq. (164)) were used to calculate the 
corrections to the nuclear radius. Since these 
constants give too strong an interaction between 
the nuclear particles, they may be expected to 
give too large results for quantities other than 
the main binding energy which depend on the 
strength of the forces. Thus B(x») =32 Mev and 
¥(xo) =17 Mev compared with the semi-empirical 
values of 20 Mev and 13 Mev, respectively. The 
constants! determined from the hydrogen and 
helium isotopes and two-particle scattering ex- 
periments give too small results; in calculating 
the c’s a suitable average was taken with less 
weight being put on the latter constants since 
they give probably too small a value for the 
inverse compressibility (xy). From the ad hoc 
constants we obtain xya’’(x») =72 Mev, in agree- 
ment with semi-empirical estimates of Feenberg.* 

The values of ¢3 and c, have also been esti- 
mated in a quasi-classical way assuming that the 
N-—Z correction energy was inversely propor- 
tional to the nuclear volume and that the surface 
energy contained a surface tension factor which 
was independent of R. Substituting semi-em- 
pirical values for 8, y and x, the resulting values 
for cs and c, agreed within 20 percent with those 
calculated from the statistical model. 

The calculations of Fliigge? have been used to 
estimate the surface corrections for light nuclei. 
Assuming that the density of protons (and of 
neutrons similarly) could be represented by the 
Gaussian distribution 

Re) expl—2(r Re)? ], (8) 
he calculated the binding energy statistically 
with the ‘‘nuclear radius’’ Rg as an adjustable 
parameter. His results for a suitable choice of 
the constants lead to Eq. (2). 

APPLICATIONS 


Since the nuclear radius can only be defined 
precisely in terms of the operations used to 


°C. v. Weisziicker, Zeits. f. Physik 96, 431 (1935). 
'"W. Rarita and R. D. Present, Phys. Rev. 51, 788 (1937). 
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measure it and since these operations are widely 
different in character and theoretically unre- 
lated except in a very approximate way at 
present, it is clear that the intercomparison of 
nuclear radii derived by different methods from 
different parts of the periodic table can be 
somewhat misleading. To obtain reliable evidence 
for the effects discussed in this paper, it is 
necessary to have at least one method for meas- 
uring nuclear radii which can be applied over 
the entire periodic table. Measurements of the 
cross sections for elastic and inelastic scattering 
of fast neutrons would provide two suitable 
methods. In view of the experimental uncer- 
taintics in the latter measurements, we shall 
consider only the evidence from elastic scattering 
experiments. As was pointed out by Placzek" 
and Bethe,” the elastic scattering of fast neutrons 
can be described by the formula for Fraunhofer 
diffraction around a spherical obstacle.’ The 
differential cross section for the seattering is 
given by 


o(0)=[R’J ARR’ sind) sind |’, (9) 


where J; is a Bessel function of the first: kind, 
k is the wave number of the incident neutrons, 
and R’ the effective radius of the nucleus acting 
as a screen. Bethe" has shown that, because of 
the finite range of the forces extending beyond 
the nuclear boundary, the effective nuclear 
radius is increased for clastic scattering by 
about cm. The difference between the 
“screening radius” of the elastic scattering ex- 
periments and the “constant density radius” of 
the statistical model is independent of mass 
number, and is approximately R’—R~2X10 © 
cm. The experimental results of Wakatuki® on 
Fe and on Pb are well fitted by the diffraction 
formula (except for the measurements at 14°) 
with R’(Fe)=7.4X10-" cm and R’(Pb)=11.3 
X10-" cm. These two values are almost exactly 
in the ratio of the cube roots of the respective 
mass numbers; since R’>R, this provides direct 


1] am indebted to Dr. Placzek for a discussion of the 
elastic scattering experiments. 

2G. Placzek and H. A. Bethe, Phys. Rev. 57, 1075.\ 
(1940). 

8M. Born, Optik (J. Springer, Berlin, 1933), p. 160. 

4H. A. Bethe, Phys. Rev. 57, 1125 (1940). 
( a Wakatuki, Proc. Phys.-Math. Soc. Japan 22, 430 
1940). 
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evidence for the calculated inflation effect 
(ro(Pb) >ro(Fe)). The earlier data of Kikuchi, 
Aoki, and Wakatuki" are best fitted with R’(Fe) 
=6.5X10-" cm and R’(Pb)=10.2X10-" cm. 
Averaging the results of the two experiments, 
subtracting 2X10~-" cm and dividing by A!, we 
obtain ro(Fe) = 1.30 10-" cm and ro(Pb) = 1.48 
<10-" cm; an increase of 13 percent compared 
with the theoretical estimate of 9 percent from 
Eq. (1). Actually it is somewhat questionable to 
apply the diffraction formula in this region where 
the neutron wave-length is still comparable with 
the nuclear radius (A~2.5X10-" cm). Much 
interesting information about nuclear radii could 
be obtained by investigating the elastic scattering 
at higher neutron energies. 

The evidence from the penetration of potential 
barriers is less certain.’ Weisskopf and Ewing'® 
have derived a value of ro>=1.310-" cm from 
proton penetrabilities into elements with A ~60. 
It is unsatisfactory to compare this value with 
the 7» derived from the spontaneous alpha-decay 
because of the uncertainties in the interpretation 
of the latter. The Weisskopf-Ewing radii, taken 
together with the trend of ro discussed in this 
paper, would indicate that the radii of the 
heaviest nuclei are close to the Gamow values; 
hence, the many-body correction and the correc- 
tion for the finite radius of the alpha-particle 
must nearly cancel one another. Another method 
for estimating ro is from the Weizsicktr-Bethe 
semi-empirical formula with the constants ad- 

‘6 Kikuchi, Aoki and Wakatuki, Proc. Phys-Math. Soc.. 
Japan 21, 410 (1939). 

17 The “penetrability radii” may be expected to be 
rather close to the ‘‘geometrical radii’: although the 
Coulomb law breaks down somewhat beyond the nuclear 
boundary, the artificially sharp cut-off of the Coulomb 
potential assumed in calculating the penetration factor 
mostly compensates for this by giving too small a break- 
down distance. 


18V. F. Weisskopf and D. Ewing, Phys. Rev. 57, 472 
(1940). 


justed to fit known mass defects. The Dempster!® 
packing fractions for s¢Pd!*, and 
lead to a value of ro>=1.2X10-" cm, whereas 
The increase in ro lies within the limits of error 
and need not be real. 

The trend of ro from A~50 down to very 
light nuclei is also rather uncertain. The data on 
elastic scattering of neutrons by Al nuclei 
cannot be used to get an accurate value of ro 
because the measurements are not much above 
the background; also the wave-length X and the 
correction R’—R are almost as large as R itself. 
Values of 79 can be obtained from the Coulomb 
energy differences between the light positron 
emitters and their stable isobars. Assuming 
either a constant density distribution of the pro- 
tons, or else the Gaussian distribution (8), whose 
electrostatic energy is given by (1/1/7)Z%e?/Re, 
or any similar distribution, the experimental 
results’? from the positron spectra of elements 
between A = 11 and A =41 indicate a remarkable 
constancy of ro or rg over the entire region.” 
If the statistical estimate of the Coulomb ex- 
change energy is included, the values of ro and r¢ 
are found to increase slightly with increasing A 
in contrast to what is predicted by Eq. (2). It 
does not seem justifiable to make comparisons 
between the radii determined in this way and 
the ‘penetrability radii’ of Weisskopf and 
Ewing. 

I wish to express my appreciation of the 
hospitality of the Harvard Physics Department. 

1” \. J. Dempster, Phys. Rev. 53, 869 (1938). 

20 M. G. White, E. C. Creutz, L. A. Delsasso and R. R. 
Wilson, Phys. Rev. 59, 63 (1941). D. R. Elliott and L. D. 
P. King, Phys. Rev. 59, 403 (1941). 

21. Any smooth distribution in density of the form 
p=pof(r/R) will have an electrostatic energy proportional 


to Ze?/R and will exhibit the same agreement with the A} 
law as the constant density distribution. 
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An Electron Lens Type of Beta-Ray Spectrometer*? 


Ciirrorp M. WitrcHEeR 
Department of Physics, Columbia University, New York, New York 
(Received April 21, 1941) 


A beta-ray spectrometer utilizing the focusing action of the homogeneous magnetic field 
produced by a long solenoid is described. The theory of the action of an instrument of this type 
is considered in some detail for the case of a point source of beta-rays, and the extension of the 
theory to cover the case of an actual source of finite area is discussed. Details of design are 
given, together with a number of experimental tests which have been made on the instrument. 
These tests include measurements of the absorption of beta-rays in the Cellophane window of 
the G-M counter used for detection, measurement of the resolution and efficiency of the 
spectrometer, calibration of the instrument, and an examination of the beta-ray spectrum of 
radium E. The beta-spectra of phosphorus ,;?® and iridium ;;Ir™ have been investigated. The 
forms of the distribution curves all agree better with the Fermi than the K-U theory. It is 
concluded that the spectrometer operates satisfactorily and possesses some definite advantages 


over the conventional designs. 


INTRODUCTION 


HE most common methods for investigating 

beta-ray energies which have been in use 
for the past three decades are: (1) absorption 
measurements, (2) measurements with a cloud 
chamber in a magnetic field, and (3) measure- 
ments with a magnetic spectrometer. essentially 
like the type developed by Rutherford and 
Robinson! in 1913. 

The suggestion that a large solenoid might be 
used as an electron lens to focus beta-rays of 
determinate energy upon some sort of detecting 
device appears to have been first made by 
Kapitza* in 1923. In 1924 an instrument of this 
type was built by Tricker.* More recently 
Klemperer has constructed an electron lens type 
of spectrometer, using the focusing properties of 
a flat coil.® 

If one considers the operation of a spectrometer 
utilizing the electron lens properties of a solenoid, 
it soon becomes evident that, if a satisfactory 
instrument of this type could be designed, it 
would possess several marked advantages com- 


* Submitted in partial fulfillment of the requirements 
for the degree of Doctor of Philosophy in the Faculty of 
Pure Science, Columbia University, New York, New York. 

7 Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 
(1913). Rutherford and H, Robinson, Phil. Mag. 26, 717 

13 

2 P. Kapitza, Proc. Camb. Phil. Soc. 22, pt. 3. 

3 R.A. R. Tricker, Proc. Camb. Phil. Soc. 454 (1924). 

40. Klemperer, Phil. Mag. 20, 545 (1935). 

5 A flat coil type of electron lens spectrometer has been 
— recently by M. Deutsch, Phys. Rev. 59, 684A 
(1941). 
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pared to the conventional spectrometer,’ among 
which might be mentioned the following: (1) The 
efficiency of the spectrometer can be relatively 
high, since a larger solid angle at the source can 
be utilized without sacrifice of resolving power. 
(2) Source material can be distributed over the 
surfaces of thin films of collodion or some similar 
substance, covering an area much larger than 
that generally possible in the conventional type 
of spectrometer. It is thus possible for an instru- 
ment of this type to furnish data on materials 
whose activity cannot be concentrated suffi- 
ciently for use in the usual type. (3) Only one 
layer of absorbing material (the window of the 
G-M counter) is present in the electron path. 
(4) The source can be located relatively far from 
any scattering material except that of the source 
holder itself. (5) The relatively large distance 
between source and detector, together with the 
large amount of absorbing material interposed by 
a central baffle system, should practically elimi- 
nate detection of gamma-rays. 


THEORY OF SPECTROMETER 


A. General considerations 


The theory of the action of magnetic electron 
lenses has been treated at length by various 
workers.* * > The present discussion is confined to 


® Witcher, O’Conor, Haggstrom and Dunning, Phys. 


Rev. 55, 1135A (1939). 
7 W. Glaser, Zeits. f. Physik 109, 700 (1938). 
8J. H. McMillen and G. H. Scott, Rev. Sci. 
288 (1937). 
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the simple treatment of a few of the relations 
which have the most direct bearing on the use 
of the uniform solenoid lens for the examination 
of beta-ray spectra. 

We consider an electron which is initially 
projected with momentum P, 


in a direction making an angle a with the 
direction of a uniform magnetic field J7. The 
electron’s path will have the form of a circular 
helix of radius 
p=(p, He) sine (1) 

and pitch 

L=(2rp, Ile) cosa. (2) 
If s denotes the displacement of the electron from 
the initial position in the direction of 77, and r its 
displacement perpendicular to //, it is readily 
shown that: 


r=(L/m) tana sin(rz L). (3) 


To examine this equation it is convenient to 
introduce the variable 


D=L/(x cosa) = 2p, He, 
whereupon (3) takes the form: 
r=D sina sin[z/(D cosa) }. (4) 


With the condition that the field JZ is initially 
fixed, this equation indicates (1) how r varies 
with s along the path of any particular electron 
(D and @ constant), or (2) how, at a specified 
distance z from the initial position of the electron, 
the displacement 7 varies with departure angle a 
for electrons of the same momentum (z and D 
constant), or (3) how, for fixed zs and a, r varies 
with electron momentum. 

The first case is obvious, and case (3) need not 
be considered in this discussion. 

As @ increases from 0 to 7/2, z and D constant, 
r starts from 0 and oscillates between a succession 
of positive and negative maxima, these being 
located at the values of @ for which dr, da=0. 
The zeros of |r| lying between successive maxima 
are at the @ values for which 


z/(D cosa)=nnz, n=1, 2, 


Physically, this behavior of |r| corresponds to 
the fact that electrons with the same momentum, 


bid 
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Fic. 1. Geometrical design functions. 


but leaving the source with larger and larger 
departure angles, may cross the plane determined 
by z during the first, second, third, etc., turns of 
their helical paths. Further investigation indi- 
cates that the successive maxima of |r| in the 
direction of increasing @ are consecutively 
greater, approaching the limiting value D as a 
approaches 7/2. Physically this expresses the 
fact that, as @ increases, the pitch of the helix 
contracts, its diameter becoming correspondingly 
greater. For a= 7/2, the helix becomes a circle of 
diameter D=2p/ He. 


B. Intensity and resolution for a point source 


It is difficult to obtain a general treatment of 
the problems of intensity and resolution when the 
source of electrons has a finite area in the plane 
perpendicular to the magnetic field. However, the 
treatment for a point source may be used to 
suggest the features of the more general problem. 

Consider a source at the point (7, z) =0 emitting 
electrons with all possible momenta from 0 to 
some upper limit, the emission being assumed 
isotropic. We first wish to examine the distri- 
bution of the current arising from electrons whose 
momenta lie in the neighborhood of a given valuc 
p in the plane determined by assigning to z the 
value Z. From the analysis given above for the 
case z and D constant, it appears that at each 
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Fic. 2. Intensity function J. 


value of ry there will be contributions to the 
current arising from an infinite set of @ values, 
extending from some definite lower limit to 7, 2. 
Let I(r, a) denote the contribution to the density 
of current at r due to any specified value of @ 
from this set. Then J(r, a) must be proportional 
to sina, and also to da/dr. It therefore follows 
that J will become infinite at the points where 
dr da=0, viz., the maxima of |r| as function of a 
considered previously. The current distribution 
due to electrons of the same momentum will 
therefore be characterized by the existence of a 
set of concentric ‘focal rings,’’ upon which its 
density becomes infinite. The smallest ring will 
occur at the smallest value of r for which dr/da 
vanishes. If this derivative vanishes concurrently 
with 7, we may speak of a “point focus.”” We 
have, from (4): 


dr /da=D cosa sin{z/(D cosa) | 
+z tan’a cos[z/(D cosa) ]=0. (5) 


On setting y=2/(D cosa), this gives: 
a=tan"/[—(1/y) tany]. (6) 


The symbol @ is used to identify the @ value as 
that corresponding to a maximum of |r|. (It is 
convenient to consider y as a parameter in this 
and the succeeding two relations.) Combining 
the last relation with (4) and (5) and using 7 to 
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denote the maximum value of 7 thus obtained, 
we have: 


Z=([(siny), y] tana, (7) 
rD=(1/n)y!, \/(y—tany). (8) 


By use of the last three relations, with y as a 
parameter, the curves of Fig. 1 are obtained. 
These curves exhibit & (corresponding to the first 
maximum of 7) and #/2Z as functions of 2/7D from 
2 rD=0.5 to 1.0. (& is plotted in radians, on the 
same scale as that of 7/2.) If Z is regarded as 
fixed, these curves indicate the variation of & and 
of the radius of the first focal ring as functions of 
electron momentum. The first point focus occurs 
at Z/rD=1, a=0. As Z/rD decreases, this opens 
up into the first ring focus, as indicated in 
Fig. 1. For Z7rD=0.9, &=15° approximately, 
and 7#/2=0.020. For Z/rD=0.8, @=22.5° and 
2=0.063. 

Intensity relations—The current flowing 
through the annular area of radii 7— Ar and 7 will 
be proportional to 


a+ 
f sinada (sin&)(Aa;+Aa2), 


a—Aa; 


= 


Fic. 3. Resolution function. 
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where 
(9) 
The intensity function 
J=sina (10) 


is plotted against Z/rD in Fig. 2. This function 
gives an indication of the variation with Z/rD of 
the current to be expected through a ring of 
constant width Ar and outer radius 7, * being the 
radius of the first focal ring. Note that J is a 
maximum at Z/rD=0.6, and that it begins to 
drop rather sharply above Z7/7D =0.85. 
Resolving power.—A convenient criterion of 
resolving power is the rate of variation of D (i.e., 
of electron momentum) with 7, expressed as a 
fraction of D itself. We accordingly consider 


(1/D)(dD/d?) = (d/d?)(logD). (11) 


In Fig. 3, logio(#D 2) is plotted against 7/2. The 
resolving power, as indicated by the slope of this 
curve, is greatest for */Z=0, decreasing rapidly 
until #/Z=0.10, and then becoming approxi- 
mately constant. 


DESIGN OF SPECTROMETER 

A. Design criteria 

The utility of any beta-ray spectrometer will 
depend in general upon the following three 
criteria: (1) the range in energy over which the 
instrument can be used; (2) the sharpness of 
resolution ; (3) the efficiency, as measured by the 
fraction actually detected of the beta-rays which 
leave the source with energies lying in the band 
passed by the instrument. 


Fic. 4. Solenoidal beta-ray spectrometer. 


"COUNTER LEADS 


B. Maximum energy 
Equation (2) in the form 


HL /cosa=2rp/e 


indicates that the maximum momentum (and 
hence the maximum energy) which it is possible 
to measure with an instrument of this type will 
depend on (1) the maximum obtainable field 
strength and (2) the length and diameter of the 
region in which the field can be considered 
uniform; i.e., the region between source and 
counter. In fact, the maximum energy increases 
with each of these factors. 


C. Efficiency and resolution 


From the theory outlined above, the electrons 
are brought to a focus at the plane Z in a ring of 
radius *. The geometry of the instrument may 
therefore be determined by a diaphragm provided 
with an annular slit of suitable radius and width 
equal to a small fraction of its radius, together 
with a_ baffle system which will allow only 
electrons with the correct angle of emergence & 
to pass through the spectrometer. 

The consideration of an extended beta-ray 
source in the form of a disk perpendicular to the 
magnetic field shows that the momentum values 
passed by the instrument are a function of (1) the 
distance from the axis to the initial position of the 
electron and (2) the azimuthal angle of its initial 
direction of motion. Since the resolving power 
decreases somewhat from the center of the source 
to its edge, the final design is a compromise 
between the resolution, the efficiency, and the 
area of the active material. Equation (4) and the 
curves in Figs. 2 and 3 permit approximate 
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selection of the actual geometrical design factors 
for a given resolution and efficiency. 


DESCRIPTION OF APPARATUS 


A. General features of spectrometer 


The general features of the spectrometer as 
finally designed® are shown in Fig. 4. Essentially, 
it consists of a long solenoid having the source of 
beta-radiation placed on the axis near one end 
and a G-M counter for detection placed 90 cm 
distant in the analogous position near the other 
end. Between source and counter are interposed a 
series of Duralumin disks and rings, which serve 
to specify, within definite limits, the dimensions 
of the helical paths which the electrons may 
follow in passing from source to counter. 

The dimensions of the instrument were pri- 
marily determined by the power supply, a 600- 
ampere 30-volt d.c. generator, which was avail- 
able before design of the spectrometer was begun. 
The theoretical considerations outlined indicate 
that for a specified power supply, the solenoid 
should be designed for (1) maximum field 
strength, (2) maximum length of homogeneous 
field region, and (3) maximum ratio of diameter 
to length. 


B. Solenoid 


The solenoid was madé by winding six layers of 
j-inch diameter copper tubing (0.035-inch wall 
thickness) on a brass tube 10.5 inches o.d., 10 
inches i.d., and 66 inches long. The actual length 
of winding is 60 inches, and there are 220 turns 
per layer. The tubing for each layer was fabri- 
cated in one piece (lengths ranging up to 700 
feet) by the manufacturers and has double 
cotton insulation 4.0 mils thick. All six layers of 
the winding are connected in parallel electrically 
by bus bars and also in parallel hydraulically by 
manifolds. The waterflow through the coil at 50 
lb./in.* pressure is adequate for cooling under full 
load conditions (18 kw). To improve the homo- 
geneity of the field inside the solenoid, correcting 
coils consisting of 90 turns of No. 12 wire were 
wound 12 inches from each end of the main 
winding. These correcting coils, together with a 
rheostat, are connected across the main winding. 

Careful comparison measurements with a 
standard solenoid showed that this solenoid gives 
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a field of 1.859 oersteds, ampere, corresponding to 
a field of 1115 oersteds with 600 amperes through 
the coil. Investigation of field distribution showed 
that the homogeneity of the field at all points 
between the positions of source and counter is 
better than 1 percent. 


C. Baffle system 


The construction of the baffle system is shown 
in Fig. 4. The positions and diameters of the disks 
forming the central baffle were calculated so that 
the edges of the disks lie on a surface of revolution 
formed by rotating about the axis of the instru- 
ment a sine curve of amplitude 7.16 cm and half- 
period 90 cm, starting from the source position. 
The middle diaphragms have openings whose 
edges lie on a similar surface of amplitude 12.5 cm 
and _ half-period 90 cm and _ intersecting the 
instrument axis at the same points as_ the 
previous surface. The over-all length of the 
baftle system is 82 cm, there being a space of 4 em 
between the source position and the annular 
opening in the first ring. The edges of the annular 
slit at the counter end also lie on the two surfaces 
of revolution previously described, and the slit is 
3 cm o.d. and 2 cm i.d. The resolution and 
efficiency of the spectrometer for a source of 2-cm 
diameter are determined almost entirely by this 
final annular slit. 

The central series of disks, together with their 
spacers, are held together by a Duralumin rod 
which passes through them. This central disk 
system is supported by the outer system of ring 
baffles through the two sets of radial rods. The 
baffle system was made of Duralumin to minimize 
electron scattering, and the further precaution 
was taken of painting all surfaces with thick 
colloidal graphite (Aquadag). 


D. Source holder 


The source holder is designed to minimize all 
possible scattering from the support and to 
permit quick introduction through an air lock. In 
this spectrometer, the radioactive material may 
be deposited on an area 2 cm in diameter without 
materially reducing the resolution, which permits 
the use of very thin sources. The thin collodion 
layer on which the source is deposited is mounted 
on a thin aluminum ring 3.5 cm in diameter, 
which is in turn supported by three thin alumi- 
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num rods attached by a plate to the end of a rod 
which enters the instrument through a Neoprene 
seal of the type described by Wilson and Kamen.’ 
The gate valve air lock is essential to facilitate 
quick operation without destroying the vacuum, 
since the vacuum must be maintained at about 
10°-*mm Hg by an oil diffusion pump in order to 
prevent discharges from occurring between the 
G-M counter leads. 


E. The G-M counter 


A single G-M counter is adequate for this type 
of spectrometer. Coincidence methods are un- 
necessary, since gamma-radiation from the source 
can contribute virtually nothing because of the 
large distance from source to counter and the 
interposed baffles. 

The design of the counter used in this 
spectrometer has necessarily been somewhat 
unconventional : (1) because the window through 
which the beta-rays were to enter had to be 
placed at one end of the counter rather than in 
the side wall, and (2) because it was desired to 
construct a window which should be at least 3.0 
cm in diameter and as thin as possible. The final 
design adopted is shown in Fig. 5. The cylindrical 
copper cathode is 3 cm i.d. and 4 cm long, and the 
central wire is 4-mil tungsten with a glass bead on 
the free end.'’ A glass sleeve covers the counter 
wire to within 1.5 cm of the end in order to 
minimize the active volume and reduce the 
background. The counter windows used thus far 
have been of Cellophane 0.5 mil thick (1.75 
mg cm*), but considerably thinner films of 
collodion may be used to facilitate more accurate 
measurements at low energies. The window is 
clamped between a brass ring and a 1.6 mm thick 
brass grille, a very thin rubber gasket being 
interposed between the ring and the Cellophane. 
This assembly is cemented to the front end of the 
counter with picein. (Note that the grille is on 
the outside of the Cellophane.) The grille was 
made by drilling holes 0.32 cm in diameter in a 
brass plate, and its open area is 73 percent. 

A mixture of 9 cm argon, and 1 cm ethyl 
alcohol" has proved satisfactory for operation 


*R. R. Wilson and M. D. Kamen, Phys. Rev. 54, 
1031 (1938). 

J]. R. Dunning and S. M. Skinner, Rev. Sci. Inst. 6, 
243 (1935). 

"A. Trost, Zeits. f. Physik 105, 399 .(1937). 


Fic. 5. G-M counter. 


over periods as long as several weeks. The counter 
filling system with appropriate argon and alcohol 
reservoirs is permanently connected to the 
counter, but a stopcock operable through a 
Neoprene seal permits the counter to be isolated 
from the filling system. 

While the geometrical proportions of the 
counter were necessarily bad, very satisfactory 
operation at high speeds was achieved by a 
capacitance coupled multi-vibrator quenching 
circuit which is a modification of that described 
by Getting.” An amplifier tube preceding a 
biased low impedance multi-vibrator circuit re- 
sulted in positive operation under all conditions. 
scale-of-cight demultiplying circuit utilizing 
three pairs of 6SJ7 tubes permitted counting at 
high rates. It was observed that much more 
satisfactory quenching of the counter pulses was 
obtained when the metal grille was held at about 
} the potential of the anode. 

When shielded by lead 2.5 cm thick to mini- 
mize stray radiation effects, the natural back- 
ground of the counter system was approximately 
10 counts min. The counter threshold was ap- 
proximately 1300 volts, and the ‘“‘plateau”’ about 
300 volts broad with a slight positive slope. 
Counter operation was regularly checked with a 
small radium source placed in a standard position, 
but over periods of many weeks no deviations 
were observed beyond those expected statistically. 

As shown in the drawing, a shutter operable 
from outside was placed in front of the counter. 
Four holes in the shutter permitted absorption 
measurements to be made on selected materials 
for electrons having any desired energy. 


F. Spectrometer power supply 


The current for the solenoid windings from the 
600-ampere 30-volt d.c. generator is stabilized by 
a voltage regulator similar to the type described 


27, A. Getting, Phys. Rev.{53,5103 (1938). 
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Fic. 6. Theoretical resolution curve and 
experimental points. 


by Anderson, Dunning and Mitchell," current 
fluctuations being not greater than 0.2 percent. 
Currents are measured by calibrated shunts con- 
nected to an ammeter or potentiometer. 


EXPERIMENTAL TESTS OF INSTRUMENT 


A. Test of resolving power—theoretical resolu- 
tion 

The absolute minimum and maximum values of 
mD (most convenient geometrical quantity for 
use in all discussions of resolution) can readily be 
calculated from Eq. (4), taking the source diame- 
ter as 2 cm and the baffle dimensions given 
previously. They are: tDpin=88.15 cm, tDmax 
=100.19 cm. The maximum and minimum 
starting angles for electrons are 23° 30’ and 
12° 36’, respectively. What may be called the 
“resolution curve” of the instrument is obtained 
by a calculation of the effective solid angle a, 
averaged over the entire source area, for electrons 
which pass through the baffle system as a func- 
tion of wD over the range of this quantity. The 
ordinates on the theoretical resolution curve 
(Fig. 6) are the values (obtained by a process of 


18 Anderson, Dunning, and Mitchell, Rev. Sci. Inst. 8, 
497 (1937). 
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numerical integration) of a. The maximum value 
of @ is approximately 0.122, giving a theoretical 
maximum efficiency of detection of approxi- 
mately 0.122/42%=1 percent. The value of at 
the peak of the curve is 94.7 cm. The estimated 
error in @ introduced by numerical approxi- 
mations does not exceed 1 percent in the neigh- 
borhood of the maximum of the curve and 4 
percent near the extreme ends of the curve. 

Experimental resolution.—To obtain an experi- 
mental test of the resolution of the instrument, 
the well-known beta-ray line of Th B+C+C’ 
+C” at 1p1387.5 (mean value from the measure- 
ments of Ellis'’ and Wang'*) was used. 

The source was prepared by collecting recoil 
nuclei from thoron on a collodion film 1 to 2 
microns thick. The deposit was confined to an 
area 2 cm in diameter by a diaphragm shield with 
a negatively charged collecting disk immediately 
behind the film, the holder itself being positively 
charged to minimize collection outside the 
desired source area. 

The plotted points of Fig. 6 show the results of 
the experimental resolution measurements on the 
Th B line. The statistical error in @ is approxi- 
mately 1 percent at the peak, and the continuous 
beta-ray spectrum has been allowed for by 
subtracting its contribution, which is approxi- 
mately 20 percent in the neighborhood of the 
line. 

The agreement between the theoretical resolu- 
tion curve and the experimental points (which 
have been matched at the peak) is seen to be very 
good. The total range of 7D is about 12.5 percent 
of the mean value (94.7 cm), but 80 percent of the 
area of the experimental curve lies within a band 
extending about 3 percent on each side of the 
peak. Higher resolution may readily be obtained 
by redesigning the baffle system as suggested 
previously with smaller annular slits. The present 
baffle system, however, represents a good com- 
promise between reasonably high resolution and 
large effective solid angle. 


B. Test of spectrometer efficiency 


The efficiency of the spectrometer was de- 
termined approximately by direct measurement 
in the following way. The number of a-particles 


4C. D. Ellis, Proc. Roy. Soc. 138, 318 (1932). 
18K. C. Wang, Zeits. f. Physik 87, 633 (1934). 
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emitted per second from a Th B source of the 
type described was first directly measured by 
counting, with a linear amplifier-ion chamber 
system, the number of particles emerging through 
a small foil-covered hole (of known diameter) at 
the end of a long evacuated tube which contained 
the source at a measured distance from the hole. 
The number of beta-particles emitted per second 
in the line /7p1387.5 could then be calculated by 
using the present accepted value for the co- 
efficient of internal conversion for that line 
(Pa=0.25). The number of electrons actually 
detected in the line when the source was placed in 
the spectrometer (after allowing for the opaque 
section of the brass counter grille) was then 
determined. These results gave a value for the 
effective solid angle in reasonably good agreement 
with the peak value of 0.12 shown in Fig. 6. The 
large solid angles for which the lens spectrograph 
should be effective were therefore achieved. 


C. Correction for absorption in counter window 


The absorption of electrons in the counter 
window was determined by making measure- 
ments of the transmission through 0, 1, 2, and 3 
additional layers of 0.5-mil Cellophane, when 
these were interposed in front of the counter by 
means of the shutter system described. Figure 7 
shows the correction factors by which the ob- 
served counting rates are to be multiplied, as 
obtained by extrapolation of the data for eight 
values of electron momentum to zero total ab- 
sorber thickness. These correction factors should 
be fairly reliable down to Hp 854, but are proba- 
bly not trustworthy below this value, although 
some electrons are still passing through the foils 
at 640. 
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Fic. 7. Correction factor for counter window absorption. 
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There is little evidence for scattering in the 
baffle system. With zero field, no counts above 
background are observed. 


D. Magnetic field measurements and calibration 


Measurements of the field within the solenoid 
at various points, both on and off the axis, were 
made by a flip coil which had been calibrated 
previously with a standard solenoid. These 
measurements gave for the field at the center 
of the solenoid a value of 1.830+2 percent 
oersteds ‘ampere. As stated previously, the homo- 
geneity of the field was better than 1 percent 
throughout the region between source and 
counter. 

If we take rD=94.7 cm at the peak of the 
resolution curve and use 1387.5 for the //p value 
of the Th B line, a comparison determination of 
the field may be obtained. If xD is assumed to be 
reliable to within +0.2 percent, the value of 
1.859+0.4 percent oersteds ampere is obtained. 
This is considered to be reasonable agreement 
with the previous value. 

The value 1.859 oersteds/ampere obtained 
from the Th B line has been used as a standard 
for the subsequent spectrometer measurements. 
The maximum field obtainable with the present 
power supply is thus 1115 oersteds at 600 
amperes, corresponding to a beta-ray energy of 
4.44 Mev. 
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It should be noted that in such an instrument 
the effect of the earth’s field (or other stray fields) 
cannot be completely neglected in the very low 
energy region. The axis of the instrument has 
been placed in the direetion of the horizontal 
component of the earth's field to minimize any 
effect due to this component. There is no evidence 
that the operation of the instrument is ap- 
preciably affected by the vertical component of 
the earth's field at the lowest energies so far used 
(100 kev.) However, a compensating coil for the 
vertical component is being installed for work 
down to 10 to 25 kev. 


INVESTIGATION OF BETA-SPECTRA 


A. Measurements on radium E 


The beta-spectrum of Ra E was selected for a 
general test of the spectrometer, since so many 
data are now available on it. 

The Ra E source was prepared from a purified 
Ra D+E+F solution by successive deposition 
on Ni foil followed by successive solutions in hot 
HCl to insure freedom from Ra D. The final 
highly purified Ra E solution was distributed as 
uniformly as possible over a 2-cm diameter area 
in the center of a thin collodion film mounted on 
the source holder. 
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Fic. 9. Fermi and K-U plots for radium E. 
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Figure 8 shows the data obtained with such a 
source after correcting for counter window 
absorption. The individual points have a sta- 
tistical error of 2 percent or less. Knowing the 
form of the resolution curve of the instrument, 
one can accurately determine the upper end point 
of any beta-spectrum by a process of successive 
approximations, provided that the source is 
reasonably strong. The value of the end point for 
the Ra E spectrum, after account is taken of the 
effect of the resolution curve, was computed to be 
IIp 5100 (1.11 Mev), corresponding to a total 
energy of 3.18 in moc? units. Plots of these data 
on the basis of the Fermi'® and Konopinski- 
Uhlenbeck" theories are shown in Fig. 9. Neither 
theory is seen to give a straight line over the 
complete range of the data, but the Fermi plot is 
much better. The extrapolated K-U end point is 
at a total energy of 3.83moc* units, in excellent 
agreement with the results of Flammersfeld,'® 
Neary,'® and other recent measurements.?®*! 
These tests with Ra E indicate that the instru- 
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ment can be regarded as reliable at least above 
1300 (135 kev). 


B. Beta spectrum of ,;;P* 
The source of 1,3?” was prepared by dis- 
tributing a few drops of dilute NaH PO, solution 


16 E. Fermi, Zeits. f. Physik 88, 161 (1934). 

17E. J. Konopinski and G. E. Uhlenbeck, Phys. Rev. 
48, 7 (1935). 

18 A. Flammersfeld, Zeits. f. Physik 112, 727 (1939). 

19G. J. Neary, Proc. Roy. Soc. 175, 71 (1940). 

20 [.. H. Martin and A. A. Townsend, Proc. Roy. Soc. 
A170, 190 (1939). 

2 J.S. O’Conor, Phys. Rev. 52, 303 (1937). 


] 
( 


( 
] 
12 t 
il 
a 
| 
1.0 2.0 3.0 4.0 


OVC 


BETA-RAY SPECTROMETER 41 


containing the active phosphorus over the source 
film, as in the case of Ra E. Some of the ;;P* was 
kindly supplied by Professor E. O. Lawrence and 
his associates during a period when the Columbia 
cyclotron was not in operation. In all sources, the 
activity of the phosphorus was such that the 
total weight of the deposit over an area 2 cm in 
diameter was less than 0.5 mg. 

The results of the measurements are shown in 
Fig. 10, where the statistical error of the points is 
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Fic. 11. Fermi and K-U plots for ,;P®. 


less than 1 percent. The observed end point is at 
Hp 7357 (1.75+0.02 Mev). This is in good agree- 
ment with values since published, of Lawson™ 
(1.72 Mev) and Lyman*™ (1.70+0.04 Mev). 

The Fermi and K-U plots for ;3P* are shown in 
Fig. 11. Again neither are straight lines over the 
entire range, although the Fermi plot is nearly so 
and shows less departure at the upper end. 


C. Beta-spectrum of ;;Ir'*' 


The source of j;Ir'™ was prepared through a 
concentration process of the Szilard-Chalmers™ 
type developed by Dr. J. Steigman.*® Approxi- 
mately 5 g of Ir in the form of a Werner complex 
compound of iridium and ethylene diamine were 
irradiated with neutrons from the cyclotron. 
After separation of the active Ir as Ir metal with 
approximately 1 mg Ir as a carrier, the active 
material was deposited as a colloidal suspension 
on the source film. 

# J. L. Lawson, Phys. Rev. 56, 131 (1939). 

** E. M. Lyman, Phys. Rev. 51, 1 (1937). 


* L. Szilard and T. A. Chalmers, Nature 134, 462 (1934). 
*® J. Steigman, Phys. Rev. 59, 498 (1941). 
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The spectrum obtained is shown in Fig. 12, 
where the points have a precision of at least +2 
percent. The end point is concluded to be at 
TIp 8870 (2.18+0.04 Mev), in fair agreement 
with previous estimates by other methods.** The 
curve strongly suggests that the spectrum is 
complex. If only one weak lower energy group is 
assumed, its end point would correspond to 3000 
to 3500 Hp. Measurements of the half-life in 
various parts of the spectrum gave the same 
value, 19.5 hours, to within the limits of error. 

The Fermi and K-U plots for Ir in Fig. 13 
indicate the Fermi plot to be very straight in the 
upper region. A complex spectrum is likewise 
suggested by these plots. 
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Fic. 13. Fermi and K-U plots for z7Ir™. 


26M_ S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
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A series of measurements of the absorption of 
the Ir beta-rays in Al were made using an ion 
chamber with 1-mil aluminum window filled with 
Freon (CCl.F.2) at 1 atmosphere, and which was 
connected to an FP-54 amplifier. The range ob- 
tained was 0.970+0.030 g cm*®. From Feather’s 
empirical range-energy relation, as recently cor- 
rected by Widdowson and Champion,”’ we have 
energy (Mev)=[0.165+range (g cm*) ]/0.536. 
This gives a value of 2.11+0.07 Mev for the end 
point, in good agreement with the spectrometer 
value. The absorption curve gives evidence for 
the presence of a soft gamma-radiation almost 
completely stopped by 3-mm of Pb. 

The author would like to express his sincere 


27 E. E. Widdowson and F. C. Champion, Phys. Soc. 
Proc. 50, 185 (1938). 
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The Scattering of One- to Three-Mev Protons by Helium 
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The scattering of protons of from one- to three- Mev energy by helium has been experimentally 
examined for the existence of resonance scattering analogous to the neutron-helium resonance 
observed for one-Mev neutrons by Staub and Stephens. If n-m and p-p nuclear forces are equal, 
such a resonance should occur in p-He scattering for protons of approximately two Mev. The 
number of protons scattered through 140° were observed as a function of their energy. This 
number passed through a maximum at two Mev as expected, but the sharpness and height of the 
maximum were several times less than for the n-He resonance. However, this greater width and 
lesser intensity can be qualitatively justified. The number of protons scattered through 76° did 
not show the presence of a resonance scattering as the energy was varied. Observations of the 
angular distribution of the scattered protons for angles of scattering between 30° and 140° for 
proton energies of 1.0, 1.5, 2.0, 2.5 and 3.0 Mev, respectively, have been obtained. 


INTRODUCTION 


ROTON-helium scattering experiments for 

protons of energy one- to  three-million 
electron volts are of especial interest because of 
the similarity of the proton-helium scattering 
problem to that of neutron-helium scattering. 
The dependence of the neutron-helium scattering 
cross section upon the neutron energy has 
been investigated experimentally by Staub and 


* Formerly Research Fellow of Carnegie Institution of 
Washington, now at the University of Illinois. 


Stephens.' They have found that the ratio of the 
n-He cross section for the production of recoils 
to the similar n-p cross section is a maximum for 
neutrons of about one-Mev energy. This maxi- 
mum they attribute to a resonance in the com- 
pound He® nucleus. The existence of this large 
amount of scattering of one-Mev neutrons by He 
has since been confirmed by a number of workers,” 


1H. Staub and W. E. Stephens, Phys. Rev. 55, 131-139 
939 


(1939). 
2 E. Hudspeth and H. Dunlap, Phys. Rev. 57, 971-975 
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most of whom agree that the half-width of the 
maximum is about 0.5 Mev and the height about 
a factor of five. Recently Staub and Tatel® have 
found a splitting of this resonance into two 
peaks, with an energy separation of 0.4 Mev. 
Now Primakoff and Goldsmith! and others have 
pointed out that if one assumes the equality of 
n-n and p-p forces, except for the Coulomb forces 
in the latter case, there should be a similar 
resonance in the scattering of protons by He but 
that the resonance should be nearer two than one 
Mev, the difference being due to the greater 
Coulomb energy of the compound Li® nucleus as 
compared to 

Previous proton-helium scattering experiments 
have been carried out by Roberts and Heyden- 
burg® but they had available only protons of 
one-Mev energy and the observations were con- 
sequently below the interesting region. In fact 
they were able to detect only a small departure 
from Coulomb scattering. Experiments on the 
scattering of a-particles by hydrogen yield the 
same information as those on the scattering of 
protons by helium and have been performed by 
several workers.* These have been done with 
a-particles of up to 8.5-Mev energy which corre- 
sponds to protons of up to 2.1 Mev. These 
energies are not sufficient to carry the observa- 
tions over the resonance peak, as found in our 
proton-helium scattering experiments. 

In the present experiments artificially ac- 
celerated protons of up to three-Mev energy have 
been used, so that an investigation of the entire 
predicted resonance region has been possible. 


APPARATUS 


The protons used in this experiment were 
accelerated with our pressure electrostatic gener- 
ator. A detailed description of the apparatus will 
be published in the future. The voltage was 
measured by means of a generating voltmeter 
calibrated in terms of the 0.867-Mev y-ray 
resonance of fluorine and the 1.368-Mev y-ray 


3H. Staub and H. Tatel, Phys. Rev. 58, 820-828 (1940). 


*H. Primakoff and H. H. Goldsmith, Phys. Rev. 55 
1117 (1939). 

5 N. P. Heydenburg and R. B. Roberts, Phys. Rev. 56, 
1092-1095 (1939). 

6 J]. Chadwick and E. S. Bieler, Phil. Mag. 42, 923-940 
(1921); C. B. O. Mohr and G. E. Pringle, Proc. Roy. Soc. 
160, 190-206 (1937) 
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Fic. 1. Scattering chamber. 


resonance of Al. The voltage on the electrostatic 
generator was manually monitored so that it did 
not fluctuate more than 50 kv during a run. 

The scattering chamber used in these experi- 
ments was the same as that used by Heydenburg, 
Hafstad and Tuve in one of their proton-proton 
scattering experiments and has been described in 
detail by them.’ A schematic diagram of the 
scattering chamber is given in Fig. 1. The 
magnetically analyzed proton beam from the 
electrostatic generator entered the scattering 
chamber through a thin aluminum window 
(0.0001-inch thick) and was defined by the indi- 
cated diaphragm system. The scattered protons 
entering slits Sz and S; were detected by means of 
a parallel-plate ionization chamber C;, connected 
to a linear amplifier and scale-of-eight counter. 
This chamber and attached defining diaphragms 
could be rotated from 0° to 140°. The guard 
shown in the diagram is to prevent protons 
scattered by the defining diaphragms from being 
counted. With the guard in, it was experimentally 
found that the number of counts in the absence of 
a scattering gas was negligible. 

The beam was monitored by observing with 
the ionization chamber C, the number of protons 


7N. P. Heydenburg, L. R. Hafstad and M. A. Tuve. 
Phys. Rev. 56, 1078-1091 (1939). 
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Fic. 2. Current calibration curve. 


which were scattered from the gold leaf foil 
through an angle of about 135°. This gold foil 
count, however, was calibrated in microamperes 
of incident protons by simultaneously determining 
the average rate of counting of the monitor and 
measuring the current with an RCA amplifier 
microammeter when the scattering chamber was 
evacuated. 

Most of the observations were made with 
helium of about ten-mm pressure in the scattering 
chamber. The helium pressure was measured with 
an Apiezon oil manometer. We used tank helium, 
purified by passing through a liquid-air-cooled 
charcoal trap. The helium in the chamber was 
tested for possible air impurity at the end of 
almost every run by a method which will be 
described in the next section. 


PROCEDURE 


At the beginning of the set of measurements, 
the scattering chamber was carefully aligned with 
respect to the beam direction. The generating 
voltmeter was calibrated in terms of the F and Al 
resonances at the beginning and end of the set of 
experiments. 

During the observations the beam intensity 
was measured in terms of the number of protons 
scattered from the gold foil, but in a separate 
experiment this count was itself calibrated in 
terms of the number of microamperes of current 
entering the chamber when it was evacuated. 
This calibration was made for several proton 
energies between one and three Mev. The results 
are given in Fig. 2 for two different gold foils 
(foil No. 1 was damaged during the course of the 


AND N. F. RAMSEY 

experiments and was replaced by foil No. 2). 
While these two curves are similar in shape they 
are displaced, probably because of a difference in 
thickness and uniformity of the two foils. A 
Coulomb (1 2) curve, adjusted to fit the gold 
scattering curve No. 1 at one Mev, is also plotted. 
It is seen that the observed scattering is too large 
by about a factor of two at 2.8 Mev. This ap- 
parent anomalous scattering from gold was not 
expected and is probably not real. However, as 
the effect was reproducible and was repeated 
several different times in the course of the 
measurements, the curves of Fig. 2 were taken as 
a calibration of the number of microampere- 
minutes of current into the chamber corre- 
sponding to a given gold scattering count at 
different’ voltages. That this calibration was 
satisfactory was further verified by observing the 
proton scattering by argon at an angle of 30°. 
With the above calibration, the proton-argon 
scattering was found to be strictly Coulomb to 
within the experimental error of four percent for 
proton energies from 1.2 Mev to 2.8 Mev, the 
highest voltage for which these measurements 
were made. 

As an additional check on the consistency of 
the apparatus, the angular distribution of protons 
scattered by He at one Mev was measured and 
found, in agreement with the earlier experiments 
of Heydenburg and Roberts,> to be purely 
Coulomb between 20° and 40°. 

After these preliminary and consistency tests 
had been completed the actual observations were 
made by counting the number » of protons 
scattered by the gas while another number nz 
was scattered from the gold. From these two 
numbers, from the measured voltage, from the 
gold scattering calibration, and from the geome- 
try of the apparatus, the number Ns of protons 
scattered from the chamber per microampere per 
minute could be calculated. The method of 
calculation is the same as that of Heydenburg, 
Hafstad and Tuve' and, since the same scattering 
chamber was used, the geometrical constants are 
identical. 

At the end of almost every run, the purity of 
the gas was tested by a simple scattering experi- 
ment. Since one-Mev protons scattered by He 


8M. A. Tuve, L. R. Hafstad and N. P. Heydenburg, 


Phys. Rev. 50, 806-825 (1936). 
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through 140° lost so much energy that they could 
not be counted, the only countable protons of 
such an energy scattered in this direction were 
necessarily due to the presence of a heavier 
scattering gas such as air. Hence, by determining 
the number of countable protons so scattered, the 
amount of impurity of the gas could be de- 
termined. All observations were rejected if the 
amount of impurity was sufficiently great for any 
scattering observations at an angle greater than 
39° to be appreciably affected by the Coulomb 
scattering of the impurity. 


TABLE I. Results of the scattering experiments. 


Ratio Ons. COULOMB Ratio Ns 


ENERGY Hle Au Ns PER Ns PER Oss. To Ns 
Mev Scat. X105 ua Min. ua MIN. CoUuLOMB 
76° 1.2 3.49 16.3 5.60 2.91 
1.8 9.70 24.3 2.49 9.75 
2.4 21.6 33.5 1.40 23.9 
3.0 25.1 33.9 0.896 37.8 
140 12 5.45 25.6 1.44 17.8 
1.4 10.7 40.0 1.06 37.7 
1.6 17.8 54.8 0.809 67.8 
1.8 28.4 71.1 0.639 111 
2.0 35.6 75.7 0.518 146 
2.2 39.8 71.5 0.428 167 
2.4 42.0 65.1 0.359 181 
2.6 35.3 49.5 0.306 162 
2.8 33.8 45.6 0.264 173 
3.0 30.3 40.8 0.230 177 
20 1.5 315.0 1065 1640 0.65 
40 12.53 42.4 57.9 0.73 
60 6.47 21.9 9.34 2.34 
80 8.70 29.4 2.97 9.90 
100 11.0 37.1 1.44 25.8 
120 9.90 33.5 0.974 34.4 
140 13.2 44.6 0.921 48.5 
20 2.0 236 500 923 0.54 
40 52.7 112 32.5 3.45 
60 18.7 39.7 5.26 7.55 
80 13.1 27.8 1.67 16.6 
100 14.9 31.6 0.808 39.1 
120 19.2 40.7 0.548 74.2 
140 32.7 69.4 0.502 138 
20 2.9 612 1055 590 1.79 
40 11.9 204 20.8 9.81 
60 35.9 61.9 3.36 18.4 
80 18.8 32.4 1.07 30.6 
100 16.9 29.2 0.517 56.5 
120 23.8 40.9 0.351 116.5 
140 34.7 59.8 0.331 180.7 
20 3.0 899 1215 410 2.96 
40 145.4 196.4 14.5 13.5 
60 53.8 72.7 2.34 31.1 
80 21.1 28.4 0.743 38.2 
100 16.5 22.3 0.359 62.1 
120 21.5 29.1 0.244 119 
140 45.4 0.230 197 
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Fic. 3. Dependence of the amount of scattering upon 
voltage at fixed angles (points shown by solid circles are 
taken from curves in Figs. 4, 5 and 6). 


RESULTS 


The results of these experiments are given in 
Table I. In column (4) the observed Nx has been 
obtained from the ratios in column (3) by use of 
the gold foil calibration curves of Fig. 2. The 
calculated Ns gives the number of scattered 
protons to be expected from the Coulomb 
scattering law of Rutherford and Darwin for our 
scattering chamber dimensions. These results are 
plotted in Figs. 3, 4, 5 and 6, both the experi- 
mentally observed scattering and the theoretical 
Coulomb scattering being plotted. The calculated 
Coulomb scattering values given in column (5) 
and plotted in the figures have been multiplied by 
a (1/sin@) factor to take into account the change 
in the effective scattering volume with angle.* 

Measurements were taken both at fixed angles 
as the voltage was varied and at fixed voltages as 
the scattering angle was varied. The results of the 
former type are shown in Fig. 3. These curves are 
at 140° and 76° in the laboratory system or at 
149° and 90° in the center-of-gravity system. 
These angles are such that a pure p wave should 
show 75 percent of its maximum effect at 149° 
and none at 90°. The curves at fixed voltages as 
the scattering angle was varied are shown in 
Figs. 4, 5 and 6. Although these two sets of 
observations were made independently, they 
obviously overlap. Therefore, points obtained 
from Figs. 4, 5 and 6 have also been plotted on 
Fig. 3 to make clear the extent of the agreement 
between the two types of observations. 


DISCUSSION OF RESULTS 


From Fig. 3 it appears that the amount of 
scattering does pass through a maximum at 
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Fic. 4. (above) Angular distribution of 1.5-Mev protons 
scattered by helium. 


Fic. 5. (below) Angular distribution of 2.0-Mev protons 
scattered by helium. 


about two Mev as expected theoretically on the 
basis of the neutron-helium scattering results and 
the assumption of the equality of m-n and p-p 
forces. However, the height and narrowness of 
the maximum are much less than in the neutron- 
helium case, being only a factor of 2 instead of 
5 in intensity, while the half-width is more than 
one Mev instead of about 0.5 Mev. Nevertheless, 
as Dr. Fano, Dr. Breit, and others have indicated, 
it is quite reasonable that the Li> resonance 
should be broader and less intense than the He®. 

This increased breadth of the resonance makes 
the interpretation of these results more difficult. 
In the first place it is more difficult to be certain 
that the maximum of the curve really is a 
resonance, since the amount of scattering could 
pass through a slight maximum merely due to the 
fact that the scattering cross section for a given 
type of scattering decreases with the increasing 
energy of the scattered particles. However, the 
maximum is probably sufficiently intense to 
justify at least the tentative assumption that the 
observed maximum is due to a resonance. 

A second difficulty arising from the increased 
breadth of the curve is the difficulty in de- 
termining if there is a doublet structure to the 
curves. In the neutron-helium scattering experi- 
ments there has been disagreement among the 
observers as to whether there is a single or 
doublet structure of the resonance.'? However, 
the more recent work of Staub and Tatel® seems 
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to show conclusively that the resonance is a 
doublet. They have analyzed their data in terms 
of the theoretical equations given by Bloch® and 
find that their resonance curve can be fitted by 
assuming the P32 and P,2 levels to be either an 
inverted or a normal level system depending on 
the values assigned to the constants occurring in 
the theoretical equation. They prefer the latter as 
being the correct interpretation when other 
considerations are taken into account. As proton 
experiments are in general less difficult than 
neutron ones, it was hoped for a time that it 
would be easier to determine the existence of a 
doublet structure in proton-helium scattering 
experiments than in 2-He experiments. However, 
the greater width of the maximum in the p-He 
scattering makes it more, rather than less, 
difficult to detect a slight doublet structure in the 
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Fic. 6. Angular distribution of protons scattered by 
helium. Curve A, 2.5-Mev protons; curve B, 3.0-Mev 
protons. 


curve. However, as far as our results yield 
information on the multiplicity of the resonance, 
there is no evidence for the existence of a 
resolvable doublet structure to the maximum in 
proton-helium scattering. 

The authors wish to thank Dr. J. A. Fleming, 
Dr. M. A. Tuve, Dr. L. R. Hafstad, R. C. Meyer, 
and others of the Department of Terrestrial 
Magnetism staff for their friendly cooperation 
and assistance in this work. 
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The Acceleration of Electrons by Magnetic Induction 


D. W. Kerst* 
University of Illinois, Urbana, Illinois 
(Received April 18, 1941) 


Apparatus with which electrons have been accelerated to an energy of 2.3 Mev by means of 
the electric field accompanying a changing magnetic field is described. Stable circular orbits are 
formed in a magnetic field, and the changing flux within the orbits accelerates the electrons. As 
the magnetic field reaches its peak value, saturation of the iron supplying flux through the orbit 
causes the electrons to spiral inward toward a tungsten target. The x-rays produced have an 
intensity approximately equal to that of the gamma-rays from one gram of radium; and, be- 
cause of the tendency of the x-rays to proceed in the direction of the electrons, a pronounced 


beam is formed 


N the past the acceleration of electrons to very 

high voltage has required the generation of 
the full voltage and the application of that 
voltage to an accelerating tube containing the 
electron beam. No convenient method for re- 
peated acceleration through a small potential 
difference has been available for electrons, al- 
though the method has been highly successful in 
the cyclotron for the heavier positive ions at 
velocities much less than the velocity of light. 

Several investigators'~* have considered the 
possibility of using the electric field associated 
with a time-varying magnetic field as an ac- 
celerating force. This is a very attractive possi- 
bility because the magnetic field can be used to 
cause a circular or spiral orbit for the electron 
while the magnetic flux within the orbit increases 
and causes a tangential electric field along the 
orbit. The energy gained by the electron in one 
revolution is about equal to the instantaneous 
voltage induced in one turn of a wire placed at 
the position of the orbit. Since the electron can 
make many revolutions in a short time, it can 
gain much energy. The comparatively small 
momentum of a high energy electron requires 
correspondingly small values of Hr for high 
energy orbits. For example, the energy of an 
electron when y~c is KE=3X10-4Hr—0.51 mil- 
lion electron volt. Thus with /7= 3000 oersteds 


*On leave at the General Electric Company Research 
Laboratory. 

'G,. Breit and M. A. Tuve, Carnegie Institution Year 
Book (1927-28) No. 27, p. 209. 

? R. Widerée, Arch. f. Electrotechnik 21, 400 (1928). 
(1925) T. S. Walton, Proc. Camb. Phil. Soc. 25, 469-81 

*W. W. Jassinsky, Arch. f. Electrotechnik 30, 500 
(1936). 


and r=5 cm, the energy of the electron would be 
about 4 Mev, and the orbit could be held between 
the poles of a smafl magnet. 

Because of the experimental experiences of pre- 
vious investigators'* with this method of ac- 
celeration, a rather detailed study of the focusing 
to be expected was made, and it is presented in 
the paper immediately following this one. With 
the results of this theoretical investigation to 
guide the design, it was possible to make an 
induction accelerator which produced x-rays of 
2.3 Mev.*® Briefly, in the focusing theory it is 
shown that: 

1. The electrons have a stable orbit, “‘equilib- 
rium orbit,’ where 


= (1) 


¢o is the flux within the orbit at ro, and //, is the 
magnetic field at ro. Both @» and // are increased 
during the acceleration process. This flux con- 
dition holds for all velocities of the electrons, and 
it shows that if a maximum flux density of 10,000 
gauss is allowed in the iron then 5000 oersteds is 
the maximum magnetic field which can be used 
at the orbit. 

2. In the plane of their orbits the electrons 
oscillate about their instantaneous circles, circles 
for which p =el/r/c, with an increasing frequency 


w,=2(1—n)}, (2) 


where @ is the angular velocity of the electron in 
its orbit, and w, is 2 times the radial focusing 
frequency. The number n is determined by the 


5D. W. Kerst, Phys. Rev. 58, 841 (1940). 
® 1D. W. Kerst, Phys. Rev. 59, 110 (1941). 
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radial dependence of the magnetic field, which we 
take to be of the form J7~1. r". For radial focus- 
ing m must be less than unity. 

3. Axial oscillations, oscillations perpendicular 
to the plane of the orbit, have 


(3) 


For axial stability 2 must be greater than zero. If 
the beam is to be smaller in an axial direction 
than it is in a radial direction then n>}. 

4. Decrease of the amplitude of both axial and 
radial vibration occurs because of the increase of 
the restoring force with increasing magnetic field. 
At nonrelativistic velocities the damping is 


da/a= —dE 4E, (4) 


where dE/E is the fractional increase of the 
kinetic, energy of an electron and da a is the 
fractional decrease in amplitude of the oscillation 
about the instantaneous circle. This holds for 
both axial and radial oscillations. 

5. Instantaneous circles not coincident with 
the equilibrium orbit shrink or expand toward 
coincidence : 


Qn 


dx x= —dE 2E, (5) 


where x is the displacement of the instantaneous 
circle from the equilibrium orbit and dx is the 
shift of the circle toward the equilibrium orbit 
while the electron’s energy increases by the 
fraction dE. E. 

Because of the shrinking or expansion of the 
instantaneous circle of a displaced electron 
toward the equilibrium orbit and the decrease of 
the amplitude of oscillation of an electron about 
its instantaneous circle, it was expected that a 
cathode or an electron injector placed outside of 
the equilibrium orbit could shoot in electrons 
which would miss the injector on successive 
revolutions around the magnet. Furthermore, 
since an instantaneous circle for electrons with a 
constant injection energy exists within the ac- 
celeration chamber for a small but finite interval 
of time while the magnetic field is increasing, a 
finite amount of charge should be captured in 
orbits not striking the walls. With constant 
potential on an injector the electrons first would 
hit the outer wall of the chamber before the 
magnetic field had grown large enough to give 
the electrons a radius of curvature less than the 
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radius of the wall. Then as the field increased it 
would reach values which give instantaneous 
circles within the vacuum tube. Eventually the 
field would be so large that the electrons coming 
out of the injector strike the walls or spiral 
around in small circles. 

From Eq. (2) it can be seen that the spreading 
rays from the injector will form injector images 
m(1—mn)~! radian apart, since if each ray of the 
beam oscillates about the instantaneous circle, it 
also oscillates about the central ray from the 
injector. 

Equations (4) and (5) indicate that for large da 
or dx it is desirable to use a small injection 
voltage -, and a high voltage per turn dE. A 
shift da in two revolutions or dx in one revolution 
of about 1 mm when the displacement of the 
injector from the equilibrium orbit is x=a=1.5 
em would require dE E=}. It is not possible to 
decrease the injection voltage E indefinitely 
since scattering increases when this is done, and 
it may be that the beam would be lost by scat- 
tering through an angle greater than that which 
magnetic focusing can handle. Conservative esti- 
mates on the scattering out of a cone of 7° half- 
angle, which is about the focusing limit, showed 
that at 10-* mm of Hg air pressure in the tube, a 
maximum //r=15,000 gauss cm, and f=600 
cycles sec. for the frequency of oscillation of the 
magnetic field, the injection voltage must exceed 
100 volts for less than a 20-percent loss of beam. 
This then requires that dE or the voltage per turn 
at the equilibrium orbit must be at least 12 volts 
at the time of injection. In the magnet which was 
built this condition was easily satisfied, for 25 
volts per revolution could be attained. There is 
no trouble with energy loss by the electrons in 
their long paths through the residual gas in the 
vacuum tube. Wideréc? was aware of the flux 
requirement (1), but apparently he did not 
realize that the voltage gain per revolution could 
not be too small compared with the injection 
voltage, or that it would be necessary to inject 
the electrons nearly tangent to the equilibrium 
orbit. 

After the injected electrons have orbits ap- 
proximately coincident with the equilibrium 
orbit, they should remain close to 7». A simple 
method was used to shift the equilibrium orbit so 
that the electrons struck a thin tungsten target. 
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\ portion of the flux through the center of the 
orbit passed through disks cemented to the center 
of the pole faces and made of material which 
saturated easily. As saturation progressed, the 
electrons had to move inward to a smaller radius 
of curvature in order that ¢=272r/7/. Eventually 
the electrons grazed past the target producing 
thin target x-radiation. 


SPACE CHARGE EFFECTS 


It is difficult to estimate the effect of space 
charge within the beam on the formation of 
images, since there are two focusing oscillations, 
axial and radial, in general having different fre- 
quencies. However, an estimate can be made of 
the upper limit of current which the magnetic 
forces can hold. 

At injection the magnetic forces are small so 
that the estimate should be made with the mag- 
netic field at this time. The electrons at the edge 
of the beam, which is assumed to be a cylinder of 
radius 4, determined by the distance between the 
equilibrium orbit and the injector, will experience 
a repulsive force of 


F, = (cQe), (10mro4) (6) 


dynes. Q is the number of coulombs in the orbit, 
and 7 is the radius of the orbit. 
For radial unbalanced magnetic forces we get 


me? (7) 


dynes acting on the electrons toward ry. If 
F.= F,, at the edge of the beam, then this balance 
will hold within the beam for all other distances 6 
smaller than A, because the Q of (6) is pro- 
portional to 6. Thus 


(8) 
Or in terms of the injection voltage £, 
(15rec). (9) 


/ is the target current in amperes, f the frequency 
of oscillation of the magnetic field, and n gives 
the dependence of the magnetic field on r. Since 
the value of E which can be used depends upon 
the voltage per turn dE, which is also dependent 
on the peak magnetic field reached or the output 
energy, and also upon the orbit shift da which is 
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required to miss the injector, we get 
(10) 


where KE is the final energy of the electrons in 
Mev. Inserting into (10) the constants of the 
machine which was built, we have A=1.5 cm, 
f=600 cycles/sec., KE=2.3 Mev, and da=0.1 
cm. The target current should be about 0.03 
microampere. This is fairly close to the lower 
limit of current as estimated from the x-ray 
output of the accelerator. This lower limit was 
found by using the thick target yield of Van Atta 
and Northrup’ in the direction of the electron 
beam for comparison with the x-ray yield from 
the induction accelerator. It is found that the 
target current in the induction accelerator must 
be greater than 0.02 microampere. It -is only 
possible to find the lower limit of this current 
since the output is thin target radiation. 

A curious behavior which was immediately 
noticed when the apparatus began to work was 
that injection voltages much greater than the 
largest voltage which would allow the orbit to 
miss the injector could be used. In fact the yield 
increased greatly as this voltage was raised. This 
seems understandable on the basis of space charge 
forces spreading the beam away from the central 
ray of the injector. 

At relativistic energies space charge forces are 
completely balanced by magnetic self-focusing of 
the beam, for the electric force on a stray electron 
at a distance A from the beam center is 


eS=2ee A, (11) 


where o is the linear charge density in e.s.u.,cm. 
The magnetic attraction due to the main current 
in the beam is 


c=(v c)*2ee A. (12) 


Thus it is evident that when vc, the magnetic 
pull of the beam for a stray electron just equals 
the electrostatic repulsion. Or, from the point of 
view of an observer on the electron, the spacing 
of the fixed number of electrons around the orbit 
will increase, since as vc his yardstick becomes a 
smaller fraction of the circumference of the orbit. 


7L.C. Van Atta and D. L. Northrup, Am. J. Roentgen. 
Rad. Ther. 41, 633 (1939). 
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THE MAGNET 


A photograph of the apparatus is shown in 
Fig. 1. 

For the production of the rapidly changing 
magnetic field, with the proper dependence on 
the radius, finely laminated iron pole faces were 
made from 0.003-inch silicon steel sheets. The 
return magnetic circuit made from the same 
material was interleaved for mechanical strength ; 
and the roughly circular central pole pieces were 


Fic. 1. The induction accelerator. The glass doughnut 
is between pole faces which are held apart by eccentric 
wedges. 


formed by stacking with different widths of 
laminations as shown in Fig. 2. Each pole piece 
was capped by a disk of radially arranged lami- 
nations so that perfect circular symmetry was 
achieved at the pole surfaces (B, Fig. 2). The 
whole pole face was held together by a thick 
Transite or asbestos board ring about its perime- 
ter, with cement of water glass and flint dust 
filling the cracks between the laminations and 
hardened by baking. The pole caps were held 
against the pole pieces by eccentric wedges be- 
tween the Transite rings. 

To supply the flux within the orbit which was 
necessary to hold the electrons out at the 
equilibrium orbit, two disks about two inches in 
diameter, each made of pressed iron dust of 
permeability about eight, were cemented onto 
the flat center portions of the pole faces. The 
thickness of these disks was chosen so that the 
equilibrium orbit was formed at about 7.5 centi- 
meters radius, and the final adjustment of the 
position of the equilibrium orbit was made by 
painting a mixture of water glass and iron filings, 
or the turnings from the compressed iron dust 
disks, onto the surface of the iron dust disks so 
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that the reluctance of the gap was correct. Since 
the iron particles in these disks are separated, the 
flux density through them is greater than the 
average flux density through the disks. This 
causes the saturation which shrinks the equilib- 
rium orbit down to the target. 

Because of the large leakage flux the main coil 
of the magnet had to be highly subdivided. Two 
hundred strands of No. 20 double enameled wire 
were made into a cable approximately 10 yards 
long which was twisted about ten times and 
which had its inside wires interchanged with its 
outside wires several times. This was formed into 
a 10-turn coil with about ,'5 inch of insulation 
between turns. Two such 10-turn coils were made 
and fitted with large lugs of 1-inch copper tubing. 
The coils were wrapped with cotton tape, dipped 
in Bakelite varnish and baked. Figure 3 shows 
the circuit with the coils connected to a total of 
eighty 5-microfarad Pyranol condensers which 
were rated at 660 volts a.c. Energy was supplied 
to this resonating circuit by a 2-turn primary of 
stranded enameled wire around the pole pieces. 
The r.m.s. electromotive force in this primary 


+ 


| 
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Fic. 2. Dimensions of the magnet. Parts B are the pole 
caps made with all laminations placed in a radial direction. 
The iron dust disks supply the central flux. Outside of 
the 14° conical surface a flat rim on the pole face tends 
to prevent too rapid decrease of the field at the edge of 
the gap. 


14° 2 > lo t 


could be run up above 100 volts, but 60 to 80 
volts were usually used in operation since this 
was sufficient to saturate the powdered iron and 
to collapse the orbit. Power was supplied by a 
4-kilowatt 600-cycle alternator driven by a 
direct-current motor with an adjustable speed. 
To determine how the magnetic field decreased 
with 7 a small search coil was arranged so that it 
could be held between the poles at different 
radial distances from the center. The e.m.f. from 
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this search coil was bucked against the e.m.f. 
from a voltage divider connected across a 1-turn 
coil about the leg of the magnet, equality of 
e.m.f. being determined by no deflection of the 
beam of an oscillograph. Since the wave form of 
the voltage on the search coil differed slightly 
from the wave form from the coil around the core, 
balance could be obtained only at one phase, and 
the phase of interest was that of zero magnetic 
field. The use of the oscillograph as a_ null 
instrument in this way proved to be sufficiently 
accurate. The pole faces gave a field following the 
1 r} law between a radius of 4.5 centimeters and 
9.25 centimeters when the separation between 
the flat central portions of the pole faces was 2.8 
centimeters. At r=10 em the law was 1 7. 

To determine the position of the equilibrium 
orbit use was made of a concentric system of 
seven 1-turn coils set in grooves in a Bakelite disk 
which fit snugly around one of the iron dust disks. 
Relative values for the e.m.f.’s induced in these 
coils were determined by the same null method 
which was used for determining the shape of the 
field. These data could be used in several ways to 
find ro. The simplest method is to find the radius 
of minimum electric field, since this will be 7o. 


THE ACCELERATION CHAMBER 


It was desirable to have the volume available 
for electron orbits as large as possible. This 
required a doughnut-shaped glass vessel with 
walls parallel to the conical pole faces. The outer 
diameter was 20 cm, and the center of the 
doughnut was 7 cm in diameter. The glass work 
for this vessel required great skill.6 A 20-cm 
spherical bulb with a wall thickness of 2.5 mm 
was flattened and dished conically on both sides 
so that it would fit between the magnet poles. 
The central 7-cm hole of the doughnut was 
formed by pushing the centers of the dished faces 
together and picking out the glass. 

The inside of this bulb was coated with a thin 
conducting layer by chemical silvering. This 
coating is necessary to prevent stray charges 
from building up potentials on the walls. Contact 
is made to the silver coating which is then 

*I am indebted to the Vacuum Tube Department of 


the General Electric X-Ray Corporation for producing a 
correctly shaped tube. 
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Fic. 3. The resonant circuit which energizes the magnet. 
Losses are supplied by a 4-kilowatt 600-cycle-per-second 
generator. 


grounded outside of the tube. Resistances of the 
silver coat between 20 ohms and 300 ohms 
between test probes put in the side arms were 
used. 

Both target and injector were mounted on the 
glass pinch seal which was waxed to a flare on the 
doughnut. The target was a piece of 0.015-inch 
tungsten sheet folded to present an edge of large 
radius of curvature to the electron orbit which 
shrinks inward toward it. Thin tungsten is used 
so that eddy currents will not be excessive. The 
injector was made of thin molybdenum sheets 
and a cylindrically spiraled tungsten wire for the 
filament. The whole assembly is shown in Fig. 4. 

Fortunately the presence of this much metal 
near the orbit does not seem to disturb the local 
magnetic field too much. However, the orbits are 
destroyed if a piece of metal is brought up to the 
side of the glass doughnut while the machine is 
running. This was observed when a block of 
beryllium was put near the tube for experiments 
with photo-disintegration. The x-ray yield disap- 
peared when the beryllium was too close. The 
original injector was arranged to send a beam of 
electrons in both directions around the magnet so 
that both directions of field could be used. Thus 
two beams of x-rays oppositely directed were 


produced. 


OPERATION 


Although the focusing theory shows that for 
the present design of the accelerator electrons 
should be injected at approximately 200 volts, it 
was found that the vield increased with voltage 
and that it was still increasing at 600 volts. 
Likewise the negative voltage on the small 
focusing plates G beside the filament gave an 
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Fic. 4. Cross section of the doughnut-shaped acceleration chamber. The equilibrium 
orbit is at ro, T is the tungsten target, A is the injector, B is a top view of the injector. 
A ribbon beam of electrons from the filament F is shot out through slots in the positive 
plates P. G are negative focusing electrodes. 


increasing vield as it was raised to 200 volts. The 
beam from the injector was allowed to pour into 
the vacuum tube continuously. 

Initially the vield obtained was approximately 
equivalent to 10 millicuries of radium when 
measured in a direction from the target at right 
angles to the electron beam as it strikes the 
target, and the yield was approximately ten 
times greater in the direction of the electron 
beam. Later the vield in the direction of the beam 
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Fic. 5. The relative intensity of x-rays from the induction 
accelerator as a function of absorber thickness. A is 
absorption in copper. The absorption coefficient is 0.454 
cm! which gives a monochromatic equivalent of 1.35 Mev 
for the x-rays. B is the curve for lead. The coefficient is 
0.62 cm™', which corresponds to 1.4 Mev monochromatic 
equivalent. C is the lead absorption curve of Van Atta 
and Northrup for x-rays produced by 2.0-Mev electrons 
striking a thick target. 


was increased to the equivalent of one gram of 
radium. 

Figure 5 shows an absorption curve of the 
radiation as it passes through lead, and the 
comparison with two-million-volt x-rays from the 
M. I. T. electrostatic machine’ shows that the 
energy reached in the induction accelerator is 
approximately 2.0 million electron volts and that 
the radiation produced is thin target radiation. 
With other iron dust disks cemented on the pole 
faces, absorption coefficients indicating 2.3-Mev 
x-rays have been obtained. By means of a colli- 
mated Geiger-Miiller counter the source of the 
X-rays was shown to be the target. A large lead 
block witha small hole initlet anabundantamount 
of x-rays through to the Geiger-Miiller tube only 
when it was pointed at the target and not when it 
was pointed at the injector or other portions of 
the vacuum tube. By connecting the vertical 
deflection plates of an oscillograph to the Geiger- 
Miller counter and the horizontal deflection 
plates to a 1-turn coil around one leg of the 
magnet, the phase of the magnetic field at which 
electrons struck the target could be determined. 
As was expected it was found that the greater the 
excitation of the dynamo which energizes the 
magnet, the earlier in the cycle the iron dust 
saturated and brought the beam in to the target. 
Lowering the primary voltage postponed suffi- 
cient saturation in the iron dust until the peak of 
the magnetic field was reached. The Geiger- 
Miiller counter then gave x-ray pulses at the 
center of the oscillograph screen. This indicated 
a path length of about sixty miles from injector 
to target. If the primary voltage was lowered 
beyond this point, the yield disappeared, for the 
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electrons were not drawn in to the target but 
were slowed down by the decreasing magnetic 
field. Fortunately the operation of the accelerator 
is not sensitive to the alignment of the pole faces. 
No difference in the output can be detected when 
the pole faces are placed off axis as far as a 
thirty-second of an inch. It is also surprising that 
vacuum requirements are not as severe as was 
expected. No rigorous outgassing is necessary and 
the apparatus has been run with a vacuum as 
poor as 10-' mm Hg. The tube can be opened for 
changes and operated three-quarters of an hour 
after sealing shut. 

At present, low flux densities have been used at 
the orbit. When these are increased, it should be 
possible to go to 5 million volts even with this 
small model. One of the promising possibilities 
for the induction accelerator as a research tool is 
that the electrons from the beam can come out 


through the glass walls of the doughnut after 
they strike the target. They should be fairly 
homogeneous in energy provided that the target 
has a high atomic number. The great increase in 
bremsstrahlung production with rising electron 
energy in addition to the concentration of this 
radiation in a cone of solid angle mc?/ £ about the 
original electron direction gives the induction 
accelerator the possibility of providing an intense 
source of x-radiation for nuclear investigations. 
Since there is no evident limit on the energy 
which can be reached by induction acceleration, 
it may soon be possible to produce some small 
scale cosmic-ray phenomena in the laboratory. 

I am indebted to Professor H. M. Mott-Smith 
and Professor R. Serber for many discussions of 
the theoretical aspects of this problem and to 
Mr. R. P. Jones for assistance in the construction 
of the magnet. 
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The first section gives a general account of the principles of operation of the electron induction 
accelerator. The second section gives the more detailed analysis of the orbits of the electrons 
which was undertaken to serve as a guide in the design of the accelerator. 


I. INTRODUCTION 


HE construction and operation of an 

induction accelerator for electrons has been 
discussed in the preceding paper.' The idea of 
using the principle of electromagnetic induction 
for the production of high energy particles has 
been entertained by a number of investigators,’ 
but has hitherto met with little success in 
application. It was therefore felt necessary to 


*On leave at the General Electric Company Research 
Laboratory. 

‘Also D. W. Kerst, Phys. Rev. 58, 841 (1940); ibid. 
59, 110 (1941). 

?G. Breit and M. A. Tuve, Carnegie Institution Year 
Book (1927-28) No. 27, 209. R. Widerée, Arch. f. Electro- 
technik 21, 400 (1928); E. T. S. Walton, Proc. Camb. Phil. 
Soc. 25, 469-81 (1929); W. W. Jassinsky, Arch. f. Electro- 
technik 30, 500 (1936). 


carry out a more careful analysis of the orbits of 
electrons in changing magnetic fields for the 
double purpose of determining whether such a 
device was practicable, and to serve as a guide 
in its design. 

The basic idea of the accelerator is a simple 
one. An electron started in a radially symmetric 
magnetic field at the proper position and with 
the right velocity will move in a circle of radius 
given by 

p=ellr_c. (1) 


If now the magnetic flux enclosed by the orbit, 
¢, is increased, a tangential electric field will be 
produced at the orbit, which will 
accelerate the electron. If the magnetic field is 
so arranged that p and J// increase proportion- 
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ately, the radius of the orbit will remain un- 
changed; the electron continues to move in the 
equilibrium orbit (r=7o), but with momentum 
constantly increasing as // is increased. 

The condition for the existence of such an 
equilibrium orbit is readily found. The rate of 
increase of momentum is 


p=ed 
which gives 
2mroc. 


Using (1), with r=ro, we find 
(2) 


Thus the change in flux through the orbit must 
be twice that which would obtain if the magnetic 
field were uniform in space. To satisfy this 
requirement one must have a strong central 
field to supply the necessary flux, and a weaker 
field to hold the electron in its circular orbit. 

The energies which can be obtained by this 
means are quite high. The Illinois accelerator 
(ro=7.5 em, Hmax=1200 gauss) should give, 
according to (1), electrons of 2.2 Mev. It will 
be observed that, since the phase of the electron 
in its orbit is immaterial to the operation of the 
accelerator, the relativistic change of mass with 
velocity causes no difficulty for a machine of 
this type. 

However, the very fact that such high veloci- 
ties are attained brings with it new difficulties 
which at first sight seem rather formidable. 
Since, during most of the time the electron is 
being accelerated, its velocity is very nearly the 
velocity of light, its path length in the machine 
is very great. The magnet of the Illinois machine 
was activated by 600-cycle a.c., and the acceler- 
ation took place during a quarter-cycle. In this 
time the electrons traveled nearly 100 kilometers 
and made about 200,000 revolutions. It is 
therefore essential, to obtain any appreciable 
beam intensity and to overcome the effects of 
scattering by the residual gas and of space 
charge, that there be strong focusing forces to 
hold the electron in its equilibrium orbit. 

The focusing requirements impose additional 
restrictions on the form of the magnetic field. 
Consider first the radial focusing, that is, the 
motion in the plane of the orbit. The equilibrium 
Eq. (1) is just the condition for the balancing of 


Il = (¢— $1), 
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the magnetic and centrifugal forces. If the 
electron is displaced from the equilibrium radius, 
an unbalanced force will act, and, since the 
centrifugal force is proportional to 1/7, this force 
will be directed towards the equilibrium radius 
provided the magnetic field falls off less rapidly 


than 1/7. The electron, when displaced, will 


then oscillate around the equilibrium orbit. We 
shall see that this is a damped oscillation ;* its _ 
amplitude is proportional to /7/~?. The focusing 
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Fic. 1. The figures represent the developed paths of 
electrons. 7o is the position of the equilibrium orbit, and r; 
is the position of the instantaneous circle. The injector is 
at 0. (A) Path of an electron injected tangentially on its 
instantaneous circle. It approaches ro without oscillation. 
(B) Path of an electron with an instantaneous circle 
coincident with the equilibrium orbit. Oscillation occurs 
about ro. (C) Path of an electron whose instantaneous 
circle does not pass through 0 at the time of injection but 
is between 0 and 7. The oscillation is about 7; while r; 
approaches ro. (D) A real beam from the injector showing 
image formation at 1, 2, 3 and 4. The instantaneous circle 
is coincident with the equilibrium orbit. 


is thus most effective at small velocities, where 
scattering and space charge effects are most 
important. Figure 1(B) shows such a damped 
oscillation about the equilibrium radius. There 
will also be vertical focusing, forces tending to 
drive the electron back to the plane of the 
orbit, provided the magnetic field strength is 
decreasing with increasing 7. The restoring force 
here is due to the curvature of the magnetic 
field near the median plane; the focusing action 
is similar to the magnetic focusing effective in 


3 This can be shown in the following way. The energy 
of oscillation is E=Jv,, where J is the action variable, 
v, the frequency of oscillation. The amplitude of oscillation 
a is given by E=}M(2zv,)’a?, with M the transverse 
mass. Thus a?=J/2x*Mv,. The frequency of rotation of 
the electron in its orbit is y= p/2x7M, so, writing n;=v,/», 
a?=Jr/xnip. It follows from the adiabatic invariance of 
theaction variable that For H~1/r", 
n; turns out to be a constant. 
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the cyclotron. The vertical oscillations are 
similarly damped. Thus both radial and vertical 
focusing can be obtained, for a magnetic field 
which in the region of the orbit is of the form 
I]~1,/r", within the limits <1. 

The injection of electrons into the machine is 
a problem which also has its difficulties. It is 
obviously not possible to inject the electrons at 
the position of the equilibrium orbit. Nor is it 
possible to shoot them into the magnetic field 
from an external source as was attempted by 
Widerée; for, unless the magnetic field were 
increasing extremely rapidly, the orbits would 
be nearly symmetric around the point where 
the radial velocity vanished: the electrons would 
shoot right out again. Thus the electrons must 
be injected at a point within the magnetic field, 
either inside or outside the equilibrium orbit 
(in the following we shall suppose it outside). 

The analysis of the paths of the electrons is 
greatly simplified by the fact that the magnetic 
field increases little during one revolution of the 
electron. The adiabatic theorem can therefore 
be applied: it states that under these circum- 
stances the motion is very nearly the same at 
each instant as if the magnetic field were held 
fixed at its instantaneous value. Consider now 


_an electron injected tangentially with a given 


momentum p, at radius 7. As the magnetic field 
is increased from zero there will come a time 
when p, H, and r satisfy the relation (1). An 
electron injected at this time would, if the 
magnetic field remained constant, travel in a 
circle. The adiabatic theorem tells us that at 
any later time, even though the field is in- 
creasing, the position of the electron will still be 
given by (1); at each instant the electron is 
moving on its ‘instantaneous circle.’ Since the 
injection radius is larger than the radius of the 
equilibrium orbit, the flux within the instan- 
taneous circle is less than that necessary to 
hold the electron at r. Momentum is gained less 
rapidly than necessary to keep the radius of the 
orbit constant, and the instantaneous circle 
shrinks towards the equilibrium orbit. The 
actual orbit of the electron will thus be a spiral 
which approaches the equilibrium orbit asymp- 
totically as shown in Fig. 1(A). 

This contraction of the instantaneous orbit in 
itself provides relatively little clearance for the 


injection electrodes. The contraction is fastest 
for small velocities, but it is inadvisable to use 
too small injection energies because of relatively 
large space charge forces and scattering. In the 
University of Illinois accelerator (at injection 
energy 180 v, injection radius 9 cm, 7>=7.5 cm) 
the instantaneous circle shrinks 1 mm during 
the first revolution of the electron. However, 
electrons injected not quite tangentially, or a 
little too early or too late, so that at the time 
of injection their instantaneous circle lies outside 
or inside the injection radius, will, as the dis- 
cussion of focusing shows, oscillate with de- 
creasing amplitude about the appropriate spiral 
orbit. This situation is shown by Fig. 1(C). The 
period of these oscillations is different from the 
period of revolution, so that the electron will, 
in general, execute several decreasing oscillations 
before the maximum displacement occurs near 
the starting electrodes. There is thus a consider- 
able range of starting angles and times which 
allow capture into the equilibrium orbit, even 
though the electrons are injected at constant 
energy. In Fig. 1(D) the paths of divergent 
rays from an injector are shown. The formation 
of an image occurs at intervals of half an oscil- 
lation. In addition to these magnetic effects, 
the space charge spreading of the beam can also 
play a role in causing electrons to clear the 
injector. 

After the electrons have been accelerated a 
variety of devices can be employed to bring 
them to strike a target, all of them based, of 
course, on destroying the relation (2) at the 
appropriate time. A simple method was used in 
the Illinois accelerator to accomplish this. A 
small amount of iron dust in an insulating 
cement made up a part of the magnetic circuit 
which fed flux through the center of the orbit. 
These dust particles saturate before the iron 
core does. When saturation sets in, p increases 
less rapidly than the magnetic field at the orbit. 
The electrons then spiral inward until they hit 
the target. 


II. THe ELECTRONIC ORBITs 


The position of the electron will be specified 
by the cylindrical coordinates (r, g, 2). We shall 
first consider the motion in the median plane 


Mis 
s of 
dr; 
ris 
its 
ion. 
rcle 
‘urs 
but 
‘ing 
rcle 
pre 
ost 
ed 
-re 
to 
he 
is 
‘ce 
tic 
on 
in 
gy 
le, 
on 
of 
/y, 
of 


56 


(s=Q0), in which the magnetic field has only a 
z component, //,=//(r, This field can be 
derived from a vector potential, A,=A(r, 4), by 
the relation 


I,=r'0(rA), dr. (3) 


There is a concomitant electric field —A 
The Hamiltonian function is 


rteA c. 


Here py is the canonical angular momentum, 
defined by the relation 


(1-7, c)}, 


and is a constant of the motion. 

For any py, the radius 7; of the instantaneous 
circle is determined by the condition p,=0, or, 
since p,=Cpp’ ie, by 


— (Pe (4) 


We use a prime to denote partial differentiation 
with respect to 7. The subscript 7 indicates that 
the quantity is to be evaluated at r=r;; the 
subscript 0 will be used for r=ro. Since an 
electron moving on the instantaneous circle has 
a momentum p;, we see, using (3), that (4) is 
another expression of the relation (1). 

The condition for motion in the equilibrium 
circle (r;=79, a constant independent of time) 
is then 

pe=9, ’ (5) 


since py ro is independent of time, and we take 
A of the form? A(r, t)=f(r)h(t). 

It can readily be shown that the condition 
A’, =0 is equivalent to (2). We have 


fH-de=S curlA-do= fA-ds=2mrAo, 


and from (3), /7>=Ao/ro. Thus on the equilibrium 
circle, ¢9 = 

It may be remarked that for a magnetic field 
everywhere of the form thus 
lacking the central field necessary to satisfy 
the flux condition (2)) the vector potential 
is A=b/(2—n)r""', and (4) gives er?/I;/c 
=(2—n)p,/(1—n), or, from (1), pr;=constant. 


* The case in which there is a time independent flux ¢; 
through the center of the orbit is also of interest. It can 
be represented by A(r, t) =f(r)h(t)+ 1/227, which necessi- 
tates only replacing py, in (5), by py+ed:/2ze. 
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In such a field the electron spirals rapidly 
towards the center as its momentum increases. 


Radial focusing 


In order to take advantage of the adiabatic 
theorem, we write r=r;+x, and suppose v<r,. 
If we introduce x and p,=p, as new canonical 
variables, and in the new Hamiltonian function, 
Ki =H retain only terms not higher than 
quadratic in x and p, we obtain 


eps e+ pip" ix (6) 


Here e=c(n’? +p?) is the energy of an electron 
moving on the instantaneous circle. 
The equations of motion are 


i= Cp, ri, 


x «6 (7) 
which, by eliminating p,, give 
er etc'pip” ix €. (8) 


The forcing terms on the right represent the 
nonadiabatic corrections to the motion. 

If we write x=(mc* )'u, and omit the forcing 
terms for the moment, (8) becomes 


titw-u=0, 


since the second and third terms of w,* are 
smaller than the first by a factor of order 
(chm €)°f? where f is the frequency of the 
magnetic field, v the frequency of revolution of 
the electron, and p,, is the maximum momentum 
attained by the electron. If we write 


rr, (10) 
Eq. (9) takes the simple form 
Nv. 


The frequency of radial oscillation is thus just 
n; times the frequency of rotation of the electron 
in its orbit. In general n; will be a slowly varying 
function of 7, but for /7~1/r", n; is independent 
of r, n;=(1—n)!. We also see that for focusing 
oscillations to occur we must have p”;>0, 
which, on the equilibrium circle, reduces to 
A”, >0, in virtue of (5). The equilibrium circle 
is thus the point of minimum electric field. 

If we choose t=0 when A =0, we require the 
solution of (9) only for w,f>>1, since this condi- 
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tion is already met at the time of injection of 
the electron. The asymptotic form of the solution 
for w,f>1 is 


u=a(c/row,)! snl f 
or 


0 


with a and y arbitrary constants. The amplitude 
of the radial oscillations is thus damped by a 


factor (p; r:)>?=(ell; 
injection of the electron 


An electron injected on its instantaneous 
circle, will, in the adiabatic limit, continue to 
move on this circle as it shrinks down to the 
equilibrium orbit. The rate of contraction of the 
instantaneous circle can be obtained by differ- 
entiating (4), 


r= —eA’; cp" A’; (A";+2A'; r;). 


kor r7;—7ry small this reduces, in virtue of (5), 
to Ayo The shift of the 
instantaneous circle toward the equilibrium orbit 
produced by a change 6// of the magnetic field 
is thus given by 


or; = — 611 71. (12) 


An electron injected not on the instantaneous 
circle, but in its neighborhood, will oscillate 
about the instantaneous circle; its motion is 
described by (11). Since the instantaneous circle 
is simultaneously approaching the equilibrium 
orbit, the resultant motion is that shown in 
Fig. 1(C). 

Both the damping of the oscillation about the 
instantaneous circle and the shrinking of the 
instantaneous circle toward the equilibrium orbit 
are used to cause the electrons to miss the 
injector on successive revolutions. In terms of 
the @eltage gain per revolution 61’, and the 

Ceei voltage I’, the damping is expressed by 
Whbtetsa,a is the fractional decrease in the 
ump) de of oscillation. Similarly the shrinking 
by (12) is 


(14) 


We may verify that the wandering from the 
instantaneous circle caused by the nonadiabatic 
correction terms in (8) is in fact very small. 
The term é7;'¢€ is of order #,v*/c?, and will be 
neglected. Let t; be the time of injection. The 
asymptotic solution of (8) for w,.>1 can be 
written down with the aid of the solutions (11). 
It is 


x= snl f wt 
ty n,(t) p(t)ri(s)} w,(s) 


(ri(tr) 


xo f wat |- 
w,(t)* 


The maximum displacement is 


Xm = —2F 
In the interesting case, near the beginning of the 
acceleration when the field is increasing linearly 
with time, differentiation of (12) gives #; 
= H. Vf 6H and 6r; are the increments 
of /7 and r; during one revolution of the electron, 
we then have 


x, = 6r 611, 
all quantities being evaluated at t=¢;. 


Vertical focusing 


Near the median plane the magnetic field has 
a radial component determined by 0//, dz 
=011, dr. The vector potential then takes the 
form 

A,=Al(r, t)-—3Ailr, D2, 
where 
For small deviations from the median plane, 


the Hamiltonian (6) must be supplemented by 
the terms 

The equations of motion are 
p:=CpiL(rp)’ €, 
which differ from (7) only in the absence of a 
forcing term, and in the replacement of p”; by 
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The solutions are thus the same 


Ji. 


as (11), but with w, and n; replaced by wa. 
=my,r; and m,, defined by 
ry re. (15) 
Choice of the field shape 
It has been shown that there must be a 


position of minimum electric field (i.e., 
A”)>0) for radial stability, and that the field 
must decrease with increasing distance from the 
axis for axial stability of the electrons. If the 
field in the region of the electron’s orbit is of 
the form 

IT, vr)", 


the vector potential is given by 
ro" 


2—atr! r 
The second term in the bracket gives an electric 
field, but no magnetic field at the position of 


the orbit; it represents the effect of the strong 


MEITNER 


central magnetic field required to satisfy the 
flux condition (2). 

Evaluating (10) and (15) for this potential 
we find 


n° =(1—n), 


The conditions for radial and vertical focusing, 
n?>0, mZ>0, thus demand 0<n<1, as we 
have already seen. 

An electron deflected from its instantaneous 
circle by a small angle scattering will execute 
oscillations of amplitude proportional to the 
transverse velocity imparted by the scattering 
divided by the initial frequency of oscillation. 
The ratio of radial to vertical amplitude will 
thus be 


s=[n (1—n) }}. 


Since in the Illinois accelerator the vertical 
clearance was smaller than the horizontal, it was 
desired to have vertical focusing stronger than 
horizontal focusing. The pole pieces were 


accordingly shaped to give a field with n= 3. 
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The Resonance Energy of the Thorium Capture Process 
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Forskningsinstitutet for Fysik, Stockholm, Sweden 


N a paper dealing with the capture cross 
section for thermal neutrons in thorium! it 
was pointed out that the pure capture process in 
thorium has a resonance character with a large 
contribution from thermal neutrons.” The follow- 
ing experiments were carried out in order to 
determine the resonance energy of this process. 
As the neutron source available was only 100 
mg Ra+Be, it was not possible to determine the 
resonance energy in the usual way, i.e., by 
measuring the absorption coefficient of boron for 
the resonance neutrons with thorium 233 as 
detector. We therefore made use of a less direct 
method by measuring the absorption in thorium 


11. Meitner, Nature 45, 422 (1940). 
2L. Meitner, O. Hahn and F. Strassmann, Zeits. f. 
Physik 109, 538 (1938). 
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for various neutron resonance groups picked out 
by suitable detectors. If there are found neutron 
groups having absorption coefficients whose ratio 
is not inversely proportional to the ratio of the 
respective velocities, one can expect that a 
resonance group of the thorium capture process is 
lying within the range of the energies under 
investigation. 

Since the large contribution of thermal. neu- 
trons to the capture process in Th suggested that 
the resonance energy might be rather low, 


Au (E,=3.5 ev),* In (E,=0.9 ev)! and Rh, which 


has the same resonance energy as In were used as 


8Q. R. Frisch, Kgl. Danske Vid. Selsk. Math. “Phys, 
Medd. 14, Nr. 5. 

he =" Baker and R. F. Bacher, Phys. Rev. 57, 1076. 
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detectors. The resonance energy of these three 
elements is rather well-determined and_ their 
resonance levels are known to be single. In 
addition, use could be made of the absorption 
coefficient of Th for thermal neutrons known 
from our earlier experiments. As already men- 
tioned, there can be no doubt that the resonance 
energy of Th is beyond the energy of the neutrons 
strongly absorbed by Cd, hence the 1, v law will 
be valid for the absorption of thermal neutrons. 

The experimental arrangement was generally 
as follows: The neutron source was placed inside 
a block of paraffin wax 3.4 cm below the surface. 
The detectors were placed 1.4 cm above the 
surface and shielded in all experiments by 0.44 
mm Cd on both faces. The face opposite to the 
source was screened by a thick foil of the material 
under investigation. Between the detectors and 
the paraffin surface the thorium absorber (13.4 
g cm? of metallic Th), or other absorbers, could 
be inserted. The contribution to the activity of 
the detectors from neutrons beyond the resonance 
region was determined by carrying out the 
irradiation with absorbers of the same material 
and sufficiently thick to cut off practically all the 
resonance neutrons. One set of experiments was 
carried out inside paraffin with the same arrange- 
ment as had been applied for the determination 
of the capture cross section of thermal neutrons 
in Th.! When account is taken of the slightly 
changed geometrical conditions, the values found 
were the same as those obtained with the 
arrangement described above. All activity meas- 
urements were followed over several half-life 
periods in order to get a higher accuracy, and 
every experiment was repeated at least three times. 

Because of the arrangement used-— the detector 
being placed directly on the absorber of nearly 
equal size—one can expect that only a small part 
of the scattering cross section enters into the 
measurements. If, however, the capture cross 
section is small itself, the effect of scattering 
might be very important. Since it was not pos- 
sible to determine this contribution for Th as 
absorber, the absorption of the resonance neu- 
trons of In in lead was measured under the same 
experimental conditions. The lead absorber had 
the same dimensions and, within 2 percent, the 
same number of absorbing nuclei per cm? as the 
thorium absorber. The capture cross section for 


thermal neutrons in lead is known® 


to be 
Pb 


= 2.5 X10-*! cm?’, and there are good reasons 
to suppose that the 1/v law is holding for lead. 
So the difference between the cross section found 
experimentally and that one computed according 
to the 1/v law can be regarded as the contribution 
from scattering within the whole range of 
energies under investigation. However, the scat- 
tering effect implies a change in the primary 
directional distribution of the neutrons. In order 
to evaluate this influence, the relative number of 
neutrons indicated behind an absorber of 23.8 
mg,cm? boron was measured with and without 
the lead absorber in position. With the Pb ab- 
sorber inserted, the apparent absorption coeffi- 
cient of boron was 10 percent larger than that one 
found without the lead. After applying the 
corresponding obliquity correction, the contri- 
bution of scattering in lead was found to be 
o,=1.2XK10-* 

From the results given below it can be seen 
that the scattering is important only for the 
absorption of the resonance neutrons of gold in 
Th. On account of the low absorption coefficient, 
these measurements are not very accurate, the 
difference in scattering to be expected for Pb 
and Th is negligible. The same holds for In and 
Rh, as the scattering effect is small compared to 
the capture process. So the absorption coefficients 
found experimentally for Th were corrected for 
the scattering contribution found with lead. 

In calculating the absolute values of capture 
cross sections, obliquity corrections (angles up 
to about 70° were involved) were made according 
to the curves for thin, compared to thick, de- 
tectors given by Frisch,’ without taking into 
account any influence of the change in the 
directional distribution caused by the scattering 
effect. We are interested only in the ratio of 
these cross sections to the cross section for 
thermal neutrons, which, as far as the latter is 
corrected in the same way, is independent of the 
applied correction. The results are given in 
Table I. The last column gives the cross sections 
calculated by assuming the validity of the 1/9 
law and by taking the capture cross section for 


Th 
thermal neutrons asa basis (oy, = 6.0K 10>*4 em®*). 
®R. Fleischmann and W. Bothe, Ergeb. d. 
Naturwiss. 16, 37 (1937). L. Meitner, reference 1. 
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TABLE I. Data on cross section values. 


In 0.9 ev 3.540.7 1.0 
Kh 0.9 ev 3.540.7 1.0 
~0.3 0.5 


Au 3.5 ev 


In the calculation we set Ey, = 0.026 ev. 

The results show that, for E=0.9 ev, the 
experimental cross section is much larger than 
the calculated cross section, while for E=3.5 ev 
the experimental cross section appears to be 
somewhat smaller than that calculated according 
to the 1,7 law. On account of the poor accuracy 
of the results it seems meaningless to calculate 
the exact resonance energy from the values ob- 
tained by assuming the existence of a_ single 
resonance level. Hence one can conclude only 


MELTNER 


that thorium has a resonance level between 0.2 ev 
(cut-off energy of Cd) and 3.5 ev, and probably in 
the neighborhood of 2 ev. 

We attempted to get additional information by 
carrying out some experiments with the resonance 
neutrons of silver (EZ, = 2.5 ev). The activities due 
to the resonance neutrons were, however, so low 
even without any absorber that it was not 
possible to obtain any reliable results. 

The absorption of the In resonance neutrons 
in thin lavers of boron was found to correspond 
to a resonance energy of 0.92 +0.06 ev, if for the 
boron absorption coefficient of thermal neutrons 


the value Pv =28 cm? g is taken. 

I wish to thank the Academy of Sciences in 
Stockholm for a grant and Professor G. Borelius 
and Doctor J. Tandberg for the indium and 
rhodium foils kindly put at my disposal. 
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Letters to the Editor 


ROMPT publication of brief reports of important dis- 

coveries in physics may be secured hy addressing them 
to this department. Closing dates for this department are, for 
the first issue of the month, the eighteenth of the preceding 
month, for the second issue, the third of the month. Because 
of the late closing dates for the section no proof can be shown 
to authors. The Board of Editors does not hold itself responsible 


for the opinions expressed hy the correspondents. 


Communications should not in general exceed 600 words 
in length. 


On a Theory of Particles with Half-Integral Spin 
WILLIAM RARITA AND JULIAN SCHWINGER 
Department of Physics, University of California, Berkeley, California 
June 18, 1941 


IERZ and Pauli! have developed a general theory of 
F particles with arbitrary spin, both integral and half- 
integral. The quantities describing particles of integral 
spin are tensors of suitable rank; half-integral spins are 
described by spinors of multiple order. It is the purpose 
of this note to suggest an alternative formulation of the 
theory on half-integral spins which avoids the complicated 
spinor formalism of Fierz and Pauli. The fundamental 
quantities of this theory Wur...0e, have the mixed trans- 
formation properties of a Dirac four-component wave 
function, and of a symmetric tensor of rank k. (The Dirac 
index is suppressed.) The equations, which describe a 
particle of spin k+ } are: 


(yrOr +x) = 0, = 0. ( 1 ) 


The usual supplementary conditions of the integral spin 
theory: 
=(), Waaus...uc =, (2) 


appear as consequences of these equations. 

The verification that such wave fields are indeed asso- 
ciated with a particle of spin k+} may proceed either by 
demonstrating that the number of independent plane wave 
solutions of definite energy and momentum is 2(+ })+1, 
or by a direct proof that the square of the intrinsic angular 
momentum has the value (k+ 3)(&+ 3) in the rest system. 
(It will not have this value in an arbitrary reference sys- 
tem.) The spin multiplicity is an invariant and is easily 
calculated in the rest system. Since each W obeys the 
Dirac equation, the ‘‘small’’ components vanish and the 
two “‘large’’ components altogether form 2(&+3)!/k!3! 
quantities. The second equation of (1) provides the infor- 
mation that the W4ue...0e=0, which are 2(k+2)!/(k—1)!3! 
in number, and that o;W jus...ue=0 (us #4) which constitute 
2(k+1)!/(k—1)!2! relations among the W’s. The total 
number of independent components is, therefore, 2(k+1), 
as desired. The operator of total spin consists of the sum 
of the & infinitesimal rotation operators associated with the 
tensor indices in addition to the spin operator }@ of the 
Dirac theory. The proof that the spin possesses the correct 


PHYSICAL 


61 


REVIEW VOLUME 60 
eigenvalue is an elementary consequence of the condition 
ips... =(), 

The special case of spin } has been treated in detail by 
Fierz and Pauli employing a method that necessitated the 
adjunction of auxiliary spinors in order to produce the 
proper subsidiary conditions from a variational principle. 
It is an advantage of the new formalism that a suitable 
Lagrangian can be constructed without the intervention 
of additional fields. One of a class of possible Lagrangians is 


+ —k)yyy. (3) 


Although the Lagrangian is not unique in the absence of 
external fields, the form (3) is the only one which permits 
a relatively simple expression of the equations of motion 
in the presence of electromagnetic fields. The four-vector 
of charge and current, obtained by the usual prescription, is 


+ 

= (4) 
in the absence of external electromagnetic fields. Although 
the charge density is not a positive definite form, the total 
charge is necessarily definite since ¥,=0 in the rest system. 
In the exceptional case of zero rest mass the wave function 
admits a gauge transformation, 

which leaves all physical quantities invariant. 

The method here presented for developing the theory of 
spin } thus contains many of the features of both the 
Proca and the Dirac theory. This is particularly evident 
in the application to B-decay which Mr. Kusaka® has dis- 
cussed in an accompanying letter. 


= 0, 


'M. Fierz and W. Pauli, Proc. Roy. Soc. A173, 211 (1939). 
2S. Kusaka, Phys. Rev. 60, 61 (1941), this issue. 


3-Decay with Neutrino of Spin } 
SHUTCHT KUSAKA 
Department of Physics, University of California, Berkeley, Calisornia 
June 13, 1941 


LL current theories of B-decay assume the neutrino 

to have the spin } although there is no a priori reason 
for this choice, except that of simplicity. Recently Oppen- 
heimer! has suggested the possibilities of neutrinos with 
different masses and spins, and it is of some interest to 
find what distribution laws, energy-lifetime relations, and 
selection rules for 8-decay are possible with a neutrino of 
spin 3. The results show that for any coupling involving 
the neutrino wave function, or its first derivative, the 
spectrum distribution is predominantly of the Konopinski- 
Uhlenbeck type, and that, correspondingly, the energy- 
lifetime relation for high energy is given by the seventh- 
power law. Further, for neutrino of non-zero mass, both 
the Fermi and the Gamow-Teller selection rules are 
possible; while for neutrino of zero mass, only the latter 
is permitted. Now, since the energy-lifetime relations for 
the observed 8-spectra agree better with the fifth-power 
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law given by the Fermi distribution rather than the 
seventh-power law, we have experimental grounds for 
ruling out the possibility of neutrino with spin 3. 

The wave equation for a particle of spin } has been 
developed by Rarita and Schwinger.? The wave function 
¥, Consists of a four-vector, each component of which is a 
Dirac spinor, and satisfies the supplementary conditions 
ty=9 and y%,=0. The plane wave solutions can be 
built up from the polarization vectors 


1 1 


where the z axis is taken in the direction of the momentum, 
and the solutions ¥,; and ¥_, of the Dirac equation for the 
spin parallel and antiparallel to the momentum. If we 
separate the space and time components of ¥, by writing 
® = and ig=y4, the four independent solutions 
for positive energy are 


=ey 
1 
H — @3y-1, = Vi p 
Vv 
1 
evitvi—evi, 


They are orthogonal, normalized, and satisfy the supple- 
mentary conditions. 

Following the usual recipe, the interaction between the 
heavy and light particles can be written down. Because of 
the supplementary conditions, there is no scalar coupling 
term. The vector coupling term is 


which gives in the non-relativistic limit for the heavy 
particle 

The tensor and the pseudo-vector coupling terms can 
similarly be written down. 

Since the plane wave solutions have terms with the 
factor E/u or p/u, the square of the matrix elements con- 
tains a factor of the form E+ay where a is a numerical 
constant depending on the type of coupling. Hence, for 
high energy, this theory gives the seventh-power law for 
the relation between the lifetime and the energy. It has 
actually been verified, by calculating the transition proba- 
bility for the most general type of interaction, that there 
is no possible choice of the constants which would lead 
to a cancellation of the energy and momentum terms and 
so give the Fermi law. 

When the mass of the neutrino is taken to be zero, the 
wave function admits a gauge transformation. Now it is 
easily seen that any interaction involving the neutrino 
wave function is not gauge invariant: in fact the only 
interactions satisfying this condition are those involving 
the curl of the wave function. Thus only the tensor and 
the pseudo-vector couplings are possible, and hence only 
the Gamow-Teller selection rule is permitted. However 
the energy-lifetime relation is still given by the seventh- 
power law since there remains only the first two of the 
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plane wave solutions given above,’ and the differentiation 
introduces the momentum factor. 

In conclusion, the author wishes to thank Professor 
J. R. Oppenheimer and Dr. Julian Schwinger for helpful 
discussions on this subject. 


' J. R. Oppenheimer, Phys. Rev. 59, 908 (1941). 

2W. Rarita and J. Schwinger, Phys. Rev. 60, 61 (1941), this issue. 
The author is greatly indebted to Drs. Rarita and Schwinger for letting 
him make use of their theory before publication. 

§ This is a special case of the general theorem due to M. Fierz, Helv. 
Phys. Acta 13, 45 (1940), that for a particle of zero mass and any spin, 
there are only two independent, non-trivial solutions. 


Observations Bearing on Stellar Collapse 


Dean B. MCLAUGHLIN 
The Observatory, University of Michigan, Ann Arbor, Michigan 
May 19, 1941 


HE recent publication, by G. Gamow and M. Schoen- 

berg,'! of an ingenious theory of the cause of nova 

outbursts makes it of interest to examine the observational 
facts concerning pre-nova and post-nova stars. 

About pre-nova stars, we know only that their photo- 
graphic luminosities are comparable with that of the sun, 
and that one and only one of them had its spectrum re- 
corded feebly before outburst. The faint record of Nova 
Aquilae 1918 was classified by Miss Cannon on the basis 
of violet extent of the spectrum, and not from spectral 
lines, as “resembling classes B and A.’ Actually class O 
would also be a possibility, since the violet extents of 
O and B spectra look alike when photographed with glass 
prisms. 

Some years after outburst, all stars which have been 
studied (sixteen in all) show continuous spectra with great 
violet extent and no absorption lines. Some have bright 
lines of hydrogen, Heu, Cit, and other elements which 
indicate temperatures similar to those of class O stars.® 
Of special importance is the fact that the photographic 
luminosities are equal to those of the pre-nova objects. 
No observational fact clearly contradicts the working 
hypothesis that the physical conditions of pre- and post- 
nova stars are nearly identical. The burden of proof thus 
rests with those who postulate that stars change from 
normal to degenerate configurations at the time of outburst. 

The collapse theory encounters grave difficulties in the 
attempt to reconcile radical changes of radius and tem- 
perature with the observed equality of pre- and post-nova 
brightness in a narrow range of wave-lengths, the common 
photographic region. 

To fix ideas, suppose a normal main sequence star 
collapses to a sub-dwarf or white dwarf configuration such 
that the bolometric luminosity is unchanged. Then, assum- 
ing blackbody radiation, we must have 


where Ro and 7» are the original radius and temperature, 
and R, and 7, refer to the collapsed star. If we take 
To=6000°K (equal to that of the sun) and 7.=50,000°K 
(estimated from spectra of old novae), we find R.=0.014Ro, 
and the density would obviously be increased by a factor 
of 3X10. These figures are characteristic of the known 
white dwarfs. Of more direct application here is the result 
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that the photographic luminosity of the collapsed star 
should be roughly three magnitudes fainter than its pre- 
nova value. Put differently, since such a large bolometric 
reduction results from the change of temperature, we 
conclude that the observed equality of pre- and post-nova 
stars requires that the end product be three magnitudes 
brighter bolometrically than before outburst. To the 
writer, at least, this appears to be a reductio ad absurdum. 
More details are given in a paper now in press in Popular 
Astronomy. 

It is, of course, possible to object that the assumption 
of equal bolometric luminosities is incorrect. This is granted 
as a likely possibility. It is, therefore, of interest to deter- 
mine in which direction it may be expected to fail. Appeal- 
ing to observation, we see that while main sequence stars 
in general fit the mass-luminosity relation, the white dwarf 
companion of Sirius is too faint for its mass. This is just 
opposite to the direction of failure which would be required 
to account for the equality of the pre- and post-nova 
photographic luminosities on the hypothesis of collapse. 

We must, of course, admit the existence of a remote 
possibility that a star could collapse, with radical changes 
of all its physical parameters except mass, and vary the 
energy distribution in its spectrum so as to keep the photo- 
graphic radiation the same. This would imply almost a 
diabolical intent to deceive the observer! And, as one of 
my colleagues has remarked, it would take a pretty ‘‘cagey”’ 
star to do the trick. 

The above remarks apply only to the common novae, 
and are not meant to include the supernovae. There are 
reasons for believing that the latter originate in a different 
way from the former, and because the energy liberated 
in these greater explosions is of the order of magnitude 
to be expected from stellar collapse, it appears quite 
possible that the theory of neutrino emission would apply 
there. The pre-outburst state of supernovae is entirely 
unknown, and the preceding arguments, therefore, cannot 
be applied to them. 

'G. Gamow and M. Schoenberg, Phys. Rev. 59, 539 (1941). 


? Harvard Annals 81, 179 (1920). 
3M. Humason, Astrophys. J. 88, 228 (1938). 


Classical and Quantum Reflections of X-Rays 


C. V. RAMAN AND P. NILAKANTAN 
Department of Physics, Indian Institute of Science, Bangalore, India 
May 7, 1941 


N a series of papers! published in April and May, 1940, 

we presented evidence proving conclusively that the 
lattice planes in a crystal give a second kind of geometric 
reflection of x-rays which we designated as the modified 
or quantum reflection, to distinguish from the unmodi- 
fied or classical reflection discovered by Laue and his 
collaborators in 1912. The process which results in the 
modified reflection was clearly established by the experi- 
mental results in the specific cases studied by us. In the 
language of classical optics, a modified reflection by the 
lattice planes results from the dynamic variation of their 
structure amplitude consequent on an oscillation, relative 
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to each other, of the interpenetrating lattices in the struc- 
ture of the crystal. In the act of such reflection, the primary 
x-ray frequency is altered by the addition or subtraction 
of one or another of the characteristic infra-red frequencies 
of the crystal. In the language of quantum mechanics, the 
modified reflection represents an inelastic collision of the 
photon with the crystal in which the two exchange energy. 
The probability of such transition, and therefore also the 
intensity of reflection, is determined by the principles of 
quantum mechanics and not by classical dynamics. The 
theory of the effect has been discussed in further papers.? 

The view has been put forward in some recent papers** 
that the phenomenon dealt with by us may be explained 
as “diffuse maxima in the scattering of x-rays by elastic 
waves of thermal origin.”” We desire energetically to con- 
tradict this suggestion. The theory of the thermal scatter- 
ing of x-rays has been fully discussed by Laue® in 1926. 
It is therefore possible to anticipate the results to be 
expected on the basis of such scattering. It may be stated 
at once that they are quite incapable of explaining the 
phenomena actually observed by us. 

In the first place, the intensity of the scattering of x-rays 
by elastic waves of thermal origin comes out proportional 
to N, the number of scattering atoms in the crystal, and 
not to N®. The absolute intensity of scattering is therefore 
very small, and the effects due to it would be perceptible 
only when large volumes of the crystal are operative. The 
present effect, on the other hand, though decidedly feebler 
than the classical reflections, can be recorded under pre- 
cisely the same experimental conditions as the latter, that 
is to say, with fine x-ray beams and the thinnest crystal 
plates. 

The so-called maxima in scattering arise from the pres- 
ence of a factor in the expression for its intensity propor- 
tional to the square of the wave-length of the clastic waves 
which are effective. However, it is only when the direction 
of scattering coincides exactly with that of the classical 
reflection, and is therefore unobservable, that this multi- 
plying factor becomes important. For other settings of the 
crystal, the “‘peak’’ disappears and is replaced by broad 
humps of relatively small intensity in the curve of intensity 
of scattering. These cannot explain the discrete reflections 
actually observed over a wide range of settings of the 
crystal. The (111) reflection by a crystal of diamond, for 
example, is remarkably sharp and may be readily photo- 
graphed with crystal settings, in which it appears dis- 
placed by as much as 10° on either side of the Laue spot. 
No significant changes are observable either in its size or 
in its definition over the whole of this range. 

Thirdly, exact measurements of the position of the (111) 
reflection by diamond over a wide range of settings have 
been made. They fit perfectly into the Raman-Nath for- 
mula, the inclination of the phase waves to the lattice 
planes coming out as half the tetrahedral valence angle. 
On the other hand, considered over the whole of this 
range, the Faxén formula is not even a rough approxima- 
tion to the truth. 

Finally, a crucial test is provided by the low tempera- 
ture results with the (111) reflections by diamond. As 
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Laue has shown, a change of frequency occurs also in the 
scattering by elastic waves. \ccording to the “scattering” 
theory, the intensity maxima correspond to elastic waves 
of maximum wave-length and minimum frequency. In such 
a case, therefore, the transition probabilities should be 
given correctly by the classical dynamics. The intensity 
maxima should fall off in proportion to the absolute tem- 
perature, the more exactly, the nearer the “scattering” 
maximum is to the Laue spot. \ctual observations show, 
however, that the intensity is practically unaffected by 
cooling the crystal down to liquid-air temperatures, and 
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this is true irrespective of the setting of the crystal. The 
observations thus clearly show that the transitions involved 
correspond to the high frequency infra-red levels and not 
to the low frequency acoustic ones. They also prove that 
the modified reflection is a quantum-mechanical effect. 


'C. Ve. Raman and P. Nilakantan, Current Science 9, 165 (1940); 
Proc. Ind. Acad. Sci. 11, 379, 389, 398 (1940). 

7C. V. Raman and N.S. N. Nath, Proce. Ind. Acad. Sci. 12, 83, 427 
(1940); Raman, Proc. Ind. Acad. Sci. 13, 1 (1941). 

IW. H. Zachariasen, Phys. Rev. 59, 207 (1941). 

1S. Siegel, Phys. Rev. 59, 371 (1941). 

ott. A. Jahn and K. Lonsdale, Nature 147, 88 (1941). 

®M. v. Laue, Ann. d. Physik 81, 877 (1926). 
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Proceedings of the Ohio Section of the American Physical Society 


HE fourth meeting of the Ohio Section of 
the American Physical Society was held at 
Muskingum College, New Concord, Ohio, on 
March 29, 1941. About fifty-five members and 
guests were present at the luncheon. Fourteen 
papers were scheduled for presentation, though 
weather conditions curtailed the programme; 
abstracts of five are appended. 

The fifth meeting of the Ohio Section was held 
at Cleveland, in the buildings of the Case School 
of Applied Science and the Western Reserve 
University, on May 10, 1941; it was a joint 
meeting with the Ohio Academy of Science. One 
hundred members and guests were present. 
Nineteen papers were scheduled for presentation ; 
of these the abstracts of five are appended. It was 


voted to cooperate with the Ohio Academy in 
sponsoring the programme of the Junior Academy 
for the creation of further interest in science in 
secondary schools, and five members (J. Albright, 
F. L. Berger, H. J. Kersten, H. P. Knauss, H. H. 
Roseberry) were chosen to collaborate in the 
various parts of the state. A symposium on 
“Light and Lighting” is planned for the meeting 
intended to be held in Cleveland in October 1941. 
Newly: elected officers of the Section for the 
season 1941-42 are the following: Chairman, 
W. E. Forsythe; Vice Chairman, P. B. Taylor; 
Secretary-Treasurer, R. H. Howe. 


RicHarD H. Howe, Secretary 
Denison University, Granville, Ohio 


ABSTRACTS OF PAPERS PRESENTED AT THE MUSKINGUM MEETING 


1. The Vibrational Spectrum of the Potassium Chloride 
Crystal Lattice. L. L. Fotpy, Case School of Applied 
Science. (Introduced by R. S. Shankland.)—The methods 
of Born, v. Karman, and Blackman for the investigation 
of the thermal vibrations of crystal lattices are applied in 
detail to an atomic model representative of a potassium 
chloride crystal. The potential interaction between ions, 
assumed of equal mass, is represented as a Coulomb inter- 
action for any ion and its eighty nearest neighbors, a van 
der Waals attractive interaction varying as the inverse 
sixth power of the distance between nearest neighbors, and 
an exponential repulsive interaction also between nearest 
neighbors. Explicit solutions of the secular equation for 
the frequencies, which is of the third order, are obtained 
for waves traveling in certain directions in the crystal. 
The calculated principal restrahlen frequency was found 
to be about 10 percent lower than the observed value. 
The dispersion for short longitudinal waves traveling 


normal to the 111 planes of the crystal was found to be 
so great that the frequency increases with wave-length 
for wave-lengths less than about four and one-half times 
the distance between nearest neighbors. This extreme be- 
havior was not manifested for transverse waves traveling 
normal to the 111 planes or for any waves traveling normal 
to the 110 and 100 planes. 


2. The Harmonic Analysis of Geiger Counter Pulses. 
Ropert S. SHANKLAND AND Ricuarp H. Case 
School of Applied Science.—The discharge of a Geiger- 
Miller counter results in a sudden drop in the potential of 
the wire electrode, followed by a more gradual return of 
this potential to its normal value. The time variation of 
potential constitutes a voltage pulse of which the shape 
and duration depend on both the discharge mechanism in 
the counter itself and upon the values of the electrical 
elements that compose the counter circuit. Thus the voltage 
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pulses observed on a cathode-ray tube may be analyzed 
to give information about the nature of the counter dis- 
charge and the behavior of the circuit. The present method 
of analysis gives the Fourier integral components making 
up the pulse in the form of a frequency spectrum of the 
amplitudes and phases. This is accomplished by charting 
the pulse to several scales of abscissae and analyzing the 
resulting curves by means of the Henrici harmonic analyzer 
according to a method previously described.! The pulses 
analyzed include an idealized pulse and those character- 
istic of counters used with the Neher-Harper circuit, the 
Neher-Pickering circuit, and a circuit employing a series 
resistance. 
'R.S. Shankland, J. Acous. Soc. Am. 12, 383-386 (1941). 


3. Heat Transfer Across Horizontal Air Spaces. Louis 
B. YounG, Marietta College—Heat transfer across hori- 
zontal air spaces is encountered in problems involving the 
calculation of heat transmission through floors and ceilings. 
Much of the transfer results from radiation exchange 
between the boundary surfaces; the remainder is due to 
convection. Convection in the case of upward heat flow, 
as in winter, is much greater than in the case of downward 
heat flow, as in summer. In both cases the actual transfer 
exceeds the possibility of ideal conduction as a transfer 
mechanism. Test apparatus was a two-foot guarded hot 
plate. Heat flow was calculated from the temperature 
gradient through a cork slab of known conductivity which 
was placed “‘in series” with the air space. The data enable 
computation of heat transmission through various floor 


and ceiling constructions incorporating different methods 
of insulation. The transfer is the sum of convection, follow- 
ing the five-fourths power law, and radiation exchange, 
following a modified Stefan’s law. 


4. Some Observations on an A.C. Resonant Pendulum. 
Davip K. WeimerR.—A condenser and an air core induc- 
tance in series are tuned to resonance at 60 cycles. .\ 
pendulum bob of iron suspended so as to swing into the 
coil is thereby kept moving at its own frequency. Measure- 
ments were made on the current and voltage as the bob 
swung in and out of the coil in an effort to determine the 
nature of the driving force acting on the pendulum. It was 
found that the current is greater as the bob entered the 
coil than it is in the corresponding position as the bob leaves 
the coil for the greater part of the cycle of the pendulum. 


5. Motion of the Walls of a Cornet. HAro_p P. KNAvss, 
Ohio State University. —Vibrations of the metal walls of a 
cornet were studied by using a crystal pick-up in contact 
with the metal. The cathode-ray oscillograph showed wave 
forms which are complex for the lower tones of the cornet 
and sinusoidal for the higher tones. By substituting for the 
sweep voltage the output of a microphone responding to 
air-borne vibrations, a kind of Lissajous figure was ob- 
tained, which was quite irregular for the lower tones and 
elliptical for the higher tones. For a range of about half an 
octave it was shaped like a figure 8, indicating that in this 
range the metal vibrates at half the frequency of the air 
column. 


ABSTRACTS OF PAPERS PRESENTED AT THE CLEVELAND MEETING 


1. Observation of “Forbidden” Bands in the Infra-Red. 
C. H. Tixpa., Ohio State University—According to the 
simple theory of tetrahedrally symmetric molecules there 
are four different normal frequencies designated as ¥;, v2, vs, 
and v4. The first two of these are optically inactive and are 
non-degenerate and twofold degenerate, respectively. The 
latter two are optically active, are threefold degenerate, 
and should appear as bands with sharp lines forming P, Q, 
and R branches. While this is essentially true for v3 it is not 
so for v4. The complex structure in v4 has been explained by 
Jahn! as due to a resonance Coriolis interaction between v4 
and inactive v2, which are only about 100 cm™ apart. 
Because of this interaction, the wave functions for »2 and 
vs become sufficiently “mixed” so that v2, becomes active in 
the infra-red. The region from 94 to 134 has been resolved 
for SiH, and GeHy, so that rotational lines are clearly 
distinguishable in two overlapping bands. The two bands 
are identified as v2 and v4, the Q branches folding in oppo- 
site directions as predicted by the theory of Jahn. The 
weaker band is identified as v2 and the other as v4. This 
identification is further verified since the rotational line 
spacing satisfies the relation presented by Teller? 

Av3+Avg= 1} Are. 
'H. A. Jahn, Proc. Roy. Soc. A168, 469, 495 (1938). 


A Tellier, Hand- und Jahrbuch der Chemischen Physik, Vol. 9, No. 2, 
p. 152. 


2. Absorption Bands of GeH, and SiH, in the Near 
Infra-Red. JosePpH W. StRALEY, Ohio State University. 
—The fundamental vibration-rotation bands Js, Js, and J, 
of GeH, and SiH, have been carefully measured and the 
moments of inertia and the X—Y distances calculated. 
Measurements of combination bands J; +J4 and in 
SiH, have also been made. The discrepancy in the appear- 
ance of the satellite structure in GeH, reported by Steward 
and Nielsen! with that reported by Lee and Sutherland? can 
now be attributed to difference in the absorbing layers of 
gas used by these observers. The recent measurements 
indicate the intensity of the satellites is not more than 10 
percent of that of the principal rotation lines. The moments 
of inertia in g-em? in the various states are as follows— 

To Te Ts 


Silla 0.986 0.974 X10 0,993 0.993 
Gels 1.030 1.021 1.032 1.039 


IW. B. Steward and Hi. H. Nielsen, Phys. Rev. 48, 861 (1935). 
tT and G. B. B. M. Sutherland, Proc. Camb. Phil. Soc. 35, 341 
(1939). 


3. Studies of Metallic Friction. k. 1. SrrouGH AND W. E. 
Rupp, Case School of Applied Science.—For the purpose of 
studying metallic friction under varying conditions of speed 
and normal load, a rotary testing apparatus was con- 
structed. In this apparatus a flat disk-shaped sample of one 
material to be tested was mounted on a faceplate ang 
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rotated in contact with three small pads of the other 
material. These pads were held in place by another 
stationary faceplate arranged so that the torque acting on 
it through friction between the samples could be measured. 
From this the coefficient of friction was computed. Friction 
phenomena of copper sliding dry on steel were studied with 
the above apparatus. Both kinetic and static coefficients of 
friction were observed to be higher for ground (crystalline) 
surfaces than for polished (non-crystalline) surfaces. Fric- 
tion was observed to increase uniformly with the number of 
rotations of the disk for both types of surfaces until a 
critical point was reached where galling occurred and the 
coefficient of friction rose rapidly to a high value. Polished 
surfaces withstood about three times as many rotations 
as ground surfaces under similar conditions before galling 
occurred, The coefficient of friction was observed to be 
independent of sliding speed from speeds of 20 cm per 
second to the limit of the apparatus (200 cm per second) 
throughout the experiments. 


4. Charge and Mass of Electronic Particles. ALFRED 
LANbDE, Ohio State University—The problem of explaining 
the empirical formula e =ch/137 is attacked from various 
angles that gradually seem to yield a reasonable picture of 
how the electronic charge is determined by quantum theory 
and relativity. From the classical point of view an electron 
can be characterized (1) by its rest energy Eo proportional 
to its mass m, (2) by its ‘“‘radius,”’ or, better, by a charac- 
teristic time period ¢) inversely proportional to the mass. 
The customary radiation theory yields the wrong result 
that the energy is infinitely large and the radius is infinitely 
small. In order to get finite results a change of the theory is 
indispensable. We have suggested a formal change of the 
energy expression, analogous to the energy reduction due to 
radiative damping of a vibrating electron. The resulting 
finite energy Eo is connected with a finite “damping radius” 
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ro and with a finite time period to = 3mc/4re® characteristic 
for the damping. From the guantum point of view the 
product Koto can have only certain selected values that are 
obtained by a quantization in the energy-momentum space 
of Einstein’s equation and in corre- 
sponding space ?—(r/c)?=te. The smallest eigenvalue of 
Eofo turns out to be ~h/210. This together with Eo= me? 
and with the value of fo vields e? =hc/140. A correction is 
necessary because fo was calculated for infinitely weak 
fields as against the actually finite self-field. 


5. Determination of Mercury Arc Temperatures by 
X-Rays. Cart W. J. Karasu, General Electric 
Co. and Case School of Applied Science—Molybdenum 
radiation at 32,000 volts, filtered through zirconium 
dioxide was used. Data for a graph of transmission vs. 
moles of Hg per sq. cm in the path of the beam were 
obtained by measuring transmissions of tubes with known 
amounts of Hg vapor. Direct-current temperature data 
were obtained for a hollow electrode, quartz lamp of 20-mm 
diameter and 29.0-cm arc length. Hg pressures, appendix 
eontrolled, were measured with a quartz pressure gauge." 
Axial transmissions were measured. End corrections were 
determined by thermocouple measurements of wall temper- 
atures, and pyrometer measurements of electrode tempera- 
tures. Axial temperatures were calculated from the perfect 
gas law. For a constant pressure of 1.0 atmosphere, the 
axial temperature varied approximately as the 20th root of 
the input, from 5950 deg. K at 30.7 watts per cm to 6100 
deg. K at 58.5 watts per cm. For a constant input of 30.7 
watts per cm, the axial temperature increased from 5950 
deg. K at 0.22 atmosphere to 6450 deg. K at 0.40 atmos- 
phere and then decreased to 5950 deg. K at 1.0 atmosphere. 
The probable error of the temperatures is estimated as 
150 deg. K. 

Carl Kenty, Rev. Sci. Inst. 11, 377 (1940). 


